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The  interaction  of  the  pressure  fluctuations  beneath  a  turbulent  boundary 
layer  (TBL)  and  a  structure  cause  that  structure  to  vibrate.  These  vibrations 
may  cause  unwanted  effects  such  as  radiated  noise  in  jet  aircraft,  spacecreift, 
naval  vehicles  and  air  conditioning  systems;  or  self  noise  in  transducer  systems. 
The  TBL  pressure  fluctuations  can  usuedly  be  described  as  rcindom  stationary 
processes  in  space  and  time.  It  is  generally  useful  to  characterize  the  spatial 
and  temporal  variations  in  terms  of  their  Fourier  trcinsforms,  wavenumber  and 
frequency.  The  vibrations  of  a  structure  excited  by  these  pressure  fields  also 
are  random  stationary  processes.  It  is  the  wavenumber-frequency  dependence  of 
these  structural  vibrations  that  is  of  interest  in  this  study. 

The  wavenumber-frequency  response  of  the  forced  transverse  vibrations  of 
a  one-dimensional  structxire,  represented  by  a  beam,  is  under  consideration  for 
this  study.  The  excitation  consists  of  discrete  shalcers  with  measured  force  levels 
instead  of  random  pressure  fluctuations.  The  velocity  is  sampled  along  the  length 
of  the  beam.  Both  the  force  imd  velocity  spatial  distributions  are  transformed 
into  wavenumber-frequency  space  via  the  Fourier  transform.  As  the  frequency 
response  of  structures  is  generally  well  understood  and  easy  to  meeisure,  the 
emphzisis  of  this  study  is  placed  on  understanding  the  wavenumber  response. 

The  general  objective  of  this  study  is  to  characterize  the  wavenumber  Altering 
characteristics  of  a  finite  beam  over  a  large  region  of  wavenumber-frequency 
space.  The  study  is  both  theoretical  and  experimental.  Results  from  theoretical 
models  on  the  bezim  response,  characteristics  of  the  shakers  used  to  excite  the 
beam  and  algorithms  used  in  taking  Fourier  transforms  over  spatially  finite 
domains  are  used  to  assist  in  the  design  of  the  experiments  and  analysis  of  the 


mezisurement  results. 


Just  as  an  impulsive  load  in  time  htis  a  white  or  flat  frequency  spectrum, 
so  does  a  spatially  impulsive  load  (i.e.  a  point  source).  Initially  it  seems,  a 
good  method  would  be  to  excite  the  beam  with  a  single  shtdcer  and  measure  the 
wavenumber  velocity  response  with  an  array  of  accelerometers.  As  the  input  force 
wavenumber  spectrum  would  be  flat,  all  velocity  wavenumbers  would  be  excited. 
However,  there  are  several  difficulties  with  this  idea.  One  is  the  fact  that  the 
spatial  location  of  the  shaker  is  important  to  the  velocity,  not  just  the  relative 
distance  between  the  excitation  point  and  the  measurement  sensor.  Hence,  even 
at  a  single  frequency,  the  output  wavenumber  velocity  spectrum  will  depend  on 
the  location  of  the  drive  force.  Two,  the  beam  is  a  fairly  selective  wavenumber 
filter.  At  any  frequency  there  is  a  preferred  wavenumber,  the  freebending 
wavenumber,  at  which  the  velocity  response  of  the  beam  will  peak.  This  is  the 
wavenumber  that  corresponds  to  the  wavelength  of  a  freely  propagating  wave 
on  an  infinite  beam  at  that  frequency.  Measurement  of  the  beam  response  at 
this  wavenumber  is  easy  because  the  signal  is  very  high.  With  only  one  point 
force  exciting  the  structure,  the  freebending  wavenumber  would  dominate  the 
response  making  measurement  at  other  wavenumbers  difficult.  In  an  attempt  to 
excite  the  beam  to  measurable  levels  at  wavenumbers  other  than  the  freebending 
wavenximber,  a  forcing  fimction  weighted  to  a  particular  wavenumber  other 
than  the  freebending  wavenumber  could  be  used.  This  requires  multiple  point 
drives  to  be  used.  In  this  way,  the  energy  of  force  excitation  as  well  as  the 
velocity  response  would  be  concentrated  in  a  region  of  wavenumber-frequency 
space  centered  around  the  drive  wavenumber.  The  velocity  response  would  be 
measured  by  the  accelerometer  EU'ray,  which  provides  a  further  concentration  of 
the  measurement  in  the  desired  wavenumber  frequency  region. 
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The  philosophy  of  driving  the  beam  at  a  peu’ticulzu:  wavenumber  is  simil^  to 
steady  state  frequency  analysis.  The  difference  is  that  in  the  temporal  domain  a 
signal  of  indefinite  length  can  be  used  while  spatially  the  signal  length  is  limited 
to  the  length  of  the  structure.  Indeed,  the  meaning  of  Fourier  transforms  over 
the  finite  length  of  the  structure  is  different  than  that  for  the  infinite  domain. 

What  is  proposed  is  the  development  of  a  system  to  measure  the  wavenumber- 
frequency  characteristics  of  a  beam  by  using  miiltiple  point  drives  and  measuring 
the  velocity  response  with  an  array  of  accelerometers.  Hence,  it  is  necessary  to 
imderstand  the  theoretical  aspects  and  heirdware  rezdization  of  the  mezisurement 
system.  Predictions  of  the  response  of  the  beam  are  derived  so  that  experimental 
results  on  both  point  admittances  and  wavenumber  spectra  may  be  validated. 
The  signal  processing  necessary  to  measure  the  waveniimber  spectra  of  the  force 
and  the  velocity  is  developed.  The  signal  processing  aspects  of  imposing  the 
desired  spatial  force  distribution  on  the  beam  is  also  derived.  The  transducers 
for  exciting  the  beam  and  measuring  the  exciting  forces  are  developed  as  well 
as  an  eqviivalent  circuit  model  to  predict  the  response  of  the  force  transducer. 
Finally,  the  equipment  and  procedure  used  for  data  acquisition  are  developed. 

1.1.  Literature  Review 

The  early  work  on  the  wavenumber  response  of  structures  seems  to  be 
motivated  by  the  need  to  imderstand  the  vibrations  induced  by  jet  and  rocket 
TBL  pressure  fluctuations  [1,2].  Some  of  the  later  work  is  more  motivated  toward 
the  measurement  of  flow  induced  TBL  due  to  vehicles  moving  through  a  fluid 
[3,4].  Essentizdly,  each  of  these  papers  studies  the  relation  between  the  applied 
pressure  field  and  the  resulting  vibration  or  radiated  pressure  due  to  the  induced 
vibration. 

Starting  with  Uberoi  and  Kovztsznay  [5],  the  wavenumber  sensitivity  of 
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sensors  with  finite  extent  has  been  studied.  Maidanik  and  Jorgensen  [6]  developed 
the  wavenumber  sensitivity  of  arrays  of  flush  moimted  pressure  transducers. 
Blake  and  Chase  [7]  examine  arrays  of  pressure  sensors  with  more  complicated 
wavenumber  pressure  sensitivity.  In  this  thesis,  the  sensitivity  of  arrays  of 
accelerometers  is  predicted  using  the  emedysis  of  Blake  and  Chase  [7]  as  well 
as  of  Gaskill  [8]. 

Hutto  [9]  investigated  the  use  of  multiple  point  drives  weighted  to  a 
wavenumber  to  excite  a  bezim.  His  theoretical  study  into  this  method  forms  some 
of  the  basis  of  this  reses^ch.  Roberts  and  Sabot  [10]  employ  a  simil2U'  technique 
utilizing  multiple  point  drives.  However,  they  attempt  to  simulate  the  TBL 
pressure  fluctuations  on  a  plate.  Hutto’s  work  on  the  response  of  Timoshenko 
beams  to  point  drives  is  extended  in  this  thesis.  The  multiple  shaker  system  is 
constructed  and  experiments  are  compared  to  theoretical  predictions  based  on 
Timoshenko  beam  theory.  A  system  for  the  measurement  of  the  wavenumber- 
frequency  pressme  response  of  structures  is  discussed  by  Strawderman  [11]. 

1.2.  Experimental  Approach 

The  basic  experimental  approach  employed  in  the  remainder  of  this  thesis 
for  measuring  the  wavenumber-frequency  response  of  a  beam  is  outlined  in  this 
section.  Also  discussed  are  the  anedytic  and  experimental  requirements  imposed 
by  the  proposed  system  which  forms  much  of  the  basis  for  this  thesis. 

1.2.1.  Experimental  Method 

The  experimented  method  consists  of  driving  the  structure  at  a  particular 
wavenumber  and  measuring  the  velocity  response  as  a  function  of  frequency. 
To  achieve  the  force  distribution,  multiple  drivers  are  placed  along  the  length 
of  the  beam.  The  voltage  input  to  the  shakers  is  weighted  to  a  particular 
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wavelength  by  setting  the  voltage  inputs  to  the  shakers  to  cos(Jtji„,),  where 
kd  is  the  drive  wavenumber  and  Xm  is  the  location  of  a  shaker.  The  resulting 
velocity  distribution  is  measured  by  an  accelerometer  placed  at  discrete  points 
along  the  beam  using  a  common  reference.  The  Fourier  transforms  of  the  force 
and  velocity  axe  performed  to  obtain  the  wavenumber  response  of  both  quantities. 
It  is  the  relationship  between  these  wavenumber  dependencies  that  characterizes 
the  spatial  filtering  action  of  the  beeim  or  the  admittance  of  the  finite  beam  in 
wavenumber-frequency  space. 

A  typical  force  amplitude  distribution  weighted  to  a  pau’ticular  wavenumber 
is  shown  in  Figure  1.1.  By  using  multiple  drives  in  a  stzinding  wave  pattern, 
the  peak  of  the  forced  wavenumber  content  is  steered  to  the  drive  wavenumber, 
kd,  without  the  need  for  continuous  control  of  the  phaise  of  each  of  the  drive 
voltages,  as  would  be  required  by  a  traveling  wave  pattern.  The  mezisurement  of 
the  velocity  distribution  by  the  array  of  accelerometers  also  allows  a  filtering  of 
the  wavenumber  dependence.  The  tise  of  the  discrete  Fourier  transform  applied 
to  the  velocity  data,  steers  the  measurement  wavenumber  to  the  different  bins  of 
the  discrete  transform  (Chapter  2).  The  hardware  and  analytic  implications  of 
this  method  are  described  next. 

1.2.2.  System  Requirements 

Clezu-ly,  one  of  the  most  importemt  requirements  is  the  need  for  multiple 
shakers  (force  exciters).  Because  the  force  is  to  be  measured  at  each  excitation 
point,  each  shaker  must  have  an  integral  force  gage.  Measurement  of  the  force 
is  required  for  two  main  reasons.  First,  the  force  measurement  allows  for  the 
determination  of  whether  or  not  the  force  wavenumber  spectrum  tracks  with  the 
applied  drive  voltage  wavenumber  spectrum.  A  feedback-control  system  could 
be  used  to  control  the  force  at  each  of  the  shakers,  but  this  is  a  much  more 
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f 


I  ' 


.4 

Non~d I man . 


.6 

L«n9th  s  X/L  s  ( 


Discrete  forcing  function  weighted  to  a  wavenumber 
(after  Hutto  [9]). 
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complicated  implementation.  Second,  the  forces  at  each  of  the  shakers  are  the 
input  to  the  analytic  models  which  predict  the  velocity  along  the  beam.  To 
prescribe  the  desired  drive  voltage  distribution,  a  multiple  channel  amplitude 
weighting  system  is  essential.  This  amplitude  weighting  system  must  have  the 
ability  to  drive  each  of  the  shakers  at  an  arbitreiry  amplitude  (read,  positive  or 
negative  amplitude).  The  requirement  of  one  channel  per  shaker  is  cut  in  half  by 
using  a  symmetric  weighting  pattern  as  in  Figure  1.1.  The  weighting  distribution 
used  is  symmetric  about  the  center  and  equal  to  cos  kd{x  —  L/2),  where  L  is  the 
length  of  the  beam.  The  velocity  of  the  beam  at  an  array  of  locations  must  be 
measured. 

Both  the  force  and  velocity  measurements  must  be  transformed  into  wave- 
number  space  via  Fourier  transforms.  Because  the  velocity  may  only  be  measured 
at  discrete  points,  the  continuous  Fourier  transform  is  simiilated  by  the  discrete 
Fourier  transform.  Because  of  the  fimte  length  of  the  beam,  the  velocity  and 
force  only  exist  over  a  finite  length. 

Validation  of  the  measurements  is  required.  The  force  gages  must  be 
calibrated  and  resonant  effects  in  the  force  gages  identified.  This  requires  zmailytic 
modeling  of  the  shzdcer/force  gage  units  used  in  the  measurements.  Also,  restilts 
firom  the  measured  response  of  the  beam  must  be  validated.  Analytic  modeling 
of  the  beam  response  is  employed  to  validate  the  measurement  results. 

1.3.  Outline  of  Thesis 

In  Chapter  2,  the  means  by  which  the  wavenumber-frequency  spectra  of 
the  force  and  velocity  eue  measured  is  described.  The  entire  experimental 
measurement  system  is  described.  The  basis  for  the  signal  processing  assump¬ 
tions  is  presented.  A  description  of  the  shaker/force  gage  unit  is  given.  In 
Chapter  3,  the  analytic  model  for  Timoshenko  beams  is  derived  and  discussed. 
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First  the  spatial-frequency  response  of  free-free  beams  to  multiple  point  drives 
is  derived  following  Hutto  [9].  The  wavenumber-frequency  response  of  free-free 
finite  beams  and  infinite  becims  is  developed.  Finally,  the  relationship  between 
the  wavenumber  content  of  the  forcing  function  and  the  velocity  response  is 
examined.  In  Chapter  4  the  combination  force  gage/shaher  units  are  modeled 
using  equivalent  circuit  analysis.  The  results  of  the  experiments  axe  compared 
to  the  various  theoretical  results  in  Chapter  5.  Lastly,  the  conclusions  and 
recommendations  for  future  research  cire  drawn  in  Chapter  6. 


Chapter  2 

DESCRIPTION  OF  MEASUREMENTS 

2.1.  Introduction 

In  this  chapter,  the  implementation  of  the  experimental  method  proposed 
in  Chapter  1  is  described.  First,  a  general  outline  of  the  experiments  is  given; 
including  descriptions  of  the  experiments,  meEisurement  system  and  procedure. 
The  design,  construction  and  calibration  of  the  shaker/force  gage  units  follow 
the  experimental  description.  Next,  the  signal  processing  used  both  to  process 
the  force  and  velocity  data  as  well  as  calculate  the  amplitudes  of  the  voltages 
sent  to  the  shakers  is  presented.  The  analysis  in  this  final  section  spells  out  the 
various  assumptions  made  in  the  signal  processing  aspects  of  handling  the  data. 

2.2.  Experimental  Description 

In  this  section,  the  overall  description  of  the  experiments  is  presented.  First, 
the  beam  geometry  and  material  properties  are  listed.  The  forcing  fimctions 
that  define  the  nine  experiments  performed  are  described  next.  The  system  and 
procedure  by  which  the  force  and  velocity  data  eu:e  obtained  are  detailed.  Also, 
the  means  of  force  excitation  is  given.  Finally,  general  aspects  of  performing  the 
experiments  are  enumerated. 

2.2.1.  Beam  Geometry 

The  steel  beam  used  in  the  experiments  is  37.9  cm  long  by  1.27  cm  wide 
by  1.92  cm  high  (or  thick).  The  material  properties  of  the  beam  are  listed  in 
Table  2.1.  The  length  of  the  beam  was  determined  by  two  constraints.  The 
number  of  channels  available  for  driving  the  shakers  is  ten;  using  a  symmetric 
weighting  fimction  with  respect  to  the  center  of  the  beam  the  number  of  shakers 
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that  could  be  driven  is  nineteen  (with  the  center  shaiker  having  no  symmetric 
pair).  The  closest  center  to  center  shaker  spacing  that  is  feeisible  is  2  cm;  hence 
a  38  cm  beam  would  be  fully  covered  by  nineteen  shakers. 

Table  2.1  Steel  material  parameters  used  in  this  thesis. 


Material  Parameter 

Symbol 

Value  (units) 

Elastic  Modulus 

E 

210x10^*^  (dyne/cm^) 

Shear  Modulus 

G 

77x10^®  (dyne/cm^) 

Density 

P 

7.8  (grams/cm*) 

Damping  Loss  Factor 

n 

This  is  the  damping  loss  factor  of  shaker-beam  system,  not  just  of  the  steel  beam. 

The  beam  is  supported  on  two  highly  porous  foam  wedges,  typicedly  used 
in  air  anechoic  rooms,  Hutto  [9]  showed  that  this  method  of  support  closely 
approximates  a  free  boundary  condition.  One  of  the  drawbacks  to  the  use  of  the 
wedges  is  that  approximately  1  cm  is  required  to  support  the  beam  at  either  end. 
Hence,  for  the  2  cm  shaker  separation  distance,  only  seventeen  and  not  nineteen 
shakers  could  be  used.  This  is  because  there  is  not  enough  room  to  place  the  last 
two  shakers  on  the  beam  and  still  properly  support  the  beam. 

The  beam  was  scored  at  1  cm  intervals  to  aide  in  the  location  of  the  shakers 
and  the  accelerometer.  The  first  mark  is  at  0.95  cm  and  the  last  mark  is  at 
36.95  cm.  These  spacings  are  symmetric  with  respect  to  the  middle  of  the  beam. 
It  is  at  each  scored  location  that  the  acceleration  is  measured. 

2.2.2.  Experimental  Test  Matrix 

The  experimental  test  matrix  is  shown  on  Table  2.2.  The  values  tmder  the 
heading  “wavelengths  across  the  beam”  are  equal  to  the  drive  wavelength  (see 
Figure  1.1)  divided  by  the  length  of  the  beam.  The  drive  wavenumber  is  kd  = 
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2-k/X.  The  nine  shakers  were  located  4  cm  apart  emanating  from  the  center  of 
the  beam;  with  four  shakers  to  the  left  of  center  shaJcer  and  four  to  the  right  of 
the  center.  The  seventeen  shalcers  were  spaced  2  cm  apart,  also  symmetrically 
from  the  center.  The  procedure  by  which  the  voltage  amplitudes  are  determined 
is  detailed  later  in  this  chapter  in  the  section  on  signal  processing. 

Table  2.2  Experimental  test  matrix. 


Experiment 

No. 

Wavelengths  Across 
the  Beam 

Drive 

Wavenumber  (cm“^) 

Number  of 
Shakers 

1 

0 

0.000 

9 

2 

1 

0.166 

9 

3 

2 

0.332 

9 

4 

4.5 

0.746 

9 

5 

0 

0.000 

17 

6 

1 

0.166 

17 

7 

2 

0.332 

17 

8 

4.5 

0.746 

17 

9 

9 

1.492 

17 

2.2.3.  Measurement  Instrumentation  and  Procedures 

The  overall  plan  for  the  measiuement  scheme  is  outlined  in  Chapter  1.  Each 
experiment  is  characterized  by  the  number  of  shakers  used  to  excite  the  beam 
and  the  voltage  distribution  applied  to  the  shakers.  For  each  of  the  experiments, 
the  spatial  distribution  of  both  the  force  and  velocity  must  be  measmed  in  order 
to  describe  the  wavenumber-frequency  filtering  characteristics  of  the  besun.  A 
description  of  the  instrumentation  and  procedure  by  which  the  force  and  velocity 
are  measured,  as  well  as  the  means  of  exciting  the  beam,  are  foimd  in  this  section. 

The  shakers  are  glued  to  the  beam  using  epoxy  (Epoxy-Patch  309  gray  made 
by  the  Hysol  division  of  Dexter).  Before  attaching  the  shakers  with  this  fast 
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setting,  high  tensile  strength  epoxy,  all  surfaces  must  be  cleaned.  First,  any  rust 
is  removed  from  the  bccim  by  using  steel  wool.  Then,  both  the  attachment  point 
on  the  beam  and  the  shaker  headmass  eire  cleaned  using  a  three  peirt  process. 
First,  a  cotton  swab  is  immersed  in  toluene  and  rubbed  on  the  surface  to  be 
cleaned  until  no  further  residue  can  be  seen.  This  procedure  is  repeated  with 
alcohol  and  then  acetone.  After  the  surfaces  are  cleaned,  each  surface  is  coated 
with  epoxy  and  held  together  either  by  hand  or  by  a  vise  until  set.  In  Figure  2.1, 
a  photograph  of  seventeen  shakers  glued  to  the  beam  is  shown. 


Figure  2.1  Photograph  of  nine  shakers  glued  to  the  beam. 


A  flowchart  of  the  measurement  system  is  shown  in  Figure  2.2.  An  IBM 
PC- AT  controls  the  experiments  by  driving  an  HP  4102A  impedance  analyzer 
through  a  GP-IB  (General  Purpose-Interface  Bus).  The  impedance  analyzer 
sends  an  oscillator  signal  at  a  specified  voltage  and  frequency  through  a  low  pass 
filter  set  at  12  kHz  to  an  HP  465A  voltage  amplifier.  The  measurements  are  all 
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tciken  using  a  step  frequency  method.  Note  that  this  figure  is  shown  specifically 
performing  a  measurement  of  transfer  admittance.  For  the  measurement  of  the 
force  ratios,  the  accelerometer  signal  is  removed  and  a  force  signals  is  sent  to  the 
impedance  analyzer. 

The  amplified  and  filtered  signed  is  feed  into  the  multi-channel  amplitude 
weighting  system  (MAWS).  The  system  heis  the  ability  to  use  the  single  input 
from  the  voltage  amplifier  to  transmit  up  to  ten  separate  voltage  signeds.  These 
signals  will  have  arbitrary  gains  and  either  a  zero  or  ISO  degree  phase  difference 
from  the  oscillator  signal.  MAWS  allows  the  amplitude  of  the  drive  signal  of 
the  array  of  shakers  to  be  tuned  to  a  particular  wavenumber  by  dialing  in  the 
voltage  signals  for  each  of  the  shakers  to  preccdculated  values  using  the  procedure 
described  in  Section  2.4.  Because  MAWS  has  only  zero  or  180  degree  phase 
control,  only  a  real  spatial  voltage  distribution  may  be  imposed  unto  the  drive 
signals.  As  previously  stated,  there  are  only  ten  channels;  a  symmetric  force 
distribution  is  used  so  that  up  to  twenty  shakers  may  be  driven  by  this  ten 
channel  system. 

The  center  shaker  is  used  as  the  reference  shaker.  The  signal  from  each  of  the 
other  force  gages  is  sent  to  a  channel  selector,  which  is  simply  a  switch  that  can 
select  any  of  the  output  signals  from  the  force  gages.  The  force  output  voltages 
from  the  reference  zmd  measurement  force  gages  are  amplified  by  Ithaco  450 
amplifiers.  The  amplitude  and  phase  of  the  ratio  of  the  two  force  gage  signals 
are  measured  by  the  impedsmce  anadyzer.  Of  course,  this  measurement  is  taken 
while  all  of  the  other  shakers  are  being  driven.  The  results  are  stored  by  the  PC 
as  a  function  of  frequency.  This  procedure  is  repeated  until  the  ratios  of  all  of  the 
shakers  to  the  center  shaker  are  meaisured.  The  stored  data  are  later  transferred 
to  the  mainframe  computer. 
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The  velocity  distribution  is  measured  in  a  manner  similar  to  the  measurement 
of  force  ratios.  An  accelerometer  glued  to  a  samarium-cobedt  magnet  serves  as 
the  velocity  probe.  The  same  epoxy  used  to  glue  the  shakers  to  the  beam  is  used 
to  attach  the  magnet  to  the  accelerometer.  The  magnet  performs  very  well  in 
the  frequency  range  of  interest  [9].  The  beam  is  made  of  magnetic  steel  so  the 
velocity  is  szimpled  simply  by  moving  the  accelerometer  down  the  length  of  the 
beam.  It  was  found  that  the  response  of  the  accelerometer/magnet  combination 
is  enhanced  by  using  Kistler  force  gage  grease  between  the  beam  surface  and  the 
magnet.  At  each  scored  location  cdong  the  beeim  the  ratio  of  the  voltage  from  the 
accelerometer  to  the  center  force  gage  is  measured  using  the  impedance  ainadyzer. 
As  part  of  the  data  acquisition  program  this  ratio  is  converted  into  admittzince 
(velocity  divided  by  force)  and  the  admittance  as  a  function  of  frequency  and 
location  is  stored  on  the  PC. 

2.3.  Shaker/Force  Gage  Units 

The  requirements  for  the  hardware  of  the  measurement  system  were  laid  out 
in  Chapter  1.  These  requirements  are  placed  into  two  categories;  treinsducer  and 
data  acquisition  requirements.  The  data  acquisition  system  was  just  described. 
In  the  sections  on  transducer  development,  the  design  and  construction  of  the 
final  shaker  is  described  along  with  interim  design  steps.  A  brief  description  of 
construction  techniques  is  followed  by  the  method  for  the  calibration  of  the  force 
gages.  In  Chapter  4,  an  equivalent  circuit  model  describing  the  overall  response 
of  the  shaker,  and  the  response  of  the  force  gage  section  is  developed.  This  model 
serves  as  a  check  on  the  measurements  made  on  these  transducers. 

The  design  and  construction  of  the  combination  force  gage  euid  driver  (or 
shaker)  is  a  major  part  of  this  study.  Because  shedcers  of  a  diameter  of  less  than 
1.27  cm  required  in  this  study  are  not  commercially  available,  it  is  necessary  to 
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fabricate  them.  The  original  design  concept  of  the  drive  section  included  the  use 
of  PZT-4  (Lead  Zirconate  Titonate)  piezoceramic  disks  as  the  active  element. 
The  use  of  two  of  these  disks  and  a  p2issive  tailmass  were  added  to  the  shaker 
design  during  the  iterative  steps  of  the  design  process.  Also  in  the  original  design, 
the  piezo-film  PVDF  (also  known  as  PVF-2)  was  used  in  the  force  gage  section. 
This  design  was  abandoned  after  many  iterations  due  to  mech2miceLl  difficulties 
and  PZT-4  was  used  in  the  final  force  gage  design.  The  force  gage  design  using 
the  PVDF  film  is  shown  in  Appendix  A. 

2.3.1.  Shaker/Force  Gage  Design 

A  schematic  of  the  final  design  is  shown  in  Figure  2.3.  The  three  sections 
of  the  shaker  are  labeled  on  the  figure.  The  sections  axe  the  drive  section,  the 
decoupler  section  and  the  force  gage  section.  The  dimensions  of  the  various 
components  and  materials  used  are  given  in  Table  2.3. 

Table  2.3  Dimensions  and  materials  used  in  the  transducer. 


Description  of  Section 

Material 

Dimensions  (cm) 

Drive 

PZT-4 

.635  long  X  1.27  Diam. 

Tailmass 

Brass 

.635  long  X  1.27  Diam. 

Decoupler 

Aluminium 

1.905  long  X  1.27  Diam. 

Force  gage 

PZT-4 

.3175  thick  x  1.27  O.D.  and  .762 

Headmass 

Aluminium 

1.59  long  X  1.27  Diam.’*^'^ 

I.D.  and  O.D.  stand  for  inner  and  outer  diameter  respectively. 

The  headmass  section  narrows  to  an  area  of  .3175  cm  x  1.27  cm  at  the  top  portion. 


The  drive  section  consists  of  two  PZT-4  disks,  three  corrugated  nickel 
electrodes,  €ind  a  brass  tailmass.  The  PZT-4  elements  are  polarized  toward  the 
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SECTION 


DRWE 

SECTION 


aluminum  head  mass 


nickel  electrode 


- P2T 

_ _ -  nickel  ELECTRODE 


•CERAMIC  INSULATOR 


Figure  2.3 


Final  shaker/force  gage  design. 
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center  electrode  and  are  electrically  in  parallel.  The  negative  terminals  (ground) 
axe  the  two  outward  facing  electrodes;  this  configuration  aides  in  shielding  the 
force  gage  signal  from  the  drive  signal.  The  contaunination  of  the  force  gage 
signal  from  the  drive  voltage  is  called  cross-talk.  Proper  grounding  techniques 
aure  discussed  later  in  this  section.  The  PZT-4  disks  are  33  coupled,  which  means 
that  the  direction  of  the  applied  electric  field  is  the  same  as  the  resulting  motion 
(or  force).  Since  both  disks  are  polarized  toward  the  center  electrode  and  are 
electrically  in  parallel,  an  applied  electric  field  causes  the  disks  to  beat  against 
each  other.  In  other  words,  both  disks  contract  toward  the  center  electrode  and 
expand  away  in  unison.  At  low  frequencies,  the  action  of  the  two  piezoceramic 
disks  is  the  same  as  a  single  PZT-4  disk  with  twice  the  length  of  the  original 
disk  [12|.  The  motivation  for  this  design  and  the  use  of  the  brass  tailmass,  is 
to  increase  the  overall  compliance  and  mass  of  the  transducer  drive  section  so 
as  to  lower  the  resonant  frequency  of  the  system.  Thus,  the  low  frequency  force 
output  of  the  shaker  is  increased. 

The  negative  terminals  are  connected  using  an  insulated  wire  to  eliminate 
the  possibility  of  shorting  out  the  positive  terminal  since  the  positive  terminal 
is  between  the  two  negative  termineds.  The  wire  is  soldered  to  the  electrodes 
and  then  to  the  black  lead  of  a  twisted  wire  with  a  shield.  The  positive  lead  of 
the  twisted  pair  (the  red  wire)  is  soldered  to  the  center  electrode  of  the  shzdcer. 
Both  wires  are  then  glued  to  the  side  of  the  shaker  with  epoxy  to  strain  relieve 
the  connection.  A  General  Radio  or  Pomona  type,  two  prong  connector  is  used 
at  the  other  end  of  the  twisted  pair  wire  with  the  outer  shield  connected  to  the 
ground  terminsd  at  this  connector. 

The  decoupler  section  consists  of  an  aluminum  cylinder.  There  are  two  main 
reasons  for  the  decoupler.  The  first  reason  is  an  extension  of  the  mechanical 
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design  of  the  drive  section,  that  is  to  lower  overall  stiffness  of  the  shaker.  So 
as  not  to  block  the  output  of  the  drive  section  at  higher  frequencies,  a  materizd 
lighter  than  brass  is  used.  By  increasing  the  compliance  of  the  shaker,  the  use 
of  the  aluminum  decoupler  increases  the  low  frequency  axied  force  output  of  the 
shaker.  Second,  the  use  of  the  decoupler  also  reduces  the  effect  of  electronic 
cross-talk.  Since  the  drive  section  acts  eis  an  antenna,  it  is  beneficied  to  place 
the  force  cell  as  fax  from  the  drive  section  eis  possible  to  reduce  the  effects  of  the 
electronic  cross-tEdk.  The  aluminum  decoupler  must  be  connected  to  ground  in 
order  to  act  as  a  shield  for  the  force  gage. 

The  force  gage  section  is  designed  to  measure  the  applied  force  to  the  test 
structure  and  to  transmit  a  line  force  across  the  width  of  the  beam.  The 
force  gage  section  consists  of  an  annular  cylinder  33  coupled  PZT-4  piece,  two 
electrodes,  eui  alsmEug  ceramic  insulating  disk  (cerEunic  insulator  in  Figure  2.3), 
and  an  aluminum  headmass.  The  positive  electrode  of  the  force  gage  is  insulated 
from  the  decoupler  by  the  alsmsdg  disk.  The  ground  side  of  the  force  gage  is 
connected  to  the  aluminum  headmsiss.  The  electrodes  are  soldered  to  low  noise 
coaxial  cable.  The  center  conductor  of  this  cable  is  the  positive  or  high  side  and 
the  outer  conductor  is  the  shield.  At  the  end  of  the  coaxial  wire  a  Micro- Tech 
connector  is  used  to  help  reduce  imwanted  cross-talk.  The  cable  is  stress  relieved 
in  a  similar  manner  to  the  drive  section;  by  gluing  the  wire  to  the  side  of  the  force 
gage.  The  Eduminum  headmass  allows  for  light  stiff  connection  of  the  shaker  to 
the  test  structure.  At  the  attachment  point  the  cross-section  of  the  headmass  is 
nEu-rowed  to  simulate  a  line  force.  The  headmass  design  enables  the  shaker  to 
easily  Euid  safely  be  attached  and  removed  from  the  test  structure.  The  effect  of 
the  headmass  on  the  response  of  the  force  gage  is  examined  in  Chapter  4. 

The  proper  grounding  technique  is  a  critical  element  of  the  shaker  design  eis 
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well  as  for  the  entire  measurement  system.  The  shakers  must  be  designed  so  that 
the  force  gage  signal  is  insulated  from  the  drive  signal.  The  ground  plane  of  the 
drive  section  and  the  force  gage  section  must  be  the  same  so  that  no  ground  loop 
exists.  All  metal  elements  of  the  shaker  are  grounded  to  prevent  charge  from 
accumulating  on  these  components  and  becoming  an  emitter  of  electro-magnetic 
noise.  The  ground  planes  of  all  of  the  equipment  must  be  the  same.  This  can 
be  accomplished  by  running  all  of  the  eunplifiers,  oscilloscopes  and  any  other 
equipment  connected  to  the  shaker  off  of  the  seime  power  line. 

2.3.2.  Construction 

The  transducer  is  constructed  by  gluing  the  surfaces  of  the  components 
together  using  epoxy.  The  surfaces  must  be  meticulously  cleaned  before  the  glue 
is  applied.  After  the  construction  is  complete,  the  strength  of  the  glue  joints  is 
greater  than  the  ceramic  in  bending.  These  procedures  were  arrived  upon  by 
experimentation  with  different  methods  as  well  as  through  discussions  with  the 
personnel  at  the  Applied  Research  Lab  (ARL)  at  Penn  State  [13]. 

The  first  step  is  to  clean  the  assembled  components.  MetsJlic  surfaces  that 
aie  to  be  glued  (electrodes,  tailmass  and  decoupler)  are  first  etched  in  a  mild  acid 
bath  to  remove  oxide  coatings.  The  mettiUic  surfaces,  the  electroded  siufaces  of 
the  PZT  and  the  alsmeiig  insulating  disk  are  next  abraded  with  a  fiberglass  brush 
for  further  cleaning.  Then  all  of  the  surfaces  to  be  glued  undergo  the  three  part 
solvent  cleaning  process  of  toluene,  alcohol,  and  acetone;  this  process  is  performed 
under  an  exhaust  hood. 

After  cleaning,  the  components  zure  glued  together  using  Shell  EA-6  epoxy. 
First,  the  glue  is  applied  in  a  thin  coat  to  both  surfaces  that  are  to  be  joined. 
The  entire  transducer  is  zissembled  sequentially  from  the  tailmass  up.  Upon 
completion  of  the  gluing,  the  transducer  is  placed  in  an  alignment  jig  which  also 
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supplies  axiaJ  compression  to  insure  a  good  bond.  The  transducer,  now  in  the 
jig,  is  placed  in  an  oven  at  SO  degrees  Celsius  for  six  to  eight  hours.  The  jig  and 
transducer  are  removed  from  the  oven.  The  residual  charge  across  the  electrode 
terminals  from  the  heating  and  cooling  of  the  piezoceramic  is  shorted  out.  All  of 
this  work  is  performed  in  the  transducer  deem  room  at  ARL. 

Next,  the  electrodes  are  bent  slightly  toward  the  body  of  the  shaker  to  prevent 
the  electrodes  from  protruding.  Wires  are  then  soldered  to  the  terminals  of  the 
PZT-4  elements  and  then  stress  relieved.  Finally,  the  outer  paurt  of  the  transducer 
is  covered  by  shrink  wrap  so  that  the  electrodes  eire  not  exposed  and  are  protected 
during  handling. 

2.3.4.  Force  Gage  Calibration 

The  method  of  calibration  is  simply  a  comparison  calibration  which  employs 
Newton’s  law.  The  shaker  is  loaded  with  a  known  mass.  The  acceleration  of  the 
mass  is  measured  with  an  accelerometer.  From  the  ratio  of  the  force  gage  voltage 
to  the  accelerometer  voltage,  the  sensitivity  of  the  force  gage  is  determined.  A 
diagram  of  the  calibration  set-up  is  shown  in  Figtire  2.4 

A  personal  computer  controls  the  calibration  by  sending  instructions  to  an 
HP  4192A  impedmce  analyzer  using  a  GP-IB.  The  impedance  analyzer  has  the 
ability  to  send  ein  oscillator  signal  at  a  prescribed  voltage  and  frequency  as  well 
as  measuring  the  ratio  of  two  signals  and  their  relative  phase.  A  data  acquisition 
computer  program  was  written  to  run  the  calibration  and  store  the  data.  The 
program  steps  the  oscillator  through  the  frequency  range  of  interest  (2-12  kHz), 
while  the  oscillator  signal  from  the  HP  4192A  is  sent  to  an  HP  465A  voltage 
amplifier.  The  amplified  drive  signal  excites  the  shaker  which,  in  turn,  excites 
the  known  calibration  mass.  The  signed  from  the  force  gage  and  the  attached 
accelerometer  are  eunplified  emd  then  sent  to  the  impedance  analyzer.  The 
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Figure  2.4  Calibration  system. 
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amplitude  and  phase  of  this  ratio  sxe  used  to  calculate  the  unknown  sensitivity 
of  the  force  gage.  The  sensitivity  data  are  sent  to  the  VAX  11/780  maiinframe 
for  plotting. 

The  sensitivity  is  calculated  from  the  following  formula, 


S(u;) 


R(u)Sa 

Me  +  Mh  m  +  Ma 


(2.1) 


where,  S{u3)  is  the  frequency  dependent  force  sensitivity,  Sa  is  the  sensitivity 
of  the  accelerometer,  Mhm  is  the  mtiss  of  the  head  mass.  Me  is  the  calibration 
mass,  Ma  is  the  mass  of  the  accelerometer,  and  R{u)  is  the  ratio  of  the  voltage 
from  the  force  gage  to  the  accelerometer  voltage.  This  formula  comes  from  the 
application  of  Newton’s  Law. 

The  maximum  and  minimum  sensitivity  for  five  different  shakers  is  shown  in 
Figure  2.5  over  the  frequency  range  of  2-12  kHz.  Since  the  value  of  the  force  gage 
sensitivity  is  very  flat  over  the  frequency  range  of  interest,  a  single  sensitivity  is 
assigned.  This  sensitivity  is  the  average  of  the  frequency  dependent  sensitivity 
over  the  frequency  range.  The  averaged  values  of  the  sensitivity  euid  phase  for 
all  of  the  force  gages  is  tabulated  in  Table  2.4.  Each  of  the  shakers  has  a  number 
associated  to  it;  these  are  the  numbers  used  to  identify  the  shedeers  and  their 
sensitivities.  A  discussion  of  the  variations  of  the  force  gage  sensitivity  is  found 


in  Chapter  4. 


2.4.  Signal  Processing  Aspects 


In  this  section,  the  signzd  processing  aspects  of  measuring  the  wavenumber 
response  of  the  beam  are  examined.  First,  the  discrete  Fourier  transform 
is  introduced.  Next,  the  wavenumber  sensitivity  of  accelerometer  arrays  is 
developed.  This  development  begins  with  a  general  treatment  of  the  problem,  so 
that  the  simplifications  made  when  using  the  discrete  Fourier  transform  can  be 
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identified  and  justified.  Also,  the  measurement  of  the  wavenumber  content  of  the 
forcing  function  is  developed  and  a  method  of  steering  the  forced  waveniimber 
content  is  developed. 

2.4.1.  Discrete  Fourier  Transform  of  the  Velocity 


The  Fourier  transform  of  /  is  defined  to  be 

OO 

F(K)=  J  ,  (2.2) 

— OO 

and  the  inverse  transform  is 

OO 

/(x)  =  i  /  F(k.)e‘'‘-’dk,  ,  (2.3) 

—  OO 

where  the  wavenumber  is  and  the  spatial  variable  is  x. 

The  discrete  Fourier  transform  (DFT)  is  a  means  of  estimating  the  wave- 
number  spectrum  from  measured  and  predicted  spatial  data.  The  DFT  may  be 
applied  to  results  of  spatial-frequency  model  or  to  experimented  data.  The  DFT 
of  the  displacement  j/(i)  is 

N-i 

Yi{k)  =  Ax  y(£Ax)e(-^'^'*)  .  (2.4) 

r=o 

The  factor  k  is  the  wavenumber  counting  index,  £  is  the  sample  number  end 
N  is  the  total  number  of  points.  The  factor  Ax  is  the  spacing  between  sample 
locations.  The  multiplying  Ax  is  the  proper  normalization  for  the  DFT  [14].  The 
transform  of  a  complex  sequence  of  N  points  results  in  N  separate  spectral  lines, 
the  first  iV/2  spectral  lines  represent  the  positive  wavenumber  spectrum  and  the 
mirror  image  of  the  second  represent  the  negative  spectrum.  The  counting  index 


k  corresponds  to  the  wavenumber  fcr  in  the  following  mmner. 


fcr  = 


2nk 

2w(k-m 


for  =  0, 1 . . .  N/2  -  2,N/2  -  1 
for  k  =  N/2,N/2  +  l...,N 


(2.5) 


In  the  discretization,  both  the  spectrum  and  the  origineil  sequence  are  considered 
periodic.  Zero  padding  of  the  data  interpolates  between  points  in  wavenumber 
space  with  a  sin(i)/x  function  [15).  Zero  padding  better  defines  the  spectrum 
by  this  interpolation,  but  does  not  increase  the  resolution  capabilities  of  the 
given  trzinsform.  The  points  are  interpolated  between  the  existing  points,  but 
the  maximum  wavenumber  (tt/Ax)  remains  the  same. 

All  of  the  processing  of  experimental  velocity  data  utilizes  the  DFT.  In  the 
next  section,  the  assumptions  that  are  implicit  in  the  use  of  the  discrete  transform 
axe  discussed.  The  discrete  transform  is  used  as  an  approximation  to  the  exact 
transform;  the  exact  transform  is  performed  analyticedly  on  the  predicted  beam 
velocity  in  Chapter  3  and  used  as  a  means  of  assessing  the  impact  of  the  spatial 
sampling  required  in  the  measurements. 

For  the  processing  of  discrete  transforms  a  Fast  Fourier  Tremsform  (FFT) 
algorithm  is  used  by  zero  padding  [14,15]  to  a  power  of  two.  This  speeds  the 
calculation  of  the  transform,  which  is  important  even  though  the  transforms 
considered  in  this  thesis  are  not  large. 

2.4.2.  Wavenumber  Frequency  Sensitivity  of  Accelerometer  Arrays 

What  is  sought  here  is  the  sensitivity  of  an  array  of  sensors  (in  this  czise 
an  array  of  accelerometers)  in  wavenumber-frequency  space.  The  objective  of 
performing  this  analysis  is  to  understand  the  wavenumber  filtering  effect  of 
the  array  and  to  determine  the  spacing  of  the  sensors  needed  to  measure  the 
wavenumber  spectra.  Specific  results  in  this  section  are  presented  for  point  or 
ideal  sensors,  but  the  analysis  may  be  applied  to  a  more  general  claiss  of  sensors 

[7]. 

The  steirting  point  of  the  anedysis  is  to  derive  the  wavenumber  sensitivity  for 
a  single  sensor.  It  is  postulated  that  there  exists  a  large  class  of  sensors  whose 


voltage  output  due  to  some  excitation  C3in  be  caiste  in  the  form  of  a  convolution 
integral  [5,6].  In  other  words,  there  exists  an  impulse  response  function  for  the 
transducer.  Additionally,  the  impulse  response  function  is  talcen  to  be  separable 
in  space  and  in  time;  in  other  words  the  space-time  impulse  response  function 
may  be  written  as 

g(x,t)  =  h{x)f{t)  '1.6) 

In  Figure  2.6  the  geometry  of  the  single  sensor  is  shown  symbolically  measuring 
some  physical  phenomenon.  The  specific  sensor  under  consideration  here  is  an 
accelerometer  so  that  the  terms  accelerometer  and  sensor  are  used  synonymously 
in  this  section.  The  measurement  center  of  the  accelerometer  is  at  xq  units  from 
the  origin  with  some  finite  width  as  shown  in  Figure  2.6.  There  is  some  voltage 
output  from  the  sensor  which  is  proportioned  to  the  velocity  (this  is  buried  in 
f(i)  as  velocity  is  the  quantity  of  interest), 


Figure  2.6  Geometry  of  single  sensor  in  a  velocity  field. 
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The  form  of  the  convolution  integral  applied  to  obtain  the  measured  velocity 
by  a  sensor  located  at  Xq  is 


OO  OO 

Vm=  J  J  v{a,T)h{xo  -  a)f{t  -  T)dTdoL. 


(2.7) 


—  OO  — OO 


Since  the  steady  state  response  is  of  interest,  the  remaining  anzdysis  proceeds  in 
frequency  domain.  Tzdcing  the  Fourier  transform  with  respect  to  time  results  in 

OO 

vm  =  F{u)  /  u(q!,u;)/i(io  —  Q)da  ,  (2.8) 

—  OO 

where  the  convolution  theorem  has  been  applied  to  the  frequency  dependence 
[15].  For  brevity,  the  frequency  dependence  is  dropped  from  the  equations.  Using 
the  asterisk  notation  for  the  convolution  operation  on  i.  Equation  (2.8)  becomes 


Vm  =  Ki)  *  ft(x)]5(x  -  Xo),  (2.9) 

where  the  use  of  the  Dirac  delta  function  in  Equation  (2.9)  is  detailed  in  Bracewell 
and  Gaskill  [8,15].  This  use  of  the  delta  function  locates  the  sensor  at  xq.  In 
order  to  examine  the  wavenumber  filtering  effect  of  a  single  sensor,  the  spatial 
Fourier  transform  is  applied  to  the  right  heind  side  of  Equation  (2.9).  This  results 
in 

v„  =  (2.10) 

here  the  asterisk  notation  indicates  convolution  on  k^.  If  xo  is  zero,  then  the 
output  measured  from  the  sensor  is 

OO 

Vm=  J  V(k,)ff{k,)dk,.  (2.11) 

— OO 

If  the  sensor  is  a  point  sensor  then  H{kx)  =  1  because  h{x)  =  S(x)  and 
the  output  of  the  sensor  is  just  the  integral  over  all  wavenumbers.  For  h(x) 
not  equal  to  a  S(x),  then  H(kx)  represents  the  filtering  action  of  the  sensor. 
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For  ail  reaJ  sensors  h{x)  is  not  a  delta  function,  but  at  low  wavenumbers  the 
response  approximates  a  delta  function;  i.e.,  the  sensor  is  equally  sensitive  to 
all  wavenumbers  below  a  certsdn  wavenumber.  An  analysis  of  a  spatial  impulse 
response  function  that  is  not  a  delta  function  is  found  in  Blaice  and  Chase  [7]. 
Their  analysis  finds  that,  for  a  sensor  with  uniform  sensitivity  to  velocity  over 
the  contact  surface,  if  the  wavelengths  meet  the  following  criterion, 

A  >  4.2R  (2.12) 

where  R  is  the  radius  of  the  sensor,  the  sensitivity  will  not  be  attenuated  by 
spatial  averaging  by  more  than  3  dB.  This  assumes  that  the  magnet  attached  to 
the  structure  has  a  uniform  velocity  sensitivity.  This  may  not  be  true  for  our 
situation  but  is  probably  a  close  approximation.  Because  of  spatial  averaging, 
the  sensor  is  a  sort  of  low  pass  wavenumber  filter.  However,  at  wavenumbers  past 
the  cutoff  wavenumber  (defined  in  Equation  (2.12))  there  axe  regions  of  increased 
sensitivity  (i.e.,  sidelobes  in  the  filter  function). 

By  examining  the  wavenumber  response  of  a  single  sensor,  the  wavenumber 
filtering  action  is  obtained.  There  is  no  ability  to  select  (other  than  the  low 
pass  filtering  action)  which  particular  wavenumber  is  to  be  examined.  Using  an 
array  of  accelerometers  allows  for  a  certmn  degree  of  wavenumber  selectivity  to 
be  achieved. 

The  velocity  sampled  by  N  identical  transducers  is  just  a  variation  of 
Equation  (2.9)  with  the  addition  of  N  sensors.  A  schematic  for  this  situation  is 
shown  in  Figure  2.7.  Mathematically,  this  array  can  be  represented  as, 

y  N-l 

Um(a:)  =  [  /  v{a)h{x  -  a)da]  ^  gn^x  -  nAz) 

-io 

=  [t;(x)  *  /i(z)],4(z). 


(2.13) 
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Figure  2.7  Multiple  sensors  in  a  velocity  field. 


where 


N-i 

n=0 


(2.14) 


>l(x)  is  the  array  function  and  the  ^„’s  su-e  arbitrary  weights  applied  to  the 
outputs  of  the  sensors  which  axe  spaced  Ax  units  apart.  Taking  the  spatial 
Fourier  trainsform  of  Equation  (2.13)  and  mzdcing  use  of  the  convolution  theorem 
again  results  in 


=  [Vik,)H{k,)]*A{kt).  (2.15) 


For  an  array  of  ideal  sensors  the  velocity  spectrum  is 


Vr„ik,)  =  Vik,)*A{k,). 


(2.16) 


The  term  A{kx)  is  the  Fourier  transform  of  the  array  function  or 


N-i 

M^t)  =  53 

n=0 


(2.15) 
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The  term  A{kx)  represents  the  filtering  action  in  wavenumber  space  of  the  array 
of  sensors.  The  DFT  of  the  velocity, 

N-l 

\Uk,)  =  Ax  ^  ,  (2.18) 

n=0 

is  equivalent  to  Equation  (2.16).  The  DFT  is  the  most  efficient  way  to  caJctilate 
the  wavenumber  response.  Additionally,  in  almost  all  cases,  this  is  the  only  way  to 
calculate  the  wavenumber  response  for,  in  order  to  evEduate  Equation  (2.15),  the 
true  spectrum  V(kx)  must  be  known  a  priori.  However,  the  DFT  representation 
tells  nothing  about  the  filtering  action  of  the  array  as  do  Equations  (2.15)  and 
(2.16).  H{kx)  represents  the  wavenumber  sensitivity  of  a  single  accelerometer 
and  A{kx)  of  the  array.  H{kx)  is  a  multiplicative  term  that  rides  over  the  entire 
spectrum.  By  assuming  the  sensor  to  be  ideal,  the  wavelengths  must  be  larger 
than  the  sensor  (i.e.  the  criterion  of  Equation  (2.12)  must  be  met). 

For  real  sensors  the  wavenumber  sensitivity  cuts  off  after  a  certain  wavenum¬ 
ber,  because,  when  the  wavelengths  in  the  structure  are  small  relative  to  the  size 
of  the  sensor,  there  is  an  averaging  effect  which  tends  to  attenuate  the  voltage 
output  of  the  sensor.  Equation  (2.15)  represents  the  wavenumber  sensitivity  of 
an  array  of  real  sensors.  If  a  single  sensor  is  unable  to  measure  the  velocity  at  a 
pzu'ticulax  wavenumber,  then  the  array  of  sensors  is  likewise  unable  measure  the 
velocity  at  that  wavenumber.  Using  an  array  of  sensors,  a  wavenumber  band-pass 
filter  is  created  by  the  summation  of  the  outputs  of  the  sensors,  but  only  below 
the  cut-off  wavenumber  of  an  individual  sensor.  The  use  of  an  array  enables  a 
selectivity  within  the  range  of  the  measurable  wavenumbers  not  available  with 
just  one  sensor.  The  convolution  integral  of  Equation  (2.16)  is  an  integration  of 
a  shifted  version  of  the  array  function  times  .the  true  velocity  spectrum  which 
results  in  the  measured  spectrum  at  kg. 
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2.4.3.  Windowing  Considerations 

The  effect  of  windowing  the  output  of  the  sensors  measuring  the  velocity 
on  a  finite  structure  must  be  viewed  from  a  different  standpoint  them  for  an 
infinite  structure.  Using  window  functions  is  an  attempt  to  smooth  the  artificial 
discontinuity  caused  by  truncation  of  the  data  sequence.  For  the  finite  beam  the 
truncation  is  not  artificial,  the  velocity  distribution  on  the  beam  truly  exists  only 
over  a  finite  length,  outside  this  domain  the  velocity  is  zero. 

It  is  useful  to  define  the  rectangle  function  as 

r  0,  for  |(x  -  xo)/i>|  >  1/2; 

rect(-  )  =  <  1/2,  for  |(x  -  xo)/i>|  =  1/2  (2.19) 

i  1,  for  |(x  -  xo)/6|  <  1/2. 

The  amplitude  of  the  Fourier  transform  of  the  rectangle  is  the  sine  function 
defined  as, 

=  5!^.  (2.20) 

There  is  an  additional  multiplicative  phase  factor  if  the  rect  function  is  not 
centered  at  x  =  0.  If  a  pure  sinosouidal  wave  exists  only  over  the  length  of  the 
structure  it  can  be  represented  as 


u(x)  =  sin(A:iox)rect(^),  (2.21) 

where  the  Fourier  transform  of  u(x)  is 

Uikr)  =  ^{6{k,  -  k,o)  -  S{k,  +  fc,o)}  *  {sinc(A:,L/2)}.  (2.22) 

2j 

Let  U{kx)  be  the  true  velocity  distribution  on  some  finite  structure,  then  the 
measured  wavenumber  spectrum  can  be  obtained  from  Equation  (2.16).  Here 
it  is  postulated  that  for  a  finite  structure  with  am  array  of  closely  spaced  point 
sensors  the  proper  weighting  is  a  uniform  one  (i.e.  all  of  the  g„’s  equal  to  one). 
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Hence,  is  the  sine  function  and  Equation  (2.16)  becomes, 

Umikr)  =  Uik,)*smc{k,L/2), 

=  -^^iHkz  -  kzo)  —  6{ki  +  Itio)}  *  {sinc(fcrI,/2)}  ♦  {sinc(fciL/2)}. 

(2.23) 

The  convolution  of  two  sine  functions  of  the  same  length  results  in  one  sine 
function  or, 

siixc(krL/2)  *  smc(kxL/2)  =  smc{kxL/2).  (2.24) 

It  follows  that, 

U„,{k,)  =  U{k,).  (2.25) 

In  this  case,  the  true  spectrum  is  measured.  Of  course  in  reality  A{kx)  only 
approximates  a  sine  function.  But  in  the  case  where  the  sampling  is  very 
fine,  the  array  function  closely  approximates  the  ideal  case.  Thus,  it  is  argued 
(heuristically)  that  the  proper  weighting  is  uniform. 

Of  coiurse  Equation  (2.25)  is  an  expected  result.  However,  if  the  spatial 
window  function  (the  gn’s)  is  other  tham  a  uniform  window,  even  if  the  sampling 
interval  is  infinitesimally  smaJl,  Equation  (2.25)  would  not  hold.  If  the  velocity 
field  is  sampled  across  the  entire  length  of  the  beam,  then  the  only  errors  inherent 
in  the  DFT  are  due  to  aliasing  [14].  If  some  spatial  weighting  is  used,  obviously, 
the  estimate  will  vary  from  the  continuous  transform  due  to  both  adiasing  and 
the  window. 

In  the  measurement  of  the  Fourier  transform  of  the  velocity  on  a  finite 
structure  using  the  DFT,  aliasing  cannot  be  avoided  ais  the  spectrum  is  not 
baind  limited.  The  effects  of  aliasing  cam  be  reduced  by  choosing  a  spacing  that 
moves  the  aliasing  lobes  of  the  array  function  to  a  region  of  wavenximber  space 
where  the  spectrum  has  very  low  values.  In  Section  5.3,  the  effect  of  varying  the 
spacing  of  the  accelerometers  is  exiunined.  The  goal  in  the  array  design  is  to 
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use  as  fine  a  spacing  eis  possible  and  to  obtain  an  array  function  whose  Fourier 
transform  is  close  to  the  proper  sine  function. 

2.4.4.  Spacing  of  Accelerometers  for  the  Experiment 

The  dispersion  relation  for  the  beam,  which  is  derived  in  Chapter  3  and 
pictured  in  Figure  3.2,  shows  that  the  highest  freebending  w'avenumber  expected 
is  about  0.6  cm“^.  The  spectrum  for  the  finite  beam  falls  off  as  a  sine  fimction 
centered  at  the  freebending  wavenumber.  In  other  words  the  spectrum  should 
fall  off  around  30  dB  after  three  octaves  or  at  kz  =  2.4  cm“^ .  Hence,  a  spacing  of 
1  cm,  which  yields  a  Nyquist  rate  of  it  cm~^  and  would  reduce  aliaising,  is  used. 

The  other  item  with  which  to  contend  is  the  criterion  of  Equation  (2.12).  The 
accelerometer  used  has  a  radius  of  0.5  cm.  The  criterion  for  this  size  transducer 
is  that  the  minimiun  wavelength  should  be  greater  than  2.1  cm.  The  measured 
wavenumber  spectra  would  exhibit  a  slight  attenuation  compared  to  an  exact 
analytic  treatment,  except  of  course  for  the  adietsing  that  would  tend  to  accentuate 
the  higher  wavenumbers.  It  is  left  to  the  results  section  to  examine  which,  if 
either,  effect  evidences  itself. 

2.4.5.  Wavenumber  Content  of  an  Array  of  Point  Forces 

Each  of  the  shakers  in  the  array  is  considered  delta  functions  in  space.  Thus, 

the  array  of  M  shakers  where  qm  is  the  complex  force  aunplitude  and  im  is  the 

shaker  location  can  be  represented  as, 

M 

q{x)  =  ^  qmS{x  -  Xm), 

”*=*  (2.26) 
M 

m=2 

where  Xm  =  mAx.  The  force  is  zero  everywhere  except  at  a  shaker.  Thus,  the 
external  forcing  function  is  completely  described.  The  Fourier  transform  may  be 
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performed  exactly  and  is 

M 

Q{k,)  =  ^  (2.27). 

m=l 

This  is  the  wavenumber  spectrum  of  the  array  of  point  forces. 

2.4.6.  Steering  the  Wavenumber  Content  of  an  Array  of  Point  Forces 

In  order  to  obtain  better  measurements  of  the  velocity  at  a  particular 
wavenumber,  it  is  helpful  to  concentrate  the  energy  of  the  forcing  function  at 
that  wavenumber.  The  objective  is  to  steer  the  peak  of  the  forced  wavenumber 
spectrum  to  the  desired  or  drive  wavenumber,  fcj.  The  peak  of  the  forced 
wavenumber  spectrum  may  be  steered  to  by  amplitude  shading  of  the  forced 
amplitudes. 

As  mentioned  previously,  this  is  done  by  setting  the  q^s  to  cos(fcd( im  —L/2)), 
which  modulates  the  spatial  signal  so  that  there  are  peaks  in  the  wavenumber 
spectrum  at  kz  =  This  holds  if  the  Nyquist  criterium  is  met,  i.e.  kd  < 

tt/Ai.  There  must  be  at  least  two  shaikers  for  each  wavelength.  Additionally, 
there  are  peaks  in  the  spectrum  at  kg  =  ±(27r/Az  —  kd)-  The  location  of  these 
peaks  can  be  predicted  by  examining  the  convolution  of  the  Fourier  transform 
of  the  cosine  function  and  the  array  of  point  forces.  The  finite  spatial  aperture 
of  the  forcing  function  will  broaden  the  peaks  of  the  force  spectrum  and  cause 
the  appearance  of  sidelobes.  The  peaks  at  kg  =  ^kd  are  called  main-lobes  while 
the  secondary  peaks  are  called  grating  lobes;  collectively,  the  main-lobes  and 
the  grating  lobes  are  referred  to  as  major- lobes.  The  sidelobes  axe  smaller  in 
amplitude  than  the  major-lobes  and  fall  between  the  major-lobes.  Lastly,  the 
force  spectrum  of  an  array  of  point  forces  is  periodic,  with  period  equal  to  2ir/ Az. 
This  characteristic  is  important  as  the  pattern  seen  in  the  region  0  <  kg  <  27r/ Az 
is  repeated  periodically  in  the  positive  and  negative  kg  directions. 
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\Vhile  it  is  easy  to  theoretically  specify  the  force  amplitudes,  it  is  a  completely 
different  matter  to  experimentally  enforce  these  amplitudes.  The  voltages  applied 
to  the  shakers  tire  controllable,  so  that  the  drive  wavenumber  Ccin  be  steered  to 
the  desired  wavenumber  by  steering  the  drive  voltages.  This  assumes  that  the 
force  applied  to  the  beam  tracks  with  the  voltage  in  wavenumber  space  and  a 
complicated  feedback- control  system  between  the  voltage  applied  and  the  force 
output  can  be  avoided.  For  the  remainder  of  this  section,  the  qm's  (the  force 
output  at  the  shtdcer)  is  assumed  to  be  controllable. 

Some  results  of  the  wavenumber  spectra  of  the  forces  seen  at  the  shakers  (i.e., 
the  applied  voltages)  are  now  presented  for  the  spacing  amd  number  of  shakers 
used  in  the  experiments.  For  the  nine  shaiker  experiments  (Ai  =  4  cm),  the 
wavenumber  spectrum  for  the  unsteered  force  distribution  (Jt^  =  0  cm”')  and  a 
steered  force  distribution  with  kd  =  0.332  cm”'  (tv/o  wavelengths  over  the  length 
of  the  beam)  are  shown  in  Figure  2.8a  and  2.8b,  respectively.  The  peaks  at  ±kd 
and  at  fc,  =  27r/Ax  —  kd  can  be  seen.  It  is  importaint  to  note  that,  due  to  the 
periodic  nature  of  the  spectrum,  by  prescribing  a  drive  wavenumber  between  0 
amd  tt/ Ai  a  peak  force  wavenumber  response  may  be  steered  to  any  wavenumber. 
There  aie  always  an  infinite  number  of  peaiks.  This  fact  is  both  am  advantage 
amd  a  drawback.  Without  the  grating  lobes,  the  only  way  to  excite  the  very 
high  wavenumbers  is  to  decreaise  the  spacing,  which  requires  smaller  shatters. 
However,  the  grating  lobes  can  complicate  the  interpretation  of  the  spectrum,  by 
exciting  the  velocity  in  an  unwamted  region. 

Results  of  the  force  wavenumber  spectra  for  seventeen  shakers  au:e  shown  in 
Figure  2.9.  For  seventeen  shakers  Ai  is  equal  to  2  cm.  Note  that  the  grating 
lobes  do  not  appear  in  these  plots  because  the  grating  lobes  aire  outside  the  chosen 
wavenumber  range.  By  examining  these  plots  it  is  interesting  to  note  that  the 


Wavenumber  (1/cm) 


(b) 

Shifted  and  unshifted  force  spectra  for  seventeen 
sheikers,  (a)  kd  =  0.  cm“‘,  (b)  kd=.332  cm“^ 
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width  of  the  major-lobes  is  a  function  of  MAx,  or  the  length  of  the  spatial  signal 
and  any  spatial  windowing  which  may  be  added  to  the  shalcer  weighting  during 
steering  (this  is  described  next).  Hence,  if  the  spatial  windowing  is  the  same  for 
the  nine  shakers  as  for  seventeen,  the  major  lobes  have  the  same  width. 

While  the  steering  of  the  voltage  input  moves  the  peak  of  the  wavenumber 
content  to  a  desired  wavenumber,  the  high  sidelobes  caused  by  rectangular 
windowing  may  produce  undesirable  effects.  These  sidelobes  may  cause  higher 
velocity  response  than  the  main-lobes  if  the  sidelobes  occur  in  a  region  of  high 
admittance.  By  applying  a  windowing  function,  the  sidelobes  are  reduced  at  the 
cost  of  a  broader  main-lobe.  The  values  of  the  voltages  for  the  windowed  aind 
shifted  distribution  become. 


qm  =  W{xm)cosiki{x^  -  1/2))  ,  (2.28) 

where  W{x)  is  the  selected  window  function.  Different  windows  may  be  applied 
to  the  shaher  forces  in  order  to  alter  the  sidelobe  structure  which  is  done  at 
the  sacrifice  of  increasing  the  width  of  the  mtiin-lobe.  A  Kaiser-Bessel  window 
with  it’s  variable  parameter  alpha  set  to  1.5  is  used  in  the  results  presented  in 
Figure  2.10.  This  window  type  is  described  in  reference  [16].  By  windowing,  the 
energy  of  the  forcing  function  is  less  likely  to  leak  into  the  nezirby  wavenumbers. 
The  periodic  mzun-lobes  are  not  aiffected  except  for  broadening.  In  Figure  2.10, 
kd  =  0  cm~^  and  a  Kaiser-Bessel  window  is  used  with  nine  and  seventeen  shal:ers. 

The  FORTRAN  program  FORCE-FILTER,  is  used  to  obtain  the  wavenum¬ 
ber  spectrum  and  the  correct  amplitude  shading  for  all  of  these  cases.  This  code 
is  found  in  Appendix  C. 


dB  re  max.  value  dB  re  max.  value 
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(b) 

Figure  2.10  Effect  of  Kaiser- Bessel  window  on  unshifted  spectra, 
(a)  for  nine  shakers,  (b)  for  seventeen  shakers. 


Chapter  3 

BEAM  THEORY 
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3.1.  Introduction 

In  the  experiments,  a  beam  is  excited  by  multiple  shakers.  The  number  of 
shakers  is  varied  and  the  response  is  measured  at  accelerometer  locations  spaced 
along  the  beam.  The  accelerometer  data  Eire  used  to  obtain  the  wavenumber 
response  via  the  discrete  Fourier  transform.  In  order  to  validate  suid  understand 
the  experimental  results,  the  theoretical  response  of  the  bcEim  to  multiple  forces  is 
developed  in  both  the  spatial  and  wavenumber  domains  according  to  Timoshenko 
beam  theory. 

The  normal  mode  solution  for  the  velocity  of  a  Timoshenko  beam  is  briefly 
developed  in  this  chapter  and  the  results  given  by  Hutto  [9]  are  restated.  These 
results  allow  the  validation  of  the  metisured  velocity  at  a  point  due  to  an  arbitrary 
forcing  function;  this  is  called  the  forced  spatial-frequency  response  of  the 
beam.  The  theoreticail  wavenumber  response,  defined  eis  the  Fourier  transform 
of  the  velocity,  may  be  obtEiined  by  tsdcing  the  discrete  Fourier  transform  of 
the  predicted  results  or  by  performing  the  Fourier  trEinsform  EuiEdytically.  The 
wavenumber-frequency  response  of  an  infinite  Timoshenko  beam  is  derived  to 
obtain  the  freebending  wavenumbers  and  as  a  trEuisition  step  in  understEinding 
the  relationship  between  force  and  velocity  in  wavenumber  space.  Finally,  the 
theory  of  the  finite  Fourier  trEinsform  is  applied  to  obtain  predictions  of  the 
wavenumber-frequency  admittance  of  the  finite  free-free  Timoshenko  besun. 

3.2.  Spatial-Frequency  Response  of  a  Finite  Beam 

The  coupled  differentiad  equations  with  internal  damping  for  a  Timoshenko 
type  beam  with  the  totEil  transverse  deflection,  y,  and  the  slope  due  to  bending. 


u\  axe  [9,17] 


pA&^yix.t)  , 

- ^ - ) - -  =q{x,t) 


dx'^ 


dx 


dt^ 


E.j3h^  +  »G-.4(^  -  =  0 
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(3.1) 


(3.2) 


dx^  •  '  dx  -V-,-./ 

where  q  is  the  forcing  function,  the  complex  elastic  modulus  is  E*  =  E(1  +  jq), 
the  complex  shear  modulus  is  G*  =  G(1  +  jt;),  /  is  the  second  moment  of  inertia 
of  the  cross  sectional  area,  ,  p  is  the  density,  ,4  is  the  cross  sectional  area,  and 
K  the  numerical  shape  factor  [18].  The  loss  factor,  q,  is  taken  to  be  equal  for 
rotation  and  translation.  Only  the  time  harmonic  case  is  considered,  with  time 
dependence  u  being  radian  frequency  and  j  =  v^— 1.  Iviciking  the  following 
substitutions  in  Equations  (3.1)  and  (3.2) 


y(x,f)  =  r(x)e>‘"‘  , 

(3.3) 

t/^(x,<)  =  ^'(x)e>‘^‘  , 

(3.4) 

9(x,t)  =  Q(x)e^‘^‘  , 

(3.6) 

and  omitting  the  time  dependence  yields. 


KG*Ar"(x)  -  «G*A^'(i)  +  pAJ^Y{x)  =  Q{x)  ,  (3.6) 

E*I^"{x)  +  kG*A{Y'{x)  -  ^'(i))  +  plu^^ix)  =  0.  (3.7) 

Where  the  primes  indicate  differentiation  with  respect  to  x  the  spatial  variable. 
These  are  the  dimensional  time  hzurmonic  differential  equations  which  are  used 
later  in  this  chapter.  In  the  development  of  the  finite  beam  response,  it  is  useful 
to  non-dimensionalize  the  above  equations.  The  non-dimensional  variable  ^  is 
used  where. 


C  =  x/I. 


(3.8) 
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Making  the  substitution  into  Equations  (3.6)  and  (3.7),  the  non-dimensional. time 
harmonic  differential  equations  are, 

+  pAu^Yd)  =  Qd)  ,  (3.8) 

^«"(0  +  KG-.4(ir'(«  -  S'CO)  +  pIwHd)  =  0.  (3.9) 

The  primes  here  indicate  differentiation  with  respect  to  the  To  obteun  the 
normal  mode  equations,  first  the  homogeneous  problem  is  solved.  In  references 
[9,17],  the  homogeneous  forms  of  Equations  (3.8)  tind  (3.9)  are  decoupled  in  Y 
and  resulting  in. 


Y""  b\r^  +  s‘^)Y'’  -  b\l  -  b^r^s‘^)Y  =  0,  (3.10) 

_  62(1  _  62^2^2)^  ^  0^  (3 


where 


62  ^ 

El 


r2  = 


AL^ 


El 

kAGL^ 


The  general  solutions  for  Y  amd  ^  aure. 


(3.12) 

(3.13) 

(3.14) 


F(^)  =  Cl  cosh(6a^)  -t-  C2  sinh(6a^)  -I-  C3  cos{bl3^)  C4  sin(6/3^)  ,  (3.15) 

^(0  =  sinh(6a^)  -h  C2  cosh(6aO  +  C'^  sin(6/?0  +  cos{b3^)  ,  (3.16) 

where 

^  =  ;^{T(r^  +  5^)  +  [(r^ +  •  (3-17) 

The  general  relationship  between  Ci  and  C'  is  given  in  [9].  Only  the  relationship 
between  Ci  and  C[  is  needed  to  describe  the  bending  slope;  that  relation  is  given 
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later.  The  boundary  conditions  for  the  free-free  besun  applied  at  the  ends  of  the 
beeim  ^  =  0  and  ^  =  1  are  that  the  bending  moment  is  zero  or, 


(3.18) 


and  the  total  shear  force  is  zero  or. 


-Y'-'H  =0 

1j 


(3.19) 


Applying  the  boundary  conditions  to  Equations  (3.15)  and  (3.16)  results  in  a  four 
by  four  matrix  to  solve  for  three  of  the  unknown  C's.  By  setting  the  determinant 
to  zero  the  chairacteristic  equation  is  found;  it  is 


2  -  2cosh(6a)cos(fc^) 


[6^r^(r^  —  -I-  (3r^  —  s^)]  sinh(6Q;)  sin(6;5)  =  0  .  (3.20) 


This  transcendent  ad  equation  is  solved  for  the  naturad  frequencies  of  the  beaun. 
In  the  computer  code  that  solves  for  the  natural  frequencies  it  is  important 
to  apply  a  high  argument  approximation  to  the  hyberbolic  functions  (that  is 
sinh(i)  =  cosh(x)  as  x  becomes  lau’ge).  Otherwise,  errors  occur. 

The  mode  shapes  are  now  known  to  within  a  multiplicative  constant.  Let 
the  subscript  i  denote  a  mode  number  corresponding  to  a  natural  frequency 


Define  the  following  constants. 


a?  -H  q|  + r2  /??  _ 


a?  +  s2  —  r2  /??  _ 

A  -  ^ 

^~0i  ’ 


aj  -h 


(3.21) 


(3.22) 


^  _  cosh(6iai)  —  cos(b,j9,) 

‘  Ajsinh(6iai)  —  Ci  sin(6i/3,) 

Using  the  notation  introduced  above,  the  mode  shapes  are 


(3.23) 


YiiO  =  Ci[cosh{biai^)+  —  cos{bi0,^)-Xi8ismh{b,ai^)-S,sm{bi0iO]  i  (3.24) 

(>• 
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'I'i('C)  =  Cj[sinh(6,a,0  -  A,  sin(6,|S,^)  -  A.(5i  cosh(6iQ,0  -  6iA,C,  cos(6,y5,^)]  . 

(3.25) 

where  the  approximations  to  the  hyperbolic  functions  of  large  tirgument  axe  used 
for  the  mode  shapes  also.  When  i,Q,  »  1,  <5,  — ►  ^  and  when  6,ai^  1 

cosh(feiQ,^)  —  sinh(6,a,^)  — +  0,  which  results  in  the  reduction  of  Equations  (3.24) 
and  (3.25)  to 

i'iCO  =  C'i[^cos(6i/3iO  -  -J-sin(i)i/9iO]  >  (3.26) 

Si  ^i 

-  2  2 

'J>,(^)  =  -Ciy — — — [A,  sin(5,/9,0  + '5.A,Ct  cos(6,/3jO]  (3.27) 

L  a 

Where  C[  has  been  substituted  for  in  terms  of  C\.  The  resontint  wavenumbers 
of  the  Timoshenko  beam  axe  bi^i.  Later,  it  will  be  seen  that  these  are  the 
free-bending  wavenumbers  of  the  solution  to  the  infinite  beam  that  also  fit  the 
free-free  boundary  conditions  to  produce  resonance  in  the  beam. 

Since  the  boundary  conditions  for  the  free-free  beam  allow  for  rigid  boc 
motion,  there  exists  non-trivial  solutions  for  the  zero  natural  frequency.  The  two 
rigid  body  modes  are 

Yoi  =  Co,  ,  (3.28) 

Yo2  =  Co2(e  -  1/2)  .  (3.29) 

The  orthogonality  conditions  for  the  free-free  Timoshenko  beam  axe  given  in 
[19].  They  axe 

=  pAMi  i  =  j 


where  M,  is  the  normalization  constant  for  the  ith  mode. 
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3.3.  Forced  Spatial-Frequency  Response  of  a  Finite  Beam 

In  the  normal  mode  solution  to  the  forced  or  non- homogeneous  Timoshenko 
beam  equations,  Y  and  are  expressed  in  terms  of  the  eigenmodes  of  the 
homogeneous  solution, 

OO 

r(0  =  x^a.v;(0  ,  (3.31) 

1=0 


1=0 


The  constants  a,  are  the  modsil  expansion  coefficients  to  be  determined.  The 
constant  Ci  is  absorbed  into  the  Oj’s.  By  substituting  Equations  (3.31)  and  (3.32) 
into  the  non- homogeneous  differential  Equations  (3.6)  and  (3.7),  eind  using  the 
orthogonality -integral  expressions  the  unknown  Oj’s  are  found.  The  details  are 
found  in  reference  [9],  The  general  form  of  the  solution  for  Oj  due  to  an  arbitrary 
forcing  function  is. 


pAMi[u>fil  +  jrj)  -  u;2] 


(3.33) 


The  integration  represents  a  meEisure  of  the  orthogonzdity  between  the  forcing 
function  and  the  I'th  normal  mode.  If  the  force  distribution  is  the  same  as  a 
mode  shape,  then  that  mode  may  have  a  very  high  modal  constant.  Of  course, 
the  steady  state  frequency  u  plays  an  important  role;  mode  shapes  whose  natural 
frequencies  are  far  removed  from  the  applied  frequency  have  little  effect  on  the 
beam  velocity. 

The  experimental  measurements  are  taken  using  shakers  that  closely  approx¬ 
imate  point  sources.  In  the  i  direction  the  extent  of  the  shakers  is  smedl,  and 
the  width  of  the  shaker  nearly  spans  the  width  of  the  beam.  So  the  shakers  are 
line  sources  in  three  dimensions,  but,  since  the  beam  is  approximated  as  a  one 
dimensional  continuum,  a  forcing  function  consisting  of  am  array  of  such  shakers 
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can  be  modeled  as  an  array  of  point  sources,  so  that  Q(0  is 


A/ 


Q(0  =  E  -  «"•)  ■ 


(3.34) 


m=l 


where  qm  is  the  complex  amplitude  of  the  force  located  a  ^  Using  the 

above  in  Equation  (3.33)  results  in, 


A/ 


*  *  IL 

at  =  Di  ^  ~[cosh(biai^tn)  +  —  cos(fc,/?,^,n)  -  A,(5,- smh(!>,o,-^m) 


mssl 


or, 


where. 


Af 


m=l 


-i  =  Di  Y, 


-«,sin(i.iftU))  ,  (3.35) 


(3.36) 


I>.= 


1 


(3.37) 


pAMi{u}j ( 1  +  ;  t; )  -  u;2 ) 
and  yi(4m)  is  value  of  the  mode  shape  evaluated  at  ^  It  should  be  noted 

that  experimentally,  the  values  q,n  are  functions  of  drive  frequency.  The  modal 
constants  for  the  translational  and  rotational  rigid  body  modes  are  respectively, 


M 


Ooi 


m=l 


(3.38) 


M 


002  =  Do2  -  -5)  .  (3.39) 

m=l 

This  is  the  analytic  model  developed  to  validate  experiments  for  the  velocity  at 
a  point  due  to  an  array  of  point  sources  exciting  a  beam  at  a  given  frequency. 

3.4.  Wavenumber  Response  of  a  Forced  Finite  Beam 


In  this  section  the  spatial  Fourier  transform  of  the  velocity  response  of  a 
Timoshenko  beam  excited  by  multiple  point  drives  is  derived.  The  spatial  Fourier 
transform  is. 

OO 

F{k,)=  j  /(x)e--'**V- 

— OO 


(3.40) 
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and  the  inverse  transform  is  . 

OO 

/(x)  =  ^  /  .  (3.41) 

—  OO 

The  wavenumber  is  kx  and  x  is  the  spatial  variable.  In  order  to  obt2iin 
the  wavenumber-frequency  response  of  the  beam,  Equation  (3.40)  is  applied 
to  Equation  (3.31).  The  velocity  distribution  of  Equation  (3.31)  exists  from 
0  >  ^  >  1  and  is  taken  to  be  zero  elsewhere.  Thus,  the  limits  of  integration 
in  Equation  (3.32)  tire  now  zero  and  one.  Applying  the  Fourier  transform  to 
Equation  (3.31)  yields 


Yik^)  =  J 

0 


=  (3.42) 

1=0 

F,(Jt^)  is  the  spatial  Fourier  trzmsform  of  each  mode  shape,  obtained  by  applying 
the  Fourier  transform  to  Equation  (3.24).  For  the  non-rigid  body  modes  ^iik^) 
is 


=  .2  r~/\.  ~^2  (c"''**  bh  cosh{biai)  -f-  6iQisinh(6iai) 

“i  )  V 

-  Xi6i{jk^  sinh(6ia<)  -f-  6,aj  cosh(6iaj))]  -  {jk^  -  (5,A,6jQ,)| 


+  t2 — ■'*Ml/Ci(;^'{Cos(6i/3i)  -  6,/?, sin(6i/?,)) 

-  Si{jk^sin{bi0i)  -H  6,/3,  cos(6i/3i))]  -f  (jfc{/Ci  “  .  (3.43) 


or. 


y;(fc?)  =  ra(A:e)  +  r.2(fc€)  ,  (3.44) 
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where  represents  the  first  six  terms  of  Equation  (3.43)  and  r.2(h) 

represents  the  second  six  terms.  There  are  two  special  cases  to  consider.  One, 
when  the  factor  r,2(t^)  becomes, 


r.2(A:^) 


.  r  1  1  e  1 

2j 


|-(  — 


e-2jfce 


) 


■  1 

Ak^Ci 


(3.45) 


Two,  when  biOi  becomes  large  with  respect  to  the  numerical  precision  of  the 
computer,  r,i(fc^)  reduces  to 


r,i{k() 


-j 

k^+j{biai) 


(3.46) 


For  the  rigid  body  modes,  the  transforms  of  the  translational  and  rotational 
modes  axe 

=  f  -  1)  .  (3.47) 

k( 

yo2(t()  =  -  1)  +  +  1)  ■  (3.48) 

The  transform  of  each  mode  shape  consists  of  one  component  due  to  the 
trigonometric  functions  zmd  another  due  to  the  hyberbolic  functions.  The 
contribution  due  to  the  hyperbolic  functions  is  basictdly  a  linear  decay  in 

• 

wavenumber.  The  sine  and  cosine  parts  produce  a  sin(x)/x  function  which  has 
a  peaJc  at  k^  =  bifii.  The  6,/3i  terms  tire  ctdled  the  resoneint  wavenumbers  of 
the  beam.  Equation  (3.43)  is  the  transform  of  a  single  mode  and  is  a  relatively 
complicated  form;  complicated  in  that  even  if  a  single  mode  was  excited  on 
the  structure  the  Fourier  spectrum  would  not  consist  of  delta  functions  at  the 
resonant  wavenumber.  The  dependence  of  the  velocity  on  the  forcing  function  is 
held  in  the  Oj's 

The  treinsform  is  in  terms  of  the  non-dimensional  wavenumber  k^.  It  is 
desirable  to  express  the  transform  in  terms  of  the  dimensionail  wavenumber, 
kx.  The  problem  of  finding  the  mapping  relation  between  the  transform  of  the 
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displacement  in  terms  of  the  dimensional  and  non-dimensional  wavenumber  is 
caste  in  terms  of  the  functions  /  and  /.  The  proof  presented  here  is  general 
for  functions  where  the  relation  between  /  and  /  is  simply,  fiO  =  /(-^O  and  ^ 
is  defined  to  be  equal  to  x/L.  These  are  the  conditions  that  relate  the  spatial 
dependence  of  the  dimensional  and  non-dimensioned  forms  of  the  displacement. 
The  relationship  between  the  Fourier  transforms  of  /  and  f  is  desired.  The 
following  definitions  are  made 

OO 

nk()=  I  . 

—  OO 

OO 

F(Jtr)=  J  /(x)e--'*=**dx  . 

—  OO 

In  order  to  find  how  the  vzdues  of  F  may  be  obtained  from  F,  first  the  change 
of  variables  x  =  is  made  in  the  second  definition  resulting  in 

OO 

F(fcx)  =  I  J  . 

—  OO 

Using  the  relationship  between  /  and  /  to  obtain, 

OO 

nk.)  =  L  j  /(e)e-'‘-“<ie  , 

—  OO 

and  finally  using  the  first  definition  yields  the  desired  result 

F(fcx)  =  IF(fcxI)  .  (3-49) 

Similarly  one  may  obtain, 

FCfcj)  =  jnk(/L)  .  (3.50) 

This  is  the  relation  by  which  the  values  of  F  may  be  obtained  from  F.  For 
the  quotient  of  two  functions  of  the  transformed  variable  the  multiplicative  L 
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term  would  cancel.  The  method  of  proof  is  suggested  by  the  deriration  of  the 
similarity  theorem  in  [15]  and  also  by  [20]. 

The  exact  Fourier  transform  of  the  spatial  response  of  the  forced  Timoshenko 
beam  serves  as  a  basis  for  comparison  for  experimental  and  theoreticed  estimates 
of  the  Fourier  spectrum. 

3.5.  Wavenumber  Response  of  an  Infinite  Beam 

In  this  section,  predictions  of  the  wavenumber  response  of  an  infinite  beam 
are  derived  using  Timoshenko  beam  theory.  Using  this  result,  the  dispersion 
relation  for  an  infinite  beam  is  obtained.  The  derivation  basically  follows  the 
development  of  Magrab  [21]. 

Equations  (3.6)  and  (3.7),  the  dimensional  differential  equations,  re-written 
here, 

kG^AY'\x)  -  kG*,4'^'(x)  +  pAu^Y{x)  =  Q{x)  , 

E'I<il”{x)  -b  kG^A{Y\x)  -  ^(i»  -b  pJa;2^(x)  =  0  , 

axe  the  starting  point  of  the  analysis.  The  Fourier  transform  as  defined  in 
Equation  (3.40)  is  used  to  transform  the  above  into  wavenumber  domain.  Let 
represent  the  Fourier  transform  operator,  the  following  identity  from  Bracewell 
[15],  J^{f'{x)}  =  jkt^{f{x)}  is  used  in  obtaining, 

Y{k^)[klKG\A  -b  pAu;^]  +  ^(it,)[-j/kr«GM]  =  Qik^)  (3.51) 

Y{k^)l-jk^KG^A]  -b  ^{kr)[-k^E*I  +  plu^  -  kG’A]  =  0  ,  (3.52) 

after  applying  the  Fourier  transform  to  Equations  (3.7)  emd  (3.8).  Define  the 
following  terms; 

fl  =  pAu)^  , 

0  =  pluj'^  ,  (3.53) 
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1/  =  kG“  A 

Also,  let  the  freebending  wavenumber  for  the  Bernoulli  Euler  beam  as  defined  in 
Junger  and  Feit  [22]  be  =  Q/E*I.  Use  the  expressions  in  Equation  (3.53) 
to  simplify  Equations  (3.51)  and  (3.52)  to 


+  n]  +  ^ik,)[-jk,i.]  =  Q(k,)  (3.54) 


y{f^r)[jk^u]  +  ^{k^)[-krE*I  +  0  -  ly]  =  0  .  (3.55) 

Equations  (3.54)  and  (3.55)  represent  two  coupled  algebraic  equations  in  ^'(ki) 
and  ’^{kj.].  By  using  Cramer’s  rule  both  unknowns  may  be  found.  Applying 
Cramer’s  rule  to  solve  for  U(^'i)  results  in 

'  +  «£•/)  + file -H  ' 


Dividing  top  and  bottom  by  E*Ii/  results  in, 

y,,  .  ^  Q(k.){-^{Q/u-l)-kl/u] 

^  kj  -  ki[e/E*i  +  n/u)  +  k%^[e/u  - 1] 


(3.56) 


The  wavenumber-frequency  admittance  is  defined  as  the  quotient  of  the  Fourier 
transforms  of  the  velocity  and  the  force,  for  the  infinite  beam 


Hooikz)  =  jw 


-^{efu-D-kllu 


kt  -  ki[e/E-i  +  Q/u]  +  k%E[Q/^  - 1] 


(3.57) 


where  the  multiplicative  ju)  comes  from  the  conversion  of  displacement  to  velocity. 
Thus,  the  wavenumber-frequency  response  of  an  infinite  beam  to  an  arbitrary 
forcing  function  has  been  derived  according  to  the  assumptions  of  Timoshenko 
beam  theory.  For  a  Bernoulli-Euler  type  beam  0  — ►  0  and  i/  —*■  oo  {u  as  defined 
in  (3.53)  not  Poisson’s  ratio)  [22],  and  Equation  (3.57)  becomes 


y{kz) 


-Q{kz)/E*I 

~  ^BE 


(3.5S) 
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From  this  it  can  be  seen  that  ksB  is  the  wavenumber  at  which  the  admittcince 
for  the  Bernoulli-Euler  beam  has  its  peadc. 

The  zeros  of  the  denominator  of  Equation  (3.58)  are 

7..2.3.4  +  +  ~  ‘ 


noting  kgE^/^  =  E*Iv.  The  7’s  axe  the  poles  of  the  solution. 

In  order  to  interpret  the  meaning  of  the  poles  of  the  admittance,  the  inverse 


Fourier  transform  of  the  solution  must  be  examined.  Subject  to. 

0  ^  17®  TTT 


(3.60) 


and  the  damping  taking  the  form  E*  =  E{1+jt]),  the  poles  of  the  Equation  (3.59) 
are  located  approximately  as  shown  in  Figure  3.1.  If  the  damping  is  zero,  then 
the  roots  of  the  denominator  are  either  purely  real  or  imaginary.  The  inverse 
Fourier  transform  solution  for  the  displacement  is, 


00 

nx)  = 

J.TZJU)  J 


Since  the  forcing  function  under  consideration  is  a  series  of  delta  functions,  it 
suffices  to  examine  the  solution  for  a  single  point  force  located  at  i  =  xo,  in 
which  case  the  preceding  equation  simplifies  to, 


00 

Y(x)  =  5^  / 


(3.61) 


This  integr£d  may  be  evaluated  using  the  theory  of  residues.  To  use  residue 
theory,  the  integramd  must  decay  at  infinity  and  consist  of  only  outgoing  waves. 
Hence  there  axe  two  caises.  One,  for  z  —  Xo  <0,  where  the  integrand  is  evaluated 
using  the  two  residues  in  the  upper  half  plane  (UHP  in  Figure  3.1).  Two,  for 
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Figure  3.1 


Location  of  the  poles  and  contours  for  the  residue 
evaluation. 


X  —  xo  >  0  the  integrand  is  evaluated  using  the  two  residues  in  the  bottom  half 
plane  (BHP  in  Figure  3.1). 

It  can  be  seen  in  Equation  (3.61)  that  the  poles  near  the  real  axis  represent 
the  propagating  part  of  the  velocity  field  and  the  poles  near  the  imaginary  axis 
the  evanescent  field.  The  freebending  wavenumber  at  a  peurticuleu'  frequency  is 
defined  as 

kf  =  7j  (3.62) 

It  can  be  easily  shown  that  the  freebending  wavenumber  should  be  the  same 
as  the  resonant  wavenumbers  at  the  resonant  frequencies.  Upon  examining  the 
definitions  of  Equations  (3.12)-(3.14)  and  Equation  (3.53),  the  following  relations 


axe  found 


6^x2  =L^Q/E'L 
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52  —L*kg£. 

Using  th  :se  in  Equation  (3.59)  the  free-bending  wavenumbers  results  in, 
71,2,3,4  =  ±  ±  \/(r2  -  s2)2  +4/62j 

-i:  b0/ L,ijba/ L 


(3.63) 


The  resonant  frequencies  are  just  those  frequencies  for  which  the  freebending 
wavenumber  satisfies  the  boundary  conditions.  It  is  very  useful  to  exzunine  the 
overall  dispersion  relation  in  order  to  better  tell  what  wavenumbers  to  expect  at 
the  frequencies  of  interest.  In  Figure  3.2,  the  resonant  wavenumbers  are  plotted 
using  solid  dots  and  the  freebending  wavenumber  for  an  infinite  Timoshenko 
beam  using  a  solid  line.  The  freebending  wavenumber  for  a  Bernoulli- Euler  beam 
is  also  plotted  in  the  same  figure.  In  Table  3.1,  the  resonant  frequencies  and  the 
resonant  wavenumbers  for  the  beam  are  listed.  The  dispersion  diagram  is  helpful 
in  two  ways.  First,  it  provides  a  guide  in  spacing  the  accelerometers  to  measure 
the  velocity  field.  Also,  this  plot  may  be  used  to  see  what  wavenumbers  might 
be  excited  in  a  beam.  In  Section  2.4.4,  the  highest  freebending  wavenumber  that 
is  expected  for  the  velocity  is  read  from  this  figure.  The  ctdculations  were  made 
using  the  parameters  for  the  beam  used  in  the  experiments  (see  Chapter  2  for 
material  constants). 


3.6.  Relation  between  Force  and  Velocity  in  Wavenumber  Space  for  a 
Finite  Beam 


.4.S  seen  in  Equation  (3.57),  the  admittance  of  an  infinite  beam  is  a 
that  function  of  only  the  physical  and  geometric  properties  of  the  beam  and 


Frequency  (kHz) 

20.0  30.0  40.0  50.0  60.0  70.0 
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*  *Resonant  wavenumbers 
~  Infinite  Timoshenko  beam 
^Infinite  Bernoull-Euier  beam 


The  dispersion  relation  for  the  beam  used  in  the 
experiment.  The  freebending  wavenumbers  for  the 
infinite  beam  are  shown  as  well  as  the  resonant 
wavenumbers  (and  frequencies)  for  the  finite  beam. 


Figure  3.2 
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independent  of  the  forcing  function.  However,  Equations  (3.33)  and  (3.42)  imply 
that  the  admittance  function  for  a  finite  beam,  formed  by  dividing  the  velocity 
wavenumber  spectrum  by  the  force  wavenumber  spectrum,  is  a  function  that 
depends  on  the  form  of  the  forcing  function  as  well  cis  the  properties  of  the  beam. 
Because  we  are  forced  to  deal  with  finite  beams,  quantification  of  the  dependence 
of  the  admittance,  computed  as  the  ratio  of  the  velocity  and  force  wavenumber 
spectra,  on  the  physical  properties  of  the  beam  as  well  as  the  properties  of  the 
forcing  function  is  required. 

For  linear,  shift  invariant  systems,  the  impulse  response  function  completely 
determines  the  system.  Both  the  finite  and  the  infinite  beam  are  linear  systems, 
but  only  the  infinite  beam  is  a  shift  invariant  system.  Physicedly,  shift  invariance 
means  that  only  the  distance  between  the  shalcer  and  the  accelerometer  2ire 
needed  to  define  the  response,  not  their  absolute  locations.  For  the  finite  beam 
the  location  of  force  and  the  sensor  are  needed. 

In  an  attempt  to  quantify  the  meaning  of  the  quantity  Y{kx)/Q{kx)  for  a 
finite  beam,  the  finite  Fourier  transform  is  applied  to  the  governing  differential 
equations  of  the  Timoshenko  type  beam.  As  defined  in  Bracewell,  the  finite 
Fourier  transform  is 


F[k.,a.b)=  f  . 

J  o 

The  inverse  finite  Fourier  transform  takes  the  form  of  a  series, 

/(i)  =  5^  E  n2nx/L.a,by^'"'‘/‘-  . 


(3.64) 


(3.65) 


The  finite  Fourier  trzmsform  of  f'{x)  the  spatial  derivative  of  /(i)  is 

j' f{x)e-^^^^dx  =  jkxF(kx,a,b)  +  -  /(a)e-^'^'‘'  .  (3.66) 

Higher  order  derivatives  may  also  be  obtained  from  the  above  formula.  This 
relation  is  different  than  relations  for  the  infinite  transform.  The  interval 
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under  consideration  is  the  non-dimensionaJ  length  of  the  beam  0  <  ^  <  1,  for 
convenience  and  brevity  the  following  notation  is  introduced. 


F/(^'rr=F(;-„0.1)  . 


(3.67) 


The  transform  variable  for  the  non-dimensiontd  equations  is  and  the  trains- 
forms  of  the  displacement,  bending  slope  and  force  are  defined  in  Equation  (3.67) 
or, 


Yfih)  = 

Jo 

(3.68) 

^fih)  = 

Jo 

(3.69) 

Qf{k^)  = 

f' . 

Jo 

(3.70) 

Applying  the  definition  of  the  finite  Fourier  transform  to  the  non-dimensional 
equations  of  motion  for  the  Timoshenko  beam  (Equations  (3.8)  and  (3.9))  yields 

-  + nY,{ko  =  Q,(k()  . 


(3-71) 


By  using  the  first  boundary  condition,  given  Equation  (3.10) 


yr'(0  - ’J'(0  =  0  for  e  =  0  and  ^  =  1 

1j 

Equation  (3.71)  simplifies  to 

=  Qf(h)  -  ■  (3.72) 

Taking  the  finite  Fourier  transform  of  the  second  differential  and  using  the  second 
boundary  condition  (^''(0)  =  ^'^(l)  =  0)  results  in, 

+ 0  -  H  =  [^no  -  . 


(3.73) 
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The  terms  that  contain  the  values  of  the  displacement  and  bending  slope  at  the 
endpoint  are  contained  in  the  following  terms. 


Siih)  =  (3.74) 

(3.75) 

Using  the  definitions  in  Equations  (3.74)  and  (3.75)  to  rewrite  the  transformed 
differential  equations  in  a  form  ready  for  solution  of  the  unknowiis  yields 


+  =  ,  (3.76) 


17(»-<)l-Y^l  +  ^'/(lf)l^^'(+e->'l=Bj(ts)  •  (3.77) 

These  equations  are  solved  using  Cramer’s  rule  for  Yf{k^),  which  results  in  , 


(ib^)  =  [0/(^-g)  +  gi(fcg)][-f/fc^  +  0  -  t/]  +  B2{hhk^^ 

-  ‘'I 


(3.78) 


Using  the  fact  that  Hoo(kz)  =  Hoo{k^/L)  the  above  equation  can  be  written  as 

:)  •  (3.79) 

The  admittance  for  the  finite  beam  is 


Y,(k,)  =  ^H^(k,IL)\Q,[k,)  +  +  E-;  .  >• 

+  7^(0  -  i') 


Hf{k^)  = 


juYf(k^) 

Qfih) 


The  wavenumber  response  of  the  finite  beam  consists  of  a  part  that  is  the  same  as 
for  an  infinite  beam  and  another  peurt  that  depends  on  the  boundary  \’alues  of  the 
displacement  and  bending  slope.  Since  the  constants  Bi{k^)  and  B2{k^)  depend 
on  the  values  of  the  spatial  variables,  to  solve  for  the  wavenumber  response  the 
spatial  response  must  be  obtained  as  well.  Although  it  is  not  done  here,  the 
constants  would  be  evaluated  by  inverting  the  Equation  (3.7S)  and  considering 
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F(0),  1^(1),  'P(O)  and  ^'(1)  as  unknown  constants.  Upon  inversion  these  constants 
could  be  evaluated  and  then  substituted  into  the  wavenumber  response. 

There  are  several  implications  from  Equation  (3. SO).  First,  the  wavenumber 
admittance  of  the  free-free  Timoshenko  beam  is  dependent  on  the  location  £ind 
amplitude  of  the  forces  applied  to  the  beam.  This  means  that  a  wavenumber 
admittance  that  is  a  function  of  only  the  physical  and  geometric  properties  of 
the  beam  is  not  obtained  simply  by  forming  the  quotient  of  the  wavenumber 
velocity  spectrum  and  force  spectrum.  The  infinite  beam  part  of  the  admittance 
is  im'ariant  with  respect  to  chcinge  in  the  forcing  function.  At  wavenumbers 
where  the  amplitude  of  the  forcing  function  is  high,  the  admittance  measured 
with  the  finite  beam  approximates  the  infinite  beam  wavenumber  admittemce. 
The  velocity  response  of  the  beam  at  a  particular  wavenumber  is  a  function  of 
how  near  this  wavenumber  is  to  the  freebending  wavenumber  and  how  well  this 
wavenumber  matches  the  boundary  conditions.  When  both  of  these  conditions 
are  met  resonance  occurs.  Since  the  impedance  of  the  beam  is  low  at  resonance, 
it  is  expected  that  the  amplitude  of  the  force  wavenumber  spectra  would  be  low  if 
the  forcing  function  matches  the  structural  mode  shape.  With  a  matched  forcing 
function,  the  finite  beam  admittance  at  resonance  is  a  poor  approximation  to 
the  infinite  beam  admittance.  It  is  important  to  note  that  even  at  a  resonant 
frequency  the  resulting  velocity  contribution  from  the  resonant  mode  shape  would 
be  small  if  the  vzilue  from  the  spatial  integration  that  occurs  in  Equation  (3.33) 
is  small,  i.e.,  if  the  force  distribution  does  not  match  the  mode  shape. 


63 

Chapter  4 

EQUIVALENT  CIRCUIT  MODELING  OF 
THE  SHAKER  FORCE/GAGE  UNITS 

4.1.  Introduction 

In  this  chapter,  the  response  of  the  combination  shaker /force  gage  units 
is  modeled  using  equivalent  circuits.  The  basic  Meison  equivalent  circuits  for 
piezoelectric  and  non-piezoelectric  materials  are  shown.  From  these  circuits  the 
low  frequency  approximations  are  derived.  First,  the  equivalent  circuit  for  the 
overall  force  output  of  the  transducer  is  developed.  From  this  model,  the  levels 
and  resonance  frequencies  of  the  force  output  due  to  an  arbitrary  termination 
impedance  are  found.  Also,  the  relationship  between  the  voltage  output  from 
the  force  gage  section  and  the  force  applied  to  the  structure  is  obtained  through 
the  use  of  equivalent  circuits.  The  effect  on  the  force  gage  output  due  to  added 
cable  capacitance,  mass  loading  and  contact  stiffness  is  obtmned.  The  method 
of  calibration  and  its  consequences  are  also  be  examined. 

4.2.  Equivalent  Circuit  Analysis  for  Transducers 

The  dynamic  response  of  shzikers  can  be  lepresented  by  using  Mason 
equivzilent  circuits.  These  circuits,  in  their  exact  form,  represent  a  solution 
to  the  electriczd  zind  elastic  differential  equations  that  govern  the  transducer 
system.  For  piezo-active  elements,  a  catalogue  of  lumped  impedance,  equivalent 
circuits  for  different  boundary  conditions  and  piezoelectric  couplings  have  been 
developed  [23].  The  analysis  of  electro-mechanical  circuits  in  general  is  found  in 
Beranek[24].  The  Mason  equivalent  circuits  of  the  piezoelectric  elements  and  the 
non-piezoactive  elements  may  be  cascaded  together  to  form  the  model  for  the 
entire  transducer  system.  In  this  study,  the  low  frequency  approximation  to  the 
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exact  solution  is  used. 

By  constructing  the  complete  circuit  representation  of  the  trzinsducer,  the 
resonance  frequency,  output  force  aind  force  sensitivity  of  the  system  cem  be 
approximated  cis  well  as  the  effect  of  changing  any  of  the  transducer  peirameters. 
The  equivalent  circuit  analysis  yields  solutions  for  the  desired  quantities,  the 
overall  force  output  and  the  force  gage  sensitivity,  in  terms  of  the  material 
properties  of  the  shaker,  excitation  frequency  and  load  impedance. 

The  equivalent  circuits  are  a  representation  of  the  solution  of  the  differential 
equations  governing  the  elements  of  the  transducer.  Certeun  boundary  conditions 
have  been  applied  to  the  non-active  areas  (in  this  case  the  radial  faces).  Basically, 
the  boundary  conditions  are  that  the  radial  faces  axe  stress  free  and  the  electric 
field  lines  do  not  fringe.  The  wave  motion  through  the  transducer  is  assumed  to 
be  a  plane  wave  traveling  in  the  zudal  direction.  Because  the  length  to  diameter 
ratio  is  large,  the  response  of  the  shaker  is  almost  entirely  in  the  axiad  direction. 
Berlincomt  et  al.  [23]  describe  in  detail  the  boundary  conditions  assumed  for 
the  equivalent  circuits  mentioned  above.  The  non-piezoactive  elements  of  the 
transducer  axe  modeled  as  longitudinally  vibrating  rods  or  bars. 

The  equivalent  circuits  utilized  axe  of  the  mechanical  impedance  type.  In 
the  mechanical  impedance  representation,  the  flow  quantity  is  velocity  and  the 
potential  quantity  is  force.  The  boundary  conditions  for  the  force  and  velocity 
transmitting  surfaces  are  left  as  axbitreury  impedances.  Each  section  of  the 
transducer  that  has  different  material  properties  will  be  modeled  using  a  similar 
equivalent  circuit,  but  with  different  material  properties.  The  rest  of  the  model 
can  be  constructed  by  cascading  the  equivalent  circuits  one  after  the  other. 

The  equivalent  circuit  for  the  33  coupled  piezoceramic  is  shown  in  Figure  4.1a 
The  equivalent  circuit  for  the  non-piezoactive  elements  is  shown  in  Figure  4.1b 


Figure  4.1 


Zo  =  pAvh, 


Equivalent  circuit  for  (a)  piezoelectric  33  coupled 
sections,  (b)  Non-piezoactive  sections  [23]. 
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[23].  In  these  figures,  U  is  the  velocity  and  F  is  the  force.  These  are  the 
solutions  to  the  differential  equations  for  the  so  called  ‘length  expander’  case. 
The  values  for  the  impedances  in  the  equivalent  circuits  are  calculated  from  the 
formulas  given  below  the  circuit  diagrams  in  Figure  4.1.  The  nomenclature  used 
to  describe  the  constants  defining  the  impedances  and  the  impedances  themselves 
is  a  combination  of  that  used  in  [23,25].  Since  piezoceramics  are  anisotropic 
materials  the  notation  used  in  the  definition  of  the  constants  use  indicia!  notation. 
Additionally,  since  there  are  both  electrical  eind  elastic  boundary  conditions  to 
-.e  denoted,  there  must  exist  a  flexibility  in  the  nomenclature  to  describe  these 
conditions. 

The  numerical  subscripts  (for  example  ^33)  have  the  following  meaning.  The 
first  index  indicates  the  axis  to  which  the  electroded  surface  is  perpendicular, 
or  the  direction  of  the  poling  of  the  crystal.  The  second  index  indicates  the 
direction  the  forces  will  be  applied  (either  piezoelectrically  induced  or  externally 
applied).  The  coordinate  system  for  a  crystal  is  shown  in  Figure  4.2  and  is  a 
local  coordinate  system  for  each  element  of  the  trainsducer.  There  are  several 
different  edphabetic  subscripts  juid  superscripts.  The  superscripts  D  and  E 
are  used  in  reference  to  elastic  constants,  E  means  that  the  constants  were 
measured  with  a  constant  electric  field,  and  D  that  the  charge  density  was 
constant.  The  superscripts  T  and  S  me  used  in  reference  to  the  electric  or 
piezoelectric  constants,  T  means  that  the  constants  were  measured  with  the 
free  elastic  boundary  conditions,  and  5  with  the  boundary  conditions  clamped. 
These  boundary  conditions  are  implicit  in  the  solution  of  the  governing  equations 
and  assumed  not  to  change.  For  exeimple,  the  radial  surface  of  the  33  coupled 
transducer  is  assumed  to  be  stress  free,  hence  the  electric  properties  are  indicated 
with  the  T  superscript.  The  subscript  h  used  in  conjunction  with  the  wavenumber 
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and  the  wave  velocity  vf  indicates  that  the  length  expander  or  bax  mode  is 
the  prominent  mode.  With  the  subscript  and  superscript  notation  described, 
the  material  constants  are  defined.  The  constants  ^33  and  dss  are  piezoelectric 
constants.  The  load  bearing  area  and  the  electroded  area  (which  axe  the  same 
for  the  piezoelectric  materials  with  the  33  coupling)  is  denoted  cls  A  and  the  axial 
length  as  /.  The  material  density  is  p.  Co  is  the  electrical  capacitcince  and  0  or 
Ip  will  be  used  in  this  chapter  to  denote  the  transformation  factor  which  converts 
from  the  electrical  domain  to  the  mechamical.  The  relative  dielectric  constant  of 
the  piezoceramic  is  C33  and  the  free  space  constant  is  eo-  The  elastic  compliemce 
is  S33  which  is  the  inverse  of  the  elastic  modulus. 


Figure  4.2  Local  coordinate  system  for  the  ceramic. 

The  representations  in  Figures  4.1a  and  4.1b  are  the  solutions  to  the 
governing  differentitd  equations  and  have  no  frequency  limitation  except  that 


68 


the  boundary  conditions  will  breaJc  down  at  high  frequency.  These  represen¬ 
tations  become  cumbersome  to  program  and  it  is  efficacious  to  use  the  smedl 
argument  approximation  for  the  sine  and  cosine  functions.  The  v’alidity  of  these 
assumptions  breaks  down  when  the  error  in  the  small  curgument  approximation  in 
any  one  of  the  elements  becomes  appreciable.  For  example,  one  would  expect  the 
small  argument  approximation  to  become  invalid  in  the  decoupler  section  first. 
The  small  argument  model  (which  amounts  to  a  low  frequency  approximation)  is 
shown  in  Figure  4.3a  for  the  piezoactive  element  and  Figure  4.3b  for  the  elastic 
element.  This  model  reduces  to  lumping  the  static  mass  and  stiffness  of  the 
elements  of  the  transducer  into  the  branches  of  the  circuit.  This  approach  allows 
some  engineering  approximations  to  be  made.  When  the  frequency  is  low,  the 
stiffness  and  mass  of  some  adjacent  elements  may  be  combined  to  simplify  the 
resulting  equivalent  circuit  without  a  significant  loss  of  accuracy.  The  overall 
transducer  model  is  obtained  by  cascading  the  equi\’alent  circuits  for  each  section 
one  after  the  other.  To  be  precise,  every  different  material  of  the  transducer 
should  be  modeled  separately.  This  could  result  in  a  very  leirge  equivalent  circuit 
with  many  branches.  Thus,  the  initial  iterations  of  the  model  are  made  by 
making  simplifying  2issumptions  about  what  are  the  important  elements  of  the 
treuisducer.  For  example,  with  respect  to  the  overall  force  output  of  the  shaker, 
the  force  gage  plays  a  small  part.  So  in  order  to  remove  a  loop  of  the  total 
equivalent  circuit,  the  mziss  and  compliance  of  the  force  gage  cerzimic  zure  lumped 
with  the  headmass  and  the  decoupler  sections. 

4.3.  Equivalent  Circuit  Model  for  the  Shaker  Force  Output 

The  transducer  described  in  Section  2.3  is  modeled  in  this  section.  The 
modeling  process  is  iterative.  Complexity  is  added  to  the  model  until  predicted 
results  match  measurements. 
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Figure  4.4  shows  the  equivalent  circuit  for  the  entire  shaker,  force  gage 
transducer  system  using  the  low  frequency  approximations  of  Figures  4.3a  and 
4.3b.  Since  the  branches  of  each  section  of  the  circuit  form  a  T,  these  circuits 
axe  called  T  circuits.  The  first  T  circuit  represents  the  drive  section  and  the 
tailmass  section.  The  compliances  of  the  two  PZT  disks  and  the  breiss  tailmeiss  eire 
included  in  Ci .  At  low  frequencies,  these  elements  act  as  a  combined  stiffness  [12]. 
The  mass  of  the  tailmass  and  one  half  the  mass  of  a  PZT  disk  zire  lumped  into  Mi . 
The  mass  of  one  and  one  half  PZT  disks  Eire  included  in  the  calculation  of  Mj. 
The  second  T  circuit  includes  the  effect  of  the  decoupler  alone.  The  compliance 
C2  is  the  compliance  of  the  decoupler.  Half  of  the  decoupler  mass  is  lumped  into 
each  of  M2  Emd  M3 .  The  third  T  circuit  includes  the  effect  of  the  headmass  and 
the  PZT  force  cell  in  a  manner  similar  to  the  first  T  circuit.  The  effect  of  the 
added  mass  due  to  the  glue  used  to  stress  relieve  the  lead  wire  attachments  is 
added  to  Mi ,  M2  and  M3  in  proportion  to  their  estimated  amounts.  The  values 
of  the  elEistic  and  piezoelectric  constants  needed  for  cEilculation  of  the  impedances 
in  the  equivalent  circuit  axe  found  in  Table  4.1  and  4.2. 

Explicitly,  the  lvalues  of  the  impedances  in  Figure  4.4  are  calculated  from; 
Ml  =  1/2  mass  of  PZT  driver  +  meiss  of  brass  teiilmass, 

M2  =  3/2  mass  of  PZT  driver  +  1/2  mass  of  decoupler, 

M3  =  l/2(  meiss  of  decoupler  +  mass  force  gage  +  mass  of  headmass), 

M\  =  l/2(  mass  force  gage  +  mass  of  headmass), 

Co  =  electrical  capacitance  of  PZT  driver. 

Cl  =  mechanical  compliance  of  drive  section, 

C2  =  compliance  of  decoupler, 

C3  =  compliance  of  headmeiss  +  force  gage, 

Zt  =  termination  impedance, 
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Figure  4.4  Equivalent  circuit  model  for  the  force  output  of  the 
transducer. 


4>  =  transformation  factor. 

These  values  are  calculated  using  the  formulas  given  in  Figure  4.2.  After 
combining  similar  impedances  and  reflecting  electrical  elements  through  the 
transformers  to  convert  them  into  the  mechtinical  domain,  the  circuit  takes  the 
form  of  Figure  4.5.  The  compliance  of  this  branch  is  C*  which  now  takes  into 
account  the  effect  of  the  electrictii  impedance. 

The  loop  equations  for  this  circuit,  using  the  currents  shown  in  Figtire  4.5 
(noting  that  current  is  velocity  tmd  voltage  is  force)  are; 


(4.1c) 
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Table  4.1  Nominal  material  properties. 


Material 

Elastic  compliance  (m*/N) 

Density  (kg/m^) 

Aluminium 

1.41x10““ 

2700 

Brass 

9.62x10“^^ 

8500 

PZT-4 

not  applicable 

7750 

Table  4.2  Matericd  properties  for  PZT-4. 


Material  Property 

Value  (units) 

g33 

26.1x10“^  (Vm/N) 

d33 

289x10““  (m/N) 

S33 

7.9  X  10“^'^  (m^N) 

S33 

15.5  xlO"^'^  (m^/N) 

^33/  ^0 

1300 

0  —  /aC  H - —  +  Zf).  (4. Iff) 

These  equations  can  be  caste  in  matrix  form,  and  solved  using  Cramer's  rule 
for  any  of  the  unknowns.  The  potential  drop  across  the  termination  impedance 
represents  the  force  applied  to  the  structure.  In  the  next  section  the  relation  of 
the  force  gage  voltage  output  to  the  force  is  examined.  In  the  current  model  the 
force  applied  to  the  structure  is, 

•Foul  =  (4.2) 


The  results  of  applying  the  equivalent  circuit  model  to  the  shaker  are 
compaured  with  theory  in  Figure  4.6.  The  results  show  good  agreement  especially 
in  the  low  frequency  range.  The  theoretical  resonance  frequency  is  29.8  kHz. 
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Figure  4.5  Equivalent  circuit  model  for  the  forced  output  of  the 
transducer  with  electrical  elements  trtinsformed  into 
the  mechauiical  domain. 

while  the  experimental  value  is  29.6  kHz.  The  fact  that  the  analytic  shaker 
model  does  not  include  any  damping  or  loss  mechtinism  accounts  for  some  of  the 
discrepancy  of  the  amplitude  between  experiment  eind  theory.  Additionally,  the 
error  associated  with  the  low  frequency  approximation  is  becoming  appreciable 
at  30  kHz.  At  this  frequency,  the  error  eissociated  with  the  low  argument 
approximation  is  8  percent. 

4.4.  Equivalent  Circuit  Model  for  the  Force  Gage  Section 

In  Section  4.3,  the  force  gage  section  was  modeled  as  having  only  high 
frequency  resonances,  hence  the  details  of  its  voltage  response  were  ignored. 
In  this  section,  the  equivalent  circuit  for  the  force  gage  section  driven  by  an 
arbitrary  force  into  an  tirbitrary  impedance  is  examined.  The  purpose  of  this 
section  is  to  determine  the  response  of  the  force  gage,  and  its  relationship  to  the 


0.0  2.S  5.0  7.S  10.0  12.S  15.0  17.S  20.0  22.5  25.0  27.5  30.0  32.5  35.0 

Frequency  (kHz) 
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^measured  force  output 


Figure  4.6  Transmitting  voltage  response  for  unloaded  shaker, 
theory  versus  experiment. 
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force  applied  to  the  structure.  Also,  calibration  procedures  are  examined  Eind  the 
normal  assumptions  regarding  the  sensitivity  of  force  gages  (that  the  sensitivity 
is  constant  with  respect  to  frequency  and  load  impedance)  axe  analyzed. 

The  model  for  the  force  gage  is  shown  in  Figure  4.7.  It  has  been  constructed 
using  the  low  frequency  approximations  in  Figure  4.3a  aind  4.3b.  The  force  source 
represents  the  rest  of  the  transducer,  which  is  modeled  as  some  arbitrary  input 
force.  In  the  calculation  of  the  sensitivity  of  the  force  gage,  if  the  input  force 
cancels  out,  then  the  sensitivity  is  independent  of  the  drive  section  and  decoupler 
section  of  the  shaker.  Note  in  Figure  4.7,  that  since  the  electrical  terminals  eire 
open,  the  positive  and  negative  capacitainces  cancel  if  the  cable  capacitance, 
C cable,  is  zero  or  small  .  Even  though  the  capacitances  are  in  series,  their  net 
effect  is  zero  because  the  positive  and  negative  impedances  of  the  two  capacitors 
cancel  when  added.  First  the  effect  of  the  cable  capacitance  is  ignored,  then 
a  correction  factor  is  listed  to  take  into  account  the  effect  of  the  finite  cable 
capacitance.  In  Figure  4.7,  Mj,  Mj,  Cj,  C2,  Co* 

Mj  =  1/2  (mass  of  PZT  force  gage  +  aluminum  headmeiss), 

M2  =  1/2  mass  of  aluminum  headmass. 

Cl  =  compliance  of  the  force  gage, 

C2  =  compliance  of  the  aluminum  headmass, 

Cq  =  electrical  capacitance  of  the  force  gage, 

Zt  =  termination  impedance, 
xj)  =  transformation  factor. 

The  sampled  voltage  is  in  the  electric  domain  and  we  must  take  into  account 
the  transformation  factor  in  calculating  the  output  voltage.  The  current  that 
flows  through  the  capacitor  across  which  the  output  voltage  is  measured  is  rp 
times  the  velocity  of  that  branch  of  the  circuit.  The  output  voltage  of  the  force 
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Figtire  4.7  Eqmvalent  circmt  for  force  gage  sensitivity. 


gage  is  then 


The  output  force  is 


Eout  —  {h  —  h)"  • 

jujCq 


Eout  —  liZt. 


(4.3) 


(4.4) 


This  is  the  force  imparted  to  the  structure.  The  first  item  of  interest  is  the  ratio 
of  the  output  voltage  of  the  force  gage  to  the  force  applied  to  the  structure.  In 
order  to  obtain  this,  the  loop  equations  must  again  be  solved.  The  equations  for 
the  loop  currents  (velocities)  are 


0  =  f^2( — 7^)  +  h{ju}M2  +  — ^  +  Z|).  (4.5c) 

JUi^C/2 


Upon  solving  the  equations  using  Cramer's  rule  for  the  velocities,  the  ratio  of 
the  output  voltage  to  the  output  force  is 


Font 


which  is  just  Equation  (4.3)  divided  by  Equation  (4.4).  A  constant  proportional¬ 
ity  between  the  two  is  desired.  The  solution  of  the  above  results  in  the  following: 


— t/)Cj  a/. 

R\  =  - (1  4-  JLJ  — 


C’ 


Zt 


K'l-y  2  \  t'  /-i'  3  A/i  A/oCn  . 


(4.6) 


The  common  factor  simplifies  to  g^z^/A  which  is  the  low  frequencv  sensitivitv 
of  the  certimic.  This  is  the  nominal  sensitivity  that  we  try  to  measure.  However, 
away  from  u  =  0  the  response  of  the  output  voltage  to  the  actual  force  appears 
to  become  more  complicated.  Actually,  most  of  the  other  terms  are  negligible 
in  the  frequency  range  exaanined  in  this  study.  If  the  force  gage  would  be  ideal, 
in  other  words  massless  and  infinitely  stiff,  C2,  Af(,and  would  be  zero  and 
the  sensitivity  would  be  independent  of  Zt  and  frequency.  The  non-zero  values 
complicate  matters.  The  terra  multiplying  ju  is  due  to  the  mass  loading  by  the 
mass  between  the  force  gage  and  structure.  The  actuail  force  that  reaches  the 
structure  is  less  by  the  velocity  times  ju;(Afi  -t-Mj).  This  factor  may  be  subtrr  'ted 
out  electronically  during  the  meeisurement  or  via  computer  software  afterward. 
At  the  lower  frequencies,  this  term  is  the  most  significant  source  of  errors  in  the 
measurement  of  the  force.  This  error  is  usually  neglected.  It  is  a  major  source 
of  error  only  when  the  impedance  of  the  structure  is  low;  i.e.,  at  a  resonance 
of  the  structure.  The  other  terms  in  higher  powers  of  u  have  8in  effect  only  at 
higher  frequencies.  If  these  effects  need  be  examined,  the  equivalent  circuit  of 
Figure  4.1a  and  4.1b  must  be  used.  When  the  admittance  of  the  structure  drops 
to  a  level  comp^able  to  the  mass  below  the  force  gage,  then  the  mass  subtraction 
should  be  used  if  precise  impedance  levels  are  required  (see  Appendix  B). 
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The  object  of  a  calibration  of  the  force  gage  is  to  obtain  the  nomintil 
sensitivity  of  the  force  gage  section.  A  standard  method  is  to  mass  load  the 
shalcer  and  mount  an  accelerometer  on  the  termination  mass.  Knowing  the 
mass,  its  acceleration  and  the  output  voltage  from  the  sensor  the  sensitivity 
can  be  inferred.  The  mass  that  loads  the  force  gage  is  M^  +  Mj  +  M2.  The 
sensitivity,  S/  is  approximately  the  voltage  outpu.  of  the  force  gage  divided  by 
the  acceleration  times  the  mass  loading.  Solving  for  this  we  obtain. 


Mt  +  Mo 


M[  +  M2  +  Mt 

this  is  the  sensitivity  measured  using  Equation  2.1.  If, 


(4.7) 


A/j  <<  M2  +  Mt', 


then  Equation  (4.7)  becomes 

=  .  (4.8) 

The  calibration  described  here  eind  in  Section  2.3.3  provides  a  good  estimate 
of  the  low  frequency  force  sensitivity.  If  the  compliance  C2  includes  only  the 
compliance  of  the  aluminum  headmass,  the  term  u!^M^C2  is  much  smaller 
than  one  even  at  100  kilohertz.  However,  the  glue  layer  attaching  the  shaker 
headmass  to  the  structure  to  be  tested  increaises  this  compliance.  The  sensitivity 
chzu’acteristically  rolls  off  at  higher  frequencies  due  to  this  contact  point  or  glue 
layer  compliance. 

Now  the  effect  of  cable  capacitance  for  the  33  coupled  force  gage  is  examined. 
The  cable  capacitance  cam  almost  double  the  mechanical  compliance  of  the  force 
gage  section  in  the  shaker  used  in  this  study.  This  is  not  a  problem  in  predicting 
the  entire  shaker’s  dynamic  response  because  the  compliance  of  the  force  gage 
will  still  be  much  less  than  the  compliance  of  the  decoupler  and  the  drive  sections. 
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Hence,  the  cable  capacitance  will  not  effect  the  dynzimics  of  the  system.  However, 
the  cable  capacitance  may  all  but  cripple  the  sensitivity  of  the  force  gage  if  the 
cable  capacitance  is  too  high.  The  cable  length  is  the  length  of  the  cable  to  the 
first  amplifier  or  buffer.  In  cases  where  the  cable  length  must  be  very  long,  a 
charge  amplifier  must  be  used.  To  account  for  the  effect  of  the  cable  capacitance 
the  following  changes  to  the  above  equations  are  made.  In  the  circuit  diagram 
in  Figure  4.7,  Cq  is  replaced  by 

Oq  =  Cq  +  Ccable  ) 

tind  mechanical  compliance  C[  is  everwhere  replaced  by  C*  which  is  found  from 

i/c*  =  i/c;-t^Vc'o'  +  i/’VCo. 


After  these  two  substitutions  are  made  in  the  circuit,  the  loop  Equations  (4.5) 
are  re-solved  eind  the  resulting  low  frequency  sensitivity  is 


5,'  = 


l+CcaWeC^-V’^IV) 


(4.9) 


Thus  if  Ccabu  is  zero  then  the  sensitivity  is  gzzt/A.  If  the  cable  capacitance  is 
very  high,  the  sensitivity  is  lowered. 

The  average  measured  sensitivity  for  the  force  gage  is  0.95  V/N  as  tabulated 
in  Table  2.3.  The  predicted  value  of  g^^tfA  is  1.02  V/N  which  matches  the 
measured  values  closely.  Average  cable  capacitzince  for  the  type  and  length  of 
cable  used  in  the  construction  of  the  force  gage  is  45  picofarads.  Plugging  this 
value  into  Equation  (4.9)  results  in  a  sensitivity  of  0.89  V/N,  again  close  to  the 
measured  value.  This  close  agreement  indicates  that  there  are  no  significant 
errors  in  the  construction  or  the  calibratioi:  of  the  transducer. 


Chapter  5 

EXPERIMENTAL  AND  THEORETICAL  RESULTS 
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5.1.  Introduction 

The  experimental  results  of  the  implementation  of  the  wavenumber  drive 
system  with  nine  and  seventeen  shakers  are  presented  in  this  chapter.  First,  the 
experimental  transfer  admittance  at  various  spaticd  locations  along  the  beeun  is 
compared  with  the  theoretically  predicted  values.  This  compcirison  allows  for  a 
point  by  point  validation  of  the  mezisurements.  Next,  the  effect  of  varying  the 
number  and  spacing  of  the  sensors  used  to  measure  the  wavenumber  spectrum 
is  analyzed  using  predicted  results.  Subsequently,  the  wavenumber-frequency 
dependence  of  the  mezisured  forcing  function  is  presented.  The  measured  and 
predicted  velocity  responses  to  this  forcing  function  are  examined.  Finally,  the 
relationship  between  the  wavenumber  content  of  the  forcing  function  and  the 
velocity  is  explored. 

5.2.  Measured  and  Predicted  Point  Admittances 

The  measured  admittance  at  various  locations  on  the  beam  is  compared  to  the 
predicted  admittance  as  a  function  of  frequency.  As  mentioned  in  Section  2.4.4, 
the  velocity  is  measured  at  37  locations  on  the  beam,  with  the  actual  locations, 
in  centimeters,  given  by 

Xi  =  ((i  -  1) -h  0.95)  , 
and,  in  non-dimensional  coordinates,  by 

=  ({i  ~l)  +  0.95)/37.9  . 

The  admittances  are  all  measured  with  respect  to  the  center  shaker,  so  all  but 
the  measurement  of  the  velocity  at  the  center  location  are  considered  transfer 
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admittances  defined  here  as 


Vi{u})  =  Vi{u>)/Fciuj)  ,  (5.1) 

where  u,  is  the  velocity  measured  at  location  i  and  Fc  is  the  force  measured  at 
the  center  shaher.  Since  the  beeim  is  being  driven  by  nine  or  seventeen  shalcers, 
the  admittance  takes  on  a  much  different  appearance  then  it  would  for  a  single 
drive  point. 

The  FORTRAN  program  VALTIM  (see  Appendix  C)  employs  the  bezun 
theory  in  Chapter  3  to  predict  the  admittance  of  the  beam  at  the  different  spatial 
locations  using  the  measured  forces  read  from  a  data  file.  VALTIM  also  reads 
the  measured  values  of  the  admittance  and  then  plots  both  the  theoretical  eind 
experimental  data  at  an  observation  point  as  a  fimction  of  frequency. 

The  forcing  condition  is  defined  by  the  experiment  number  given  in  Table  2.2. 
As  the  agreement  between  experiment  and  theory  is  generally  very  good, 
Experiment  7  has  been  chosen  as  a  representative  measurement.  Several  locations 
for  the  comparison  of  the  predicted  amd  measured  admittances  for  Experiment  7 
2u:e  presented  while  only  a  single  plot  is  shown  for  the  other  experiments.  The 
seventeen  shedcer  experiments  are  more  complicated  (i.e.,  more  sources  for  error), 
so  the  representative  data  set  was  chosen  from  the  last  five  experiments.  Also, 
many  of  the  features  for  the  nine  shalcer  experiments  are  the  same  as  for  the 
seventeen  shaher  experiments,  so  it  is  not  necessary  to  specifically  examine  in 
detail  all  of  the  data  sets  from  the  nine  and  seventeen  shaker  experiments. 

The  admittance  for  the  center  location  (^  =  0.5)  is  shown  in  Figure  5.1  from 
Experiment  7.  The  good  agreement  between  theory  and  experiment  is  seen.  The 
first  peak  around  3.1  kHz  is  due  to  the  system  resonance  caused  by  the  interaction 
between  the  beam  and  shakers.  At  this  frequency  the  force  and  velocity  at  nearly 
every  location  on  the  beam  hais  a  peak  in  amplitude  with  respect  to  the  force 
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Figure  5.1.  Comparison  of  the  theoretical  and  experimental  driv¬ 
ing  point  admittance,  magnitude  (dB  re  1  cm/dyne-s) 
and  phase  for  Experiment  7  at  the  location  ^  =  .5. 
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seen  at  the  center  shaJcer.  In  other  words,  the  driving-point  impedance  at  the 
center  of  the  beam  is  nearly  zero,  consequently  the  driving  force  is  nearly  zero. 
This  resonance  peaJc  occurs  around  3.1  kHz  for  all  of  the  experiments,  however 
the  location  of  this  peak  is  a  function  of  the  number  of  shakers  and  the  weighting 
pattern  of  the  drive  voltages.  The  shifting  of  this  peak  from  experiment  to 
experiment  may  be  due  to  mass  loading  of  the  force  gage  (for  a  brief  discussion 
of  mass  loading  see  Appendix  B),  however,  the  reasons  behind  the  shifting  of  this 
peak  are  not  investigated  in  detail. 

The  predicted  admittance  shows  a  peaJ:  at  the  resonance  frequencies  of 
3.6  kHz  and  8.4  kHz,  while  these  peaJcs  are  not  seen  in  the  experimental 
data.  These  are  even  or  symmetric  modes  (with  respect  to  the  center  of  the 
beam)  and  therefore  are  excited  more  by  the  symmetric  forcing  function  than 
the  antisymmetric  modes.  The  absence  of  these  predicted  resonance  effects  in 
the  experimental  data  could  be  due  to  the  use  of  a  damping  coefficient  in  the 
predictions  (0.05  is  used)  that  is  lower  them  the  actual  damping  coefficient.  This 
damping  is  due  to  having  shakers  glued  to  the  beam. 

Figures  5.2  and  5.3  show  the  trjinsfer  admittance  for  Experiment  7  at  a 
location  5  cm  to  the  left  and  to  the  right  of  the  center  of  the  beam  (^  =  0.368  and 
0.632).  If  the  forcing  function  were  perfectly  symmetric  these  admittances  would 
be  exactly  the  same.  While  they  are  quite  simileir  there  are  differences,  indicating 
that  the  forcing  function  was  not  perfectly  symmetric.  In  Figure  5.4,  the  results 
for  position  number  3  at  ^  =  0.078  are  presented.  While  the  magnitude  and  phase 
of  the  predicted  and  meeisured  drive  point  admittzinces  match  at  frequencies 
below  6  kHz  except  at  the  resonance  frequencies,  in  Figures  5. 2-5. 4  there  are 
deviations  between  the  predicted  and  measured  results  above  6  kHz. 
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Fieure  5.4.  Comparison  of  the  theoretical  and  experimental  trans¬ 
fer  admittance,  magnitude  (dB  re  1  cm/dyne-s)  and 
phase  for  Experiment  7  at  the  location  ^  =  .078. 
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The  transfer  admittances  at  location  14  =  0.368)  for  Experiments  1-6  and 

8-9  are  shown  in  Figures  5.5-5.12.  From  these  plots,  the  excellent  agreement 
between  theory  and  experiment  can  be  seen.  Again,  the  resonance  peaks  appeair 
in  the  predicted  results  that  are  not  found  in  the  experimental  results.  These 
peaks  are  found  around  3.6  kHz  (e.g.  Figure  5.7)  or  8.4  kHz  (e.g.  Fig\ire  5.10). 

The  suppression  of  the  resonance  peaks  in  the  meaisurements  is  notable. 
The  use  of  multiple  shakers  changes  the  usual  interpretations  of  the  admittance 
especially  when  the  forcing  function  is  weighted  to  a  particular  wavenumber. 
With  only  one  or  two  shakers  on  the  bezim,  the  mode  of  the  nearest  resonance 
frequency  ailmost  always  contributes  the  most  to  any  measured  response.  The 
exception  is  when  the  measurement  or  excitation  point  is  on  a  node  of  a  mode. 
With  multiple  shakers,  the  orthogonality  of  the  forcing  function  to  a  particular 
mode  may  result  in  that  mode  having  at  least  as  much  influence  on  the  velocity 
response  as  the  mode  whose  resonance  frequency  is  actually  closer  to  the  drive 
frequency.  Also,  the  gluing  of  nine  or  seventeen  shakers  to  a  beam  increases  the 
damping  of  the  beam  due  to  the  presence  of  the  glue,  which  also  suppresses  the 
measured  resonance  responses.  Thus,  w'ith  the  array  of  shakers  the  resonance 
may  be  suppressed  by  a  mismatch  between  the  spatial  response  of  the  resonance 
mode  and  the  spatial  distribution  of  the  applied  forces. 

Finally,  the  drive  wavenumber  for  Experiments  4,  8  and  9  corresponds  to 
a  freebending  wavenumber  whose  frequency  falls  far  outside  the  measurement 
frequency  remge  (compare  Tables  2.1  2md  3.2).  The  overall  admittance  for  these 
experiments  is  much  lower  than  for  the  other  experiments. 
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Figure  5.9  Compaxison  of  the  theoretical  and  experimental  trans¬ 
fer  admittance,  magnitude  (dB  re  1  cm/dyne-s)  and 
phase  for  Experiment  5  at  the  location  ^  =  0.368. 
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Figure  5.10  Comparison  of  the  theoretical  and  experimental  trzins- 
fer  admittance,  magnitude  (dB  re  1  cm/dyne-s)  and 
phzise  for  Experiment  6  at  the  location  ^  =  0.368. 
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Figure  5.12 


Comparison  of  the  theoretical  and  experimental  trans¬ 
fer  admittance,  magmtude  (dB  re  1  cm/dyne-s)  and 
phase  for  Experiment  9  at  the  location  ^  =  0.368. 
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5.3.  Eflfects  of  the  Spatial  Measurement  Array 

In  this  section,  the  effect  of  varying  the  accelerometer  spacing  on  the 
estimation  of  the  normalized  wavenumber  velocity  response  is  examined.  The 
normalized  wavenumber-frequency  velocity  spectrum,  is 

,  (5.2) 

where  V{kx,uj)  is  the  Fourier  trcinsform  of  the  velocity.  To  assess  the  effect  of 
spatial  sampling,  the  DFT  is  performed  on  the  theoretically  predicted  spatial 
velocity  response  and  compared  to  the  an£ilytic  or  continuous  transform  as  given 
in  Section  3.4.  The  DFT  is  calculated  using  a  128  point  FFT;  the  number  of  the 
actual  data  points  varies  but  the  data  length  is  always  zero-padded  to  a  length  of 
128.  Only  the  64  bins  corresponding  to  the  positive  spectrum  axe  plotted.  The 
FORTRAN  program  ARRAY -EFFECT  implements  the  calculation  of  both  the 
analytic  and  the  discrete  transforms  then  plots  the  results. 

Before  examining  the  results  of  these  calculations,  the  Fourier  treinsforms  of 
the  array  pattern  or  A{kx),  as  defined  in  Chapter  2,  is  plotted  for  spacings  of 
1  cm,  2  cm  and  4  cm  in  Figures  5.13-5.15.  A{kx)  also  is  called  the  array  pattern. 
Only  the  positive  part  of  array  fimction  spectrum  is  plotted,  the  negative  part  of 
the  spectrum  is  symmetric.  Re-writing  Equation  (2.12)  in  terms  of  the  current 
variables, 

V^{kx)  =  V{kx)*A{kx)  .  (5.3) 

This  is  the  measured  wavenumber  velocity  spectrum  (after  dropping  the  fre¬ 
quency  dependence).  In  this  equation,  the  convolution  integral  for  the  measured 
response  at  a  wavenumber  kx  can  be  interpreted  as  the  integration  of  a  shifted 
version  of  A{kx)  times  V{kz)  over  all  wavenumbers.  In  Figures  5.13-5.15,  the 
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Figure  5.15  The  Fourier  transform  of  the  array  function  for  nine 
accelerometers  with  a  spacing  of  4  cm. 
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mziin-lobe  of  A(kx)  is  located  at  kx  =  0cm“^,  the  lobes  appearing  at 

k'x  =  ±27rl/Ax  for  /  =  1, 2, ... 

axe  called  aliasing  lobes.  As  the  measurement  wavenumber  increases,  the  array 
pattern  ,ets  shown  in  Figures  5.13-5.15,  is  shifted  to  the  right  and  the  convolution 
integration  is  performed  to  produce  Vm(kz).  So,  as  the  measured  wavenumber 
neaxs  tt/Ax  the  first  negative  aliasing  lobe  neaxs  — tt/Ax  from  the  left.  Since 
V(kx)  is  almost  symmetric,  if  we  see  peaks  pjist  tt/Ax  in  the  positive  part  of  the 
spectrum  there  axe  peaks  also  in  the  negative  part  of  the  spectrum.  It  is  these 
peaks  in  the  negative  part  of  the  spectrum  that  alias  into  the  positive  part  of 
the  spectrum.  The  ideal  array  pattern  is  a  sinc(fciT/2)  function  ais  discussed  in 
Chapter  2;  the  array  pattern  for  the  1  cm  spacing  most  closely  resembles  that 
sine  function. 

For  the  comparison  of  different  array  spacings,  the  forcing  function  from 
Experiment  9  in  Table  2.2  is  used.  Because  the  drive  wavenumber  for  this 
experiment  is  the  highest  of  the  experiments  conducted  (fcj  =  1.49  cm“^),  more 
energy  should  be  seen  in  the  higher  wavenumbers  of  the  velocity  spectrum. 
With  more  energy  input  into  the  high  wavenumber  end  of  the  spectrum  than 
any  of  the  other  experiments,  the  velocity  should  be  the  most  challenging  to 
raeeisure.  The  wavenumber  dependence  at  8.0  kHz  is  exeimined.  In  Figure  5.16, 
an  array  of  37  accelerometers  equally  spaced  &oin  the  center  of  the  beam  and 
covering  the  entire  length  of  the  beam  (in  this  case  Ax  =  1.02  cm)  is  used  to 
estimate  the  wavenumber  spectrum  of  the  velocity.  As  shown  in  Figure  5.16, 
the  agreement  between  the  continuous  and  discrete  transforms  of  the  predicted 
velocity  is  excellent  except  at  the  higher  wavenumbers.  There,  the  aliasing  effects 
tend  to  cause  the  levels  calculated  by  the  FFT  to  fall  off  less  rapidly  then  the 
continuous  transform. 
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Figure  5.16  Comparison  of  the  magnitude  of  the  continuous  and 
discrete  Fourier  tramsforms  of  the  velocity  for  Exper¬ 
iment  9  at  a  frequency  of  8  kHz.  The  array  for  the 
FFT  consisted  of  37  elements  (Ax  =  1.02  cm). 


102 

Using  a  nineteen  element  array  consisting  of  accelerometers  located  at  2  cm 
intervals  equally  spaced  away  from  the  center  of  the  beam,  also  produces  good 
agreement  between  the  discrete  and  continuous  transforms  for  Experiment  9  as 
presented  in  Figure  5.17.  At  the  Nyquist  wavenumber  (defined  here  as  tt/Az) 
the  local  peak  in  the  response  causes  the  discrete  transform  to  overestimate  once 
again.  The  results  for  an  array  using  nine  accelerometers  is  shown  in  Figure  5.18. 
Considering  that  the  Nyquist  wavenumber  for  this  array,  0.75  cm”^,  is  one  half 
the  drive  wavenumber,  the  results  are  remarkably  good.  It  seems  that  the  aliasing 
effects  have  somehow  canceled  out  in  the  convolution  integral.  Generally,  so  long 
as  the  largest  amplitude  seen  in  the  velocity  spectrum  is  contained  within  the 
Nyquist  rate  of  the  array,  the  FFT  approximates  the  genereil  features  of  the 
spectrum  very  well. 

Returning  to  the  nineteen  accelerometer  array,  the  effects  of  a  slight  decrease 
in  the  inter-element  spacing  is  now  examined.  By  decreasing  Az  the  aliasing 
wavenumber  is  increztsed  slightly.  By  doing  so,  the  local  peak  in  the  true 
wavenumber  spectnun  now  falls  below  the  aliasing  lobe  of  the  array  pattern  if  the 
Az  is  chosen  correctly.  The  result  of  using  an  Eirray  of  21  accelerometers  is  shown 
in  Figure  5.19.  While  the  discrete  transform  still  overestimates  the  value  of  the 
spectrum  slightly  around  7r/2  cm”^,  the  shape  of  the  true  spectrum  is  faithfully 
followed.  It  is  noteworthy  that  merely  by  using  two  more  accelerometers  the 
results  may  be  noticeably  improved. 

The  importance  of  knowing  the  highest  wavenumber  at  which  a  peak  is 
to  be  seen  has  been  demonstrated.  While  examining  the  dispersion  relation 
for  the  beam  (see  Figure  3.2)  yields  an  idea  of  what  the  highest  freebending 
wavenumber  to  expect  for  a  pairticulair  frequency  range,  the  content  of  the  forced 
wavenumber  content  must  be  known  as  well.  At  8.0  kHz  the  forced  wavenumber 


dB  re  1  cm*/dyne-s 


103 


-80 

-90 

-100 

-110 

-120 

-130 

-140 

-ISO 

•160 

-170 


-180-+- 


1 

1 

1 

' 

i 

_ i _ 

: 

- ! - , 

1  1  1 

- 

' 

• 

’ 

_ 

.  i  : 

ai 

,  ^  i 

^ .  k 

V 

■1 

■ 

H 

nn 

9 

|H 

■ 

H - 

— 

— 

— 

— 1 

0.00  0.28  O.SO  0.7S  1.00  1.2S  1.50  1.7S  2.00  2.2S  2. SO  2.75  3.00  3.25 

Wavenumber  (l/cm) 


analytic  transform 

^fft  of  theoretical  results.  19  accelerometers 


Figure  5.17  Comparison  of  the  magnitude  of  the  continuous  and 
discrete  Fourier  tremsforms  of  the  velocity  for  Exper¬ 
iment  9  at  a  frequency  of  8  kHz.  The  array  for  the 
FFT  consisted  of  19  elements  (Ax  =  2  cm). 
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Figure  5.18  Comparison  of  the  magnitude  of  the  continuous  and 
discrete  Fourier  transforms  of  the  velocity  for  Exper¬ 
iment  9  at  a  frequency  of  8  kHz.  The  array  for  the 
FFT  consisted  of  nine  elements  (Az  =  4  cm). 
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Comparison  of  the  magnitude  of  the  continuous  and 
discrete  Fourier  transforms  of  the  velocity  for  Exper¬ 
iment  9  at  a  frequency  of  8  kHz.  The  array  for  the 
FFT  consisted  of  21  elements  (Ax  =  1,8  cm). 


Figure  5.19 
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content  is  shown  in  Figure  5.20.  This  shows  that  for  this  frequency  the  force 
wavenumber  content  has  been  successfully  steered  to  the  desired  wavenumber. 
However,  the  main-lobe  of  the  wavenumber  pattern  is  fairly  broad  and  extends 
past  the  Nyquist  wavenumber  for  a  spacing  of  2  cm.  By  realizing  the  extent 
of  the  force  wavenumber  spectrum,  the  meeisurement  array  can  be  adjusted  for 
proper  measurement. 

The  final  array  effect  to  be  ex2imined  is  that  of  truncation  error.  The  length 
of  a  discrete  transform  is  Lt, 

Lt  =  ,VAx  ,  (5.4) 

where  N  is  the  number  of  samples.  As  was  mentioned  previously,  it  is  desirable  to 
have  the  Fourier  transform  of  the  array  function  match  the  sine  function  with  the 
correct  length,  Lf  The  location  of  the  actual  measurements  is  given  in  the  very 
beginning  of  this  chapter.  The  length  of  the  FFT  is  A’Ai  =  37  cm;  the  length 
should  have  been  forced  to  be  38  cm  either  by  using  a  different  spacing  or  by  using 
a  different  starting  point  (which  would  2illow  the  use  of  another  accelerometer). 
The  effect  of  truncation  error  is  noticeable  eis  a  shifting  of  the  location  of  the 
sidelobe  structure  at  wavenumbers  approaching  the  Nyquist  wavenumber.  The 
spacing  used  in  the  experiment  is  input  to  ARRAY  JlFFECl  to  cadculate  the 
beam  response  at  8.0  kHz  for  Experiment  9  and  the  wavenumber  response  is 
obtained  via  the  FFT.  The  results  are  sliown  in  Figure  5.21.  The  shifting  of  the 
sidelobe  pattern  is  clearly  visible.  The  array  used  to  calculate  the  spectrum  in 
Figure  5.16  used  a  data  window  whose  length  is  nearly  38  cm,  and  the  results 
2ire  much  better.  The  truncation  error  causes  the  width  of  the  main-lobe  of  the 
array  pattern  to  be  wider  than  that  of  the  proper  sine  function.  The  sidelobes 
at  the  higher  wavenumbers  are  misaligned. 
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Figure  5.20  Wavenumber  content  of  the  force  relative  to  the  center 
shaker  for  Experiment  9  at  a  frequency  of  8  kHz. 
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Figure  5.21  Comparison  of  the  magnitude  of  lUe  continuous  and 
discrete  Fourier  transforms  of  the  velocity  for  Exper¬ 
iment  9  at  a  frequency  of  8  kHz.  The  array  for  the 
FFT  consisted  of  37  elements  (Ai  =  1.0  cm). 
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Although  the  error  of  1  cm  causes  only  a  small  deviation,  the  truncation 
problem  can  be  more  severe.  For  example,  for  an  accelerometer  spacing  of 
2  cm  it  is  tempting  to  place  accelerometers  only  at  seventeen  shaker  locations. 
This  spacing  yields  a  Nyquist  wavenumber  which  is  the  same  as  the  Nyquist 
wavenumber  for  the  array  that  produced  the  plot  in  Figure  5.16;  but  Lt,  the 
aperture  length,  for  this  case  is  only  34  cm.  The  results  for  this  truncated  array 
cire  shown  in  Figure  5.22.  This  result  clearly  points  out  the  danger  of  truncation 
error. 

In  this  section,  the  effects  of  the  spacing  zind  length  of  the  measurement  array 
on  the  estimated  wavenumber  spectrum  have  been  investigated.  The  spacing 
must  be  close  enough  to  enable  the  peaks  at  the  higher  wavenumbers  to  be 
measured.  A{kt)  must  be  as  close  as  possible  to  a  sinc(fciI-/2),  hence  the  length 
of  the  transform  must  be  almost  that  of  the  structure.  The  length  should  fall 
within  0.5Ax  of  the  length  of  the  beam. 

5.4.  Wavenumber-Frequency  Dependence  of  the  Force  and  Velocity 

The  main  results  of  this  thesis  are  presented  in  this  section.  The  wavenumber- 
frequency  dependence  of  the  measured  force  as  well  as  the  measured  and 
predicted  velocity  are  presented  in  surface-contour  plots.  The  surface  portion 
of  the  plots  allow  for  a  quzilitative  examination  of  the  mean  features  of  the 
spectrum.  The  contours  of  the  surface  axe  shown  in  a  pleine  beneath  the  surface. 
From  these  contours,  more  detailed  information  about  the  values  of  the  spectrum 
can  be  obtaiined.  From  these  plots,  two  of  the  most  importeint  questions  with 
respect  to  the  success  of  the  wavenumber-frequency  measurement  system  are 
answered.  The  first  is  whether  or  not  the  force  wavenumber  content  of  the  array 
of  point  forces  can  be  steered  to  a  desired  wavenumber  by  amplitude  shading  of 
the  drive  voltages  sent  to  the  shakers.  Second,  contingent  upon  steering  the  force 
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Figure  5.22  Comparison  of  the  magnitude  of  the  continuous  and 
discrete  Fourier  transforms  of  the  velocity  for  Exper¬ 
iment  9  at  a  frequency  of  8  kHz.  The  array  for  the 
FFT  consisted  of  seventeen  elements  (Ax  =  2.0  cm), 
placed  at  the  shaher  locations. 
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wavenumber  content,  what  enhancement  in  the  measurement  of  the  velocity  is 
obtained  by  exciting  the  beam  in  this  fashion.  In  other  words,  are  we  able  to 
measure  the  presence  of  the  higher  wavenumber  velocity  fields,  which  have  so  far 
been  unmeasureable. 

For  each  of  the  nine  experiments  (listed  in  Table  2.2),  the  spatial  Fourier 
transform  of  the  measured  force  relative  to  the  force  at  the  center  shaker  is 
plotted.  The  units  of  this  function  are  centimeters.  Additionally,  the  predicted 
normalized  wavenumber  velocity  spectra  (using  the  analytic  Fourier  transform 
as  formulated  in  Section  3.4)  are  presented.  As  before,  the  velocity  is  normalized 
to  the  force  at  the  center  shaker,  hence  the  Fourier  transform  has  the  units  of 
admittance  times  centimeters.  The  experimental  velocity  data  are  processed  in 
the  following  manner:  the  37  normalized  velocities  are  read  into  a  computer 
program  which  zero  pads  the  data  to  a  length  of  128  points  and  then  performs 
an  FFT.  The  general  processing  of  the  data  is  discussed  in  Chapter  2.  Since 
all  of  the  quantities  present  in  this  section  are  normalized  to  the  center  shaker's 
force,  the  designation  “normalized”  is  dropped  for  brevity.  Also,  any  peculiarities 
associated  with  any  of  the  data  sets  (e.g.,  missing  accelerometer  locations)  ran 
be  found  in  Appendix  D. 

5.4.1.  Nine  Shaker  Experiments 

For  the  first  experiment,  the  drive  wavenumber  is  Jtj  =  0  cm~'.  The 
wavenumber  dependence  of  the  force  spectra  is  shown  in  Figure  5.23.  The  force 
wavenumber  is  characterized  by  four  mmn  features,  the  forced  wavenumber  ridge, 
the  freebending  wavenumber  influence,  system  resonance  and  the  periodicity  of 
the  spectrum. 

The  forced  wavenumber  ridge  is  the  ridge  apparent  at  the  drive  wavenumber 
in  Figure  5.23.  it  is  a  high  amplitude  region  for  constant  k^.  The  forced 
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wavenumber  ridge  is  also  called  the  main-lobe  of  the  force  wavenumber  spectrum. 
A  coherent  ridge  occurs  when  the  amplitude  shading  of  the  input  voltages  has 
succeeded  in  imposing  a  similzir  distribution  on  the  forces  at  the  shahers.  The 
grating  lobes  are  due  to  the  inherent  periodicity  of  the  spectrum  of  discrete 
point  forces  (the  period  is  2Tr/Ax).  Hence,  if  the  forced  wavenumber  ridge 
has  been  successfully  enforced  there  also  are  grating  lobes  present  at  kz  = 

±27r//Ax  ±  kd,l  =  ±1,±2 .  For  the  first  experiment,  the  main-lobe  and  the 

grating  lobes  can  be  clearly  seen  at  0,  ir/2  and  ir  cm~^ . 

The  system  resonance  is  described  in  Section  5.2.  On  the  surface  plots,  the 
system  resonance  appears  as  a  sharp  ridge  for  constant  frequency.  This  is  a 
consequence  of  the  low  impedance  seen  by  the  center  shaker  eind  the  low  force 
output  due  to  this  small  impedance.  At  resonance,  all  of  the  normalized  forces 
ana  velocities  appear  high  because  they  eure  measured  with  respect  to  the  force 
at  the  center  shaker.  The  system  resonance  is  not  labeled  on  Figure  5.23,  but 
can  be  seen  around  3.0  kHz. 

The  freebending  wavenumber  for  the  Timoshenko  beeim  has  approximately 
a  square  root  dependency  on  frequency  (see  Figure  3.2)  for  frequencies  below 
12  kHz.  For  a  Bernoulli-Euler  beeun,  the  relationship  for  the  freebending 
wavenumber  is  exactly  a  square  root  dependency.  Due  to  the  high  admittance 
of  the  beam  at  the  freebending  wavenumber,  it  is  expected  that  there  is  a  region 
around  a  parabolic  line  where  the  influence  of  the  freebending  wavenumber  can 
be  seen.  The  pairabola  hzis  a  positive  kx  branch  and  a  negative  kx  branch.  The 
influence  of  the  positive  branch  on  the  spectrum  can  be  seen  in  Figure  5.23. 
Regions  where  the  effect  of  the  higher  admittance  at  the  freebending  wavenumber 
can  be  seen  are  denoted  by  the  label  ‘freebending  wavenumber’  on  the  plots  (in 
Figure  5.23).  Even  though  the  drive  voltage  is  not  steered  to  the  freebending 
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wavenumbers,  the  force  spectrum  exhibits  a  ridge  at  these  wavenumbers  because 
the  betim  response  is  high  enough  to  drive  the  shakers.  For  all  of  the  plots, 
not  every  local  ridge  associated  with  this  effect  is  pointed  out  explicitly.  The 
periodicity  of  the  spectrum  means  that  there  is  a  repetition  of  this  parabola 
centered  at  =  ±27:11  Ax,  I  =  ±1,  ±2, ... ;  these  periodic  repetitions  can  be  seen 
at  7r/2  cm~^  and  tt  cm~^  (labeled  in  Figure  5.23  as  periodic  repetitions  of  the 
freebending  wavenumber). 

The  wavenumber- frequency  velocity  spectrum  predicted  using  the  measured 
forcing  function  from  Experiment  1  is  shown  in  Figures  5.24  and  5.25.  The  first 
plot  is  the  analytic  Fourier  transform  and  the  second  is  the  FFT  of  the  predicted 
velocity.  As  these  plots  are  almost  identicad,  only  the  analytic  Fourier  transform 
is  shown  for  the  rest  of  the  experiments.  The  meaisured  wavenumber-frequency 
velocity  response  is  shown  in  Figure  5.26. 

The  wavenumber-frequency  velocity  surfaces  are  chau-acterized  by  three  maun 
features,  the  forced  wavenumber  ridge,  the  freebending  wavenumber  and  the 
system  resonance.  Note  that  there  is  no  inherent  periodicity  in  the  trauisform  of 
the  continuous  velocity  field.  The  freebending  wavenumber  ridge  is  due  to  the 
high  admittance  in  this  region  and  the  forced  wavenumber  ridge  is  due  to  the  large 
amplitude  of  the  force  wavenumber  spectrum  in  this  region.  In  Figures  5.24-5.26, 
the  forced  wavenumber  ridge  and  the  freebending  wavenumber  peadcs  can  be  seen 
clearly  as  marked  on  the  plots.  The  agreement  between  theory  and  experiment  is 
quite  good.  The  system  resonamce  can  be  seen  at  3.1  kHz  on  all  of  the  plots.  Both 
of  the  theoreticadly  produced  plots  have  a  ridge  at  the  first  grating  lobe  of  the 
force  wavenumber  spectrum,  at  —  r/2  cm~^  (labeled  as  “forced  wavenumber 
ridge”  on  the  plot).  In  order  to  show  this  pesdc,  a  slice  is  made  through  the  surface 
at  a  frequency  of  11.5  kHz.  This  plot  of  the  analytic  and  measured  wavenumber 


Wavenimiber-firquency  surface-contour  plot  of  the  i)rc(.licted  velocity  spectrum  (spec 
trum  via  the  FFT)  for  Experiment  1.  Velocity  relative  to  Ff(u>). 
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velocity  is  shown  on  Figure  5.27.  The  peaJc  can  be  seen  in  the  predicted  results. 
The  measurements  do  not  reflect  this  however,  and  it  may  be  that  the  large 
amplitude  of  the  spectrum  at  the  lower  wavenumbers  somehow  overwhelms  the 
spectrum  at  the  higher  wavenumbers.  A  possible  mechanism  for  the  burying  of 
the  higher  wavenumber  components  in  the  measured  data  is  the  following;  the 
accelerometer  is  more  sensitive  at  the  lower  wavenumbers,  the  peak  at  the  zero 
wavenumber  ajid  the  associated  sidelobes  are  given  more  weight  than  the  smaller 
peak  at  the  higher  wavenumbers.  In  this  way,  the  stronger  component  at  the 
lower  wavenumbers  buries  the  weaker  components  at  the  higher  wavenumbers  in 
the  measurement;  but  the  theory  has  no  such  sensor  error  and  is  able  to  measure 
the  higher  wavenumber  peaJc. 

In  Experiment  2,  the  drive  wavenumber  is  0.16  cm“^.  Due  to  the  width 
of  the  main-lobe  in  the  force  spectrum,  the  peaks  at  ±.kd  overlap  as  shown 
in  Figure  5.28.  The  grating  lobes  and  the  main-lobe  create  force  wavenumber 
passbzind  regions.  The  freebending  wavenumber  influence  is  not  as  noticeable 
as  in  Experiment  1,  because  the  forced  wavenumber  ridge  is  present  where  the 
freebending  wavenumber  effect  would  be  seen.  The  system  resonance  is  seen 
again  at  3.1  kHz. 

The  meaisured  and  predicted  wavenumber  velocity  responses  are  shown  in 
Figures  5.29  and  5.30.  Again,  the  forced  wavemunber  ridge  is  clearly  present 
and  theory  and  experiment  agree  nicely,  especially  below  kx  =  1.5  cm“'.  The 
grating  lobe  effect  around  kx  —  ?r/2  cm”*  can  be  seen  in  the  predicted  velocity 
spectrum  but  not  in  the  meeisured  spectrum.  This  ridge  in  the  predicted  velocity 
is  present  above  9.0  kHz.  The  system  resonance  occurs  at  3.1  kHz.  Only  slight 
local  peahs  can  be  seen  in  the  region  where  the  freebending  wavenumber  influence 
is  expected. 


forced  waveniiiiiLer 
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Wavenumber-frequency  surface-contour  plot  of  the  measured  velocity  spectrum  for 
Experiment  2.  Velocity  relative  to  Fc{u}). 
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The  force  wavenumber  spectrum  for  Experiment  3  is  displayed  in  Figure  5.31. 
The  drive  wavenumber  for  this  experiment  is  kd  =  0.332  cm“^  The  forced 
wavenumber  ridge  for  this  experiment  is  the  least  defined  for  any  of  the  nine 
shaker  experiments.  From  2-6.5  kHz,  the  forced  wavenximber  ridge  is  fairly  well 
defined  and  properly  located.  From  6.5-9  kHz,  the  forced  wavenvunber  ridge 
(and  the  grating  lobes  as  well)  have  shifted  and  become  less  defined.  A  forced 
wavenumber  ridge  once  again  appears  from  9-12  kHz.  For  this  experiment  and 
for  Experiment  7,  kd  corresponds  to  a  freebending  wavenumber  whose  frequency 
is  within  the  measurement  range  of  2-12  kHz.  The  frequency  of  a  freebending 
wavenumber  equal  to  0.332  cm~^  is  5.0  kHz.  It  is  not  necessarily  expected  that 
there  would  be  a  peak  in  either  the  force  wavenumber  content  or  the  velocity 
wavenumber  content  at  this  frequency.  This  is  not  expected  for  two  reasons. 
One,  this  is  not  a  resonance  wavenumber  of  the  beam.  Two,  even  if  it  were  a 
resonance  wavenumber,  the  force  distribution  must  match  that  particular  mode 
in  order  to  produce  a  highly  resonance  response. 

Figures  5.32  and  5.33  are  the  surface  contour  plots  showing  the  dependence 
of  the  velocity  on  wavenumber  and  frequency  for  Experiment  3.  The  velocity  in 
general  tracks  with  the  force  wavenumber  main- lobe.  As  the  main- lobe  of  the 
force  shifts  up  or  down  in  wavenumber  so  does  the  velocity  response.  Since  the 
forced  wavenumber  ridge  wanders,  the  velocity  ridge  does  not  stay  at  a  constant 
wavenumber.  The  freebending  wavenumber  effect  is  not  seen.  Perhaps,  it  is 
the  proximity  of  the  drive  wavenumber  to  the  freebending  wavenumbers  in  this 
frequency  range  that  is  the  explzmation  for  the  lack  of  a  nicely  formed  forced 
wavenumber  ridge.  The  system  resonance  for  the  measured  values  is  higher  than 
the  predicted  value.  However,  if  the  mass  loading  would  be  removed,  the  two 
curves  should  compare  more  closely. 
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Figure  5.31  Wavenumber-frequency  surface-contour  plot  of  the  measured  force  spectrum  for  Experi¬ 
ment  3.  Force  relative  to  F’c(u;). 
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The  final  nine  shaker  experiment  is  Experiment  4.  Here,  the  shakers,  zu-e 
driven  in  a  pattern  such  that  every  other  shaker  is  out  of  phase  with  the  center 
shaJcer  (i.e.  the  first  shaker  is  in  phase,  the  second  out  of  phase).  The  force 
wavenumber  content  in  Figure  5.34  displays  the  strong  main  jind  grating  lobe 
pattern  near  the  desired  wavenumber  until  about  10  kHz.  The  peak  of  the  forced 
wavenumber  spectrum  is  found  at  ki  =  ir/4  cm~^  instead  of  the  value  used  to 
generate  the  weights  in  FORCE-FILTER,  0.746  cm~^.  It  is  difficult  to  control 
the  wavenumber  spectra  to  this  fine  a  degree.  The  influence  of  the  freebending 
wavenumber  can  once  again  be  seen  and  is  labeled  on  the  plot. 

Figures  5.35  and  5.36  are  the  surface  contour  plots  for  the  predicted  and  the 
measured  velocity  spectra  for  Experiment  4.  The  forced  wavenumber  ridge  is 
seen  to  become  more  and  more  defined  as  frequency  increases.  At  frequencies 
below  3  kHz,  the  admittance  of  the  beam  at  kd  is  so  low,  that  even  with  the  well 
defined  forced  wavenumber  ridge,  the  velocity  spectrum  does  not  exhibit  a  ridge 
at  the  drive  wavenumber  .  At  these  low  frequencies,  the  freebending  wavenumber 
dominates  the  response  of  the  structure.  The  overall  spectrum  of  the  velocity 
response  is  much  lower  than  the  other  experiments  because  most  of  the  energy  of 
the  force  spectrum  falls  into  a  low  admittance  region.  The  system  resonance  is 
very  much  subdued  as  compared  to  some  of  the  other  experiments,  and  cannot 
be  seen  in  the  force  wavenumber  content.  The  grating  lobe  appears  at  such  a 
high  wavenumber  that  the  velocity  content  is  not  disturbed  at  that  wavenumber, 

5.4.2.  Seventeen  Shaker  Experiments 

The  results  for  the  five  seventeen  shaker  experiments  are  now  presented. 
The  force  wavenumber-frequency  spectra  for  the  nine  and  seventeen  shaker 
experiments  are  very  similar.  The  major  difference  between  the  two  is  due 
to  the  shaker  spacing;  the  periodicity  of  the  seventeen  shaker  force  spectra  is 


Figure  5.35  Wavenumber-frequency  surface-contour  plot  of  the  predicted  velocity  spectrum  (spec¬ 
trum  via  the  continuous  transform)  for  Experiment  4.  Velocity  relative  to  Fc(u>). 
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rr  cm“^  Hence,  the  grating  lobes  are  separated  farther  from  the  main-lobe 
and  the  periodic  repetition  of  the  freebending  wavenumber  parabola  is  centered 
further  away  from  the  origin.  The  means  by  which  the  results  from  the  first  four 
seventeen  sheiker  experiments  are  to  be  presented  is  to  point  out  any  differences 
from  the  corresponding  nine  shaker  experiment  and  any  other  noticeable  features. 

The  first  four  seventeen  shaker  experiments  are  driven  at  the  same  wavenum¬ 
ber  as  the  nine  shaker  experiments  (see  Table  2.2).  The  force  and  velocity  spectra 
for  the  matching  experiments  axe  very  similcir.  The  force,  predicted  velocity 
and  measured  velocity  wavenumber-frequency  distributions  for  Experiment  5  axe 
shown  in  Figures  5.37-5.39  respectively,  for  Experiment  6  (one  wavelength  over 
the  structure)  in  Figures  5.40-5.42,  for  Experiment  7  (two  wavelengths  over  the 
structure)  in  Figures  5.43-5.45  and  for  Experiment  8  (4.5  wavelengths  over  the 
structure)  in  Figures  5.46-5.48. 

The  major  differences  seen  between  the  nine  shaker  experiment  of  a  paxticulzu' 
force  distribution  and  the  seventeen  shaker  experiment  with  the  same  force 
distribution  axe  summarized  in  this  peiragraph.  In  Experiments  5  and  6,  there 
axe  no  grating  lobes  present  at  kg  =  itj2  cm“'  (Figures  5.37  and  5.40).  Thus, 
the  velocity  spectra  for  these  two  experiments  (see  Figures  5.38  auid  5.41)  do 
not  exhibit  a  ridge  at  this  wavenumber  as  do  the  predicted  spectra  for  the  nine 
shaker  experiments.  In  ail  of  the  seventeen  shaker  experiments,  the  grating  lobes 
of  the  force  wavenumber  distribution  occur  at  a  wavenumber  that  is  too  large 
to  noticeably  excite  the  velocity.  In  Experiment  6,  the  system  resonauice  is  very 
strongly  excited,  more  so  in  the  mezisurement  than  the  theory.  For  Experiment  7, 
the  forced  wavenumber  ridge  is  even  more  poorly  defined  than  for  Experiment  3 
as  can  be  seen  in  Figures  5.31  zmd  5.43.  The  velocity  spectrum  for  this  experiment 
looks,  however,  very  similEir  to  the  corresponding  nine  shaker  experiment.  For 
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Wavenumber-frequency  surface-contour  plot  of  the  predicted  velocity  spectrum  (spec 
trum  via  the  continuous  transform)  for  Experiment  5.  Velocity  relative  to  Fc(w). 
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Figure  5.40  Wavenumber-frequency  surface-contour  plot  of  the  measured  force  spectrum  for  Experi¬ 
ment  6.  Force  relative  to  Fc{w). 
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Figure  5.43  Wavenumber-frequency  surface-contour  plot  of  the  measured  force  spectrum  for  Experi 
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Figure  5.44  Wavenumber-frequency  surface-contour  plot  of  the  predicted  velocity  spectrum  (spec¬ 
trum  via  the  continuous  transform)  for  Experiment  7.  Velocity  relative  to  Fc(a;). 
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freebending  wavenumber  repetitions  of  the 


Wavenumber-frequency  surface-contour  plot  of  the  predicted  velocity  spectrum  (spec 
trum  via  the  continuous  transform)  for  Experiment  8.  Velocity  relative  to  Fc(u;). 


Figure  5.48  Wavenumber- frequency  surface-contour  plot  of  the  measured  velocity  spectrum  for 

Experiment  8.  Velocity  relative  to  Fc(u>). 
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the  Experiment  S,  again  the  nine  shaJcer  experiment  seems  to  have  had  more 
success  in  prescribing  a  coherent  forced  wavenumber  ridge  at  kd  as  shown  in 
Figures  5.34  and  5.46.  Especially  at  the  higher  frequencies  (above  9  kHz),  the 
ridge  for  the  seventeen  shaker  experiment  breaks  up.  The  resulting  velocity 
response  is  characterized  by  a  forced  wavenumber  ridge  that  is  less  coherent 
than  for  the  nine  shaker  experiment. 

Experiment  9  is  unique  to  the  seventeen  shaker  experiment  as  the  shakers 
axe  excited  at  the  Nyquist  rate.  The  peak  of  the  force  wavenumber  spectrum 
occurs  at  kz  =  tt/2  cm~*;  in  Figure  5.49  the  success  in  prescribing  the  drive 
wavenumber  can  be  seen.  Because  the  drive  wavenumber  corresponds  to  a 
freebending  wavenumber  at  60  kHz,  the  main-lobe  of  the  force  wavenumber 
content  does  .not  interfere  with  the  freebending  wavenumber  region.  So  the 
influence  of  the  freebending  wavenumber  along  with  the  periodic  repetition  of 
the  negative  wavenumber  branch  emanating  from  kg  =  it  cm~^  can  be  seen  on 
the  force  wavenumber  spectrum.  The  predicted  velocity  spectrum  display  the 
forced  wavenumber  ridge  beginning  at  7  kHz  (see  Figure  5.50).  The  measured 
velocity,  presented  in  Figure  5.51,  also  shows  a  ridge  at  the  forced  wavenumber, 
but  at  higher  frequencies.  A  slice  of  the  predicted  and  measured  velocity  surfaces 
is  presented  in  Figure  5.52  at  a  frequency  of  11.9  kHz.  The  drive  wavenumber 
for  this  experiment  is  the  same  £is  the  grating  lobe  for  Experiment  1.  In  the 
Experiment  1,  the  measured  wavenumber  velocity  spectrum  did  not  result  in 
any  ridge  at  the  high  wavenumber,  while  in  Experiment  9,  in  the  absence  of 
the  strong  low  wavenumber  velocity  field,  the  high  wavenumber  variations  were 
detected  experimentally. 

In  summary,  several  important  results  were  discovered.  The  ability  to  steer 
the  wavenumber  content  of  the  array  of  point  sources  merely  by  steering  the  drive 
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Figure  5.49  Wavenumber-frequency  surface-contour  plot  of  the  measured  force  spectrum  for  Experi¬ 
ment  9.  Force  relative  to  Fc{u}). 
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Figure  5.50  Wavenumber-frequency  surface- contour  plot  of  the  predicted  velocity  spectrum  (spec¬ 
trum  via  the  continuous  transform)  for  Experiment  9.  Velocity  relative  to  f’c(uj). 
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Figure  5.52  Comparison  of  the  magnitude  of  the  predicted  spec¬ 
trum  versus  the  measured  spectrum  for  Experiment  9 
at  a  frequency  of  11.9  kHz. 


voltages  was  shown.  In  general,  the  velocity  spectra  in  the  low  wavenumber 
region  was  accurately  measured.  The  accuracy  of  the  measurement  system 
(using  the  theoretical  values  for  comparison)  in  the  high  wavenumber  region 
for  this  system  (above  1.5  cm"*)  was  improved  by  using  an  array  of  shakers  to 
provide  excitation  concentrated  in  the  high  wavenumber  region.  The  ability  in 
Experiment  9  to  measure  the  forced  wavenumber  ridge  while  in  Experiment  1  the 
system  was  not  able  to  measure  velocity  components  at  the  same  wavenumber 
demonstrates  the  improved  capability. 

5.5.  Relationship  between  Force  and  Velocity  in  Wavenumber  Space 


In  this  section,  the  relationship  between  the  wavenumber-frequency  spectrum 
of  the  force  and  that  of  the  velocity  field  is  investigated.  The  finite  beam 
admittance  as  defined  in  Chapter  3  is 

V(k,,u;) 


(5.5) 


Q{kx,u) 

Both  the  normalized  force  and  the  normalized  velocity  spectra  are  normalized 
to  the  force  metisured  at  the  center  of  the  betim.  The  units  of  this  function 
axe  those  of  admittance  (centimeters  per  dyne  seconds).  The  examination  of 
this  function  centers  tiround  the  postulation  in  Chapter  3  that  the  finite  beam 
admittance  should  approach  the  infinite  beam  admittance  when  the  amplitude  of 
the  force  wavenumber  spectrum  is  high.  To  this  end,  the  frequency  variation  of 
the  predicted  and  measured  finite  beeim  admittances  axe  compzired  to  the  infinite 
beam  admittance  as  a  function  of  frequency  at  the  drive  wavenumber.  The 
FORTRAN  program  COMPARE  (see  Appendix  C),  computes  the  admittance  of 
the  various  functions  at  a  wavenumber  as  a  function  of  frequency  and  plots  the 
results.  The  wavenumber  at  which  the  results  axe  presented  corresponds  to  the 
bin  of  the  FFT  that  is  closest  to  the  drive  wavenumber.  The  notation  on  the 
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j)lots.  "oxact  transform",  indicates  that  the  finite  beam  admittance  is  calculated 
usinjf  the  analytic  Fourier  transform  solution  to  the  wavenumber  velocity. 

In  Fijrures  5.53a  and  5.53b.  the  three  admittances  for  each  of  Experiment  1 
and  5  are  displayed.  The  drive  wavenumber  is  0  cm~*.  All  six  of  the  plots  are 
virtually  identical.  The  success  in  prescribing  a  high  force  output  at  the  drive 
wavenumber  resulted  in  the  finite  beam  admittance  being  the  same  ais  the  infinite 
beam  admittance.  At  the  zero  wavenumber,  the  beam  acts  as  a  lumped  mass 
and  the  admittance  looks  masslike. 

The  admittance  for  Experiment  2  and  G  are  shown  in  Figure  5.54a  and  5.54b 
(respectively).  For  the  finite  beam  admittance  a  peak  is  seen  around  3  kHz 
for  both  the  nine  and  seventeen  shaker  experiments.  This  peak  occurs  at  a 
lower  frequency  for  the  seventeen  shaker  experiment;  otherwise  the  finite  beam 
admittance  for  nine  and  seventeen  shakers  are  the  same.  Note  also  the  excellent 
agreement  between  the  predicated  and  measured  values  for  the  finite  beam 
admittance.  Also,  the  infinite  beam  admittance  is  very  close  to  the  finite  beam 
admittance  except  around  3  kHz.  The  drive  wavenumber  for  this  experiment  is 
O.IG  cm“‘  and  the  center  of  the  nearest  bin  of  the  FFT  (for  which  the  data  are 
presented)  is  0.147  cm~*. 

The  drive  wavenumber  for  the  third  and  seventh  experiments  is  0.332  cm~’. 
A  freebending  wavenumber  of  the  same  value  occurs  around  5  kHz.  The 
admittances  are  shown  at  kz  =  0.344  cm~*  on  Figure  5.55a  and  5.55b.  The 
maximum  of  the  infinite  beam  admittance  occurs  at  the  freebending  frequency  for 
that  wavenumber.  The  theory  and  experiment  generally  agree  in  the  estimation 
of  the  admittance  for  the  finite  beam.  For  the  finite  beam  admittance,  the  first 
peak  in  the  admittance  of  the  nine  and  seventeen  shaker  experiments  takes  place 
at  nearly  the  same  frequency.  The  second  peak  of  the  finite  beam  admittance 
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Figure  5.53  Plots  of  the  ansJytic  and  experimental  finite  beam 
admittance  along  with  the  admittance  of  the  infinite 
beam  at  =  0.  cm“*  for  (a)  Experiment  1  and  (b) 
Experiment  5. 
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=  0.147  cm~^  for  (a)  Experiment  2  and 
;nt  6. 
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(b) 

Figure  5.55  Plots  of  the  einalytic  and  experimental  finite  beam 
admittance  along  with  the  admittamce  of  the  infinite 
beam  at  =  0.344  cm“*  for  (a)  Experiment  3  and 
(b)  Experiment  7. 
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for  Experiment  7  occurs  at  a  lower  frequency  than  does  the  second  peak  for 
Experiment  3.  Above  the  second  peak,  the  admittcinces  of  the  finite  beam  for 
both  experiments  collapses  to  the  infinite  beam  admittance.  Recall  that  the 
force  wavenumber  spectrum  for  these  experiments  was  the  least  defined  and  this 
results  in  part  of  the  discrepancy.  However,  further  investigation  is  required  to 
determine  what  effect  the  proximity  of  the  freebending  wavenumber  htis  on  the 
finite  beam  admittance. 

For  the  nine  and  seventeen  shaker  experiments  whose  drive  wavenumber 
corresponds  to  4.5  wavelengths  across  the  length  of  the  beam,  the  finite  beam 
admittance  generally  stays  close  to  that  of  the  infinite  beam;  as  pictured  in  Figure 
5.56a  and  5.56b.  Upon  examination  of  the  force  wavenumber  spectra  for  botti 
experiments  (Figures  5.34  and  5.46),  the  main-lobe  for  the  nine  shaker  spectra 
is  better  defined  and  this  is  probably  why  the  finite  beam  admittance  for  this 
experiment  stays  closer  to  the  infinite  beam  admittance. 

The  finite  beam  admittance  for  Experiment  9  is  shown  in  Figure  5.57.  There 
are  some  fairly  large  deviations  between  the  predicted  and  measured  values  of  the 
finite  beam  admittance.  This  is  expected  as  the  wavenumber  here  is  1.57  cm”* 
and  at  these  wavenumbers  the  meeisurements  are  more  difficult.  Note  that  the 
experimental  results  generally  lie  below  those  for  the  predicted  admittance;  this 
may  be  attributed  to  the  decrease  in  the  sensitivity  of  the  accelerometer  at  the 
higher  wavenumbers.  The  finite  beam  admittances  still  follow  the  trend  of  the 
infinite  beam  admittance. 

For  completeness,  similar  plots  for  Experiment  2  at  the  wavenumber  of 
0.69  cm”*  where  the  magnitude  of  the  force  wavenumber  spectrum  is  low  are 
shown  in  Figure  5.5S.  The  finite  beam  admittances  show  no  similarity  to  the 
infinite  beam  admittance  except  perhaps  at  the  higher  frequencies. 


dB  re  1  cm/dyne-8 
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oxact  tranaform  Frequency  (kHz) 

°^fft  of  experimental  results 
^admittance  of  Infinite  beam 


(b) 

Figure  5.56  Plots  of  the  analytic  and  experimental  finite  beam 
admittance  along  with  the  admittance  of  the  infinite 
beam  at  =  0.74  cm~*  for  (a)  Experiment  4  and  (b) 
Experiment  8. 


Frequency  (kHz) 


exact  transform 
o^fft  of  experimental  results 
^admittance  of  infinite  beam 


Figure  5.57  Plot  of  the  analytic  and  experimental  finite  beam 
admittance  along  with  the  admittance  of  the  infinite 
beam  at  kx  =  1.57  cm~^  for  Experiment  9. 


dB  re  1  cm/dyne- 
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exact  transform 
'^'^fft  of  experimental  results 
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Figure  5.58  Plot  of  the  analytic  and  experimental  finite  beam 
admittance  along  with  the  admittance  of  the  infinite 
beam  at  =  0.69  cm"^  for  Experiment  2. 
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The  postulation  has  been  made  that  the  finite  bezun  admittance  should 
approach  the  infinite  beam  admittance  at  a  wavenumber  where  the  force 
amplitude  is  high.  The  analysis  and  results  presented  in  this  section  certainly 
point  in  this  direction. 
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CHAPTER  6 

CONCLUSIONS  AND  RECOMMENDATIONS 
6.1.  Conclusions 

In  this  thesis,  the  experimental  and  theoretical  aspects  of  exciting  and 
measuring  the  wavenumber-frequency  spectrum  of  a  finite  structure  have  been 
examined.  A  system  that  steers  the  force  amplitude  of  an  array  of  shakers 
to  a  particular  wavenumber  have  been  developed.  The  response  of  a  free- 
free  beam  to  this  excitation  have  been  measured  experimentally  and  predicted 
theoretically.  Experiment  and  theory  have  been  shown  to  match  very  well,  both 
in  point  transfer  admittances  and  in  the  comparison  of  wavenumber  spectra.  The 
relationship  between  the  wavenumber  dependence  of  the  force  and  the  velocity 
have  been  examined  and  caste  in  terms  of  the  finite  beam  admittance.  The 
finite  beam  admittance  have  been  shown  experimentally  to  approximate  the 
infinite  beam  admittance  in  regions  of  wavenumber-frequency  space  where  the 
force  amplitude  is  high.  Also,  in  these  regions  the  finite  beam  admittance  for 
experiments  using  nine  shakers  and  seventeen  are  nearly  the  same. 

In  order  to  measure  the  relation  between  the  force  and  the  velocity,  one  must 
be  reasonably  sure  that  it  is  the  actued  force  and  velocity  spectra  that  are  being 
measured.  Hence,  it  have  been  a  major  emphasis  of  this  study  to  understand 
the  means  by  which  the  two  quantities  has  been  measured.  To  this  end,  the 
wavenumber  sensitivity  of  an  array  of  sensors  has  been  developed.  The  effects 
of  using  different  types  of  measurement  arrays  has  been  examined.  Equivalent 
circuit  modeling  has  been  used  both  to  predict  the  overall  response  and  the 
output  of  the  shaker /force  gage.  The  predicted  force  output  and  force  gage 
sensitivity  have  been  shown  to  be  commensurate  with  the  measured  values.  These 
predictions  have  been  a  used  as  a  check  on  the  validity  of  the  measurements  and 
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the  calibration  procedure. 

G.2.  Recommendations 

Further  theoretical  and  experimental  study  into  the  relationship  between 
the  wavenumber-frequency  content  of  the  force  and  the  velocity  spectrum  is 
indicated.  A  more  complete  analysis  of  the  results  of  Section  3.6  would  yield 
a  theoretical  model  for  analysis  of  any  possible  exceptions  to  the  empirical 
observations  made  relating  the  force  to  the  velocity.  Also,  further  examination 
of  the  data  obtained  in  this  thesis  would  aide  in  this  analysis. 

This  method  of  force  excitation  and  velocity  meeisurement  could  be  applied 
to  structures  coated  with  compliant  layers.  By  using  the  measurement  system, 
the  filtering  action  of  these  coatings  might  be  pairameterized.  Hence,  this  system 
could  be  a  useful  tool  in  the  better  design  of  compliant  layers. 

Finally,  the  use  of  other  spectral  estimation  methods  in  examining  the 
response  of  the  structure  would  be  useful  [26].  By  using  so-ceJled  high  resolution 
methods  (e.g.,  meiximum  likelihood,  maximum  entropy  and  the  extended  Prony 
methods),  the  individual  mode  shapes  might  be  resolved.  With  this  information, 
the  relation  between  the  force  and  velocity  would  become  more  straightforward. 
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APPENDIX  A 

PVF-2  FORCE  GAGE  DESIGN 

A  schematic  of  the  force  gage  design  using  PVF-2  as  the  piezo-active 
component  is  shown  in  Figure  A.l.  This  design  proved  to  be  mechzmically 
unsuitable.  Upon  repeated  attachment  cind  removal  from  the  structure,  the  force 
gage  would  delaminate.  This  problem  motivated  the  use  of  the  PZT-4  force  gage 
which  is  more  sturdy. 

The  components  and  their  dimensions  are  shown  in  Table  A.l.  The  force 
gage  is  constructed  in  the  following  manner.  First,  two  pieces  of  PVF-2  eire 
sandwiched  around  the  inner  brass  shimstock  (which  is  the  positive  lead).  The 
surface  of  positive  polarity  of  each  of  the  pieces  of  piezo- film  are  glued  facing  each 
other  toward  the  inner  brass  using  Devcon  five- minute  epoxy.  Next,  the  longer 
and  wider  strip  of  brass  shimstock  is  folded  in  half.  The  sandwiched  combination 
is  inserted  into  the  the  now  folded  outer  brass  shield.  At  this  point,  the  area 
over  the  PVF-2  is  glued  to  the  outer  brass  shield  and  held  in  a  vise  while  a  bend 
is  made  at  the  point  1.27  cm  from  the  fold  point  of  the  outer  breiss.  This  bend 
allows  for  the  force  gage  to  fit  nicely  over  the  headmass  of  the  sheiker.  After  this 
gluing  is  complete,  the  center  conductor  of  a  low  noise  coaxial  cable  is  soldered 
to  the  inner  brass  shimstock.  This  is  the  positive  electric  terminal.  The  exposed 
positive  connection  (the  combination  inner  braiss  shimstock  and  coaxial  cable) 
is  covered  by  flexible  insulating  material.  The  outer  brsiss  is  then  glued  to  the 
insulating  material.  Finally,  the  shield  wire  of  the  coaxi^Ll  cable,  which  has  been 
dressed  back,  is  soldered  or  glued  using  conductive  epoxy  to  the  outer  brass 
shimstock. 
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LOW  NOISE  COAXIAL  INSIJIATIMH  MATTPIAl 
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matirial 


•  TRANSDUCER  THICKNESS  OVER  ACTIVE 
AREA  S  0.010  INCH 


Figure  A.l  Schematic  of  the  PVF-2  force  gage  design. 
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Table  A.l  Materials  used  in  PVF-2  transducer. 


Description  of  Section 

Material 

Dimensions  (cm) 
length  X  width  x  thickness 

Outer  shield 

Brass  shim 

3.49  X  .4  X  .00254 

Force  cell 

PVF-2 

1.43  X  .35  X  .0052 

Positive  lead 

Brass  shim 

1.75  X  .3  X  .00254 

Flexible  insulation 

Kapton 

cut  tc  size 

The  PVF-2  may  be  cleaned  using  a  Freon  bzised  solvent.  The  other 
components  of  the  force  gage  are  cleaned  before  gluing  using  the  procedures 
mentioned  in  Chapter  2. 

The  force  gage  is  glued  to  the  top  of  the  headmass.  The  force  gage/shaJcer 
is  then  glued  to  the  test  structure.  Initially,  the  force  gage  sensitivity  (measured 
using  the  same  procedure  as  in  Section  2.3.4)  amplitude  and  phase  was  very 
flat  as  a  function  of  frequency.  However,  as  the  force  gage  was  removed  and 
re-attached  to  the  structure,  the  sensitivity  was  altered.  It  was  obvious  upon 
inspection,  the  the  force  gage  layers  were  becoming  unglued  aind  the  force  gage 
was  rendered  useless. 


APPENDIX  B 

MASS  LOADING  OF  THE  FORCE  GAGE 
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In  this  appendix,  the  manner  by  which  the  effect  of  the  mass  below  the  force 
gage  may  be  removed  from  the  measured  data  is  presented.  There  are  three 
quantities  which  require  the  mass  loading  effect  to  be  removed,  the  drive  point 
admittance,  transfer  admittance  and  the  normalized  force  measurements.  The 
mass  loading  effect  causes  an  error  in  the  force  measured  by  the  force  gage,  as 
described  in  Chapter  4.  Essentially,  the  force  applied  to  the  structure  differs 
from  that  mecisured  by  the  force  gage  by  the  impedance  of  the  mass  loading 
times  the  velocity  at  that  point.  The  results  are  presented  in  the  same  order  as 
the  cdgorithm  used  to  correct  the  errors  must  be  written.  This  is  because  the 
first  result  is  used  to  correct  the  second  mass  loaded  quantity  and  the  second  to 
correct  the  third. 

The  first  quantity  examined  is  the  measured  drive  point  admittance  at  the 
center  of  the  beam,  denoted  as  Ac.  Note  that  all  quantities  eire  understood  to  be 
functions  of  frequency,  u.  The  mezisured  drive  point  admittance  in  the  center  is. 


Ac  = 


{B.l) 


Fc  +  ZlVc 

where  Zi  =  Fc  is  the  force  applied  to  the  structure  at  the  center  of  the 

bezun,  and  Vc  is  the  velocity  measured  at  the  center  of  the  beam.  The  quantity 
Ml  is  known  from  the  analysis  of  Section  4.4.  The  denominator  of  Equation  (B.l) 
is  the  measured  force.  The  desired  quantity  is  Ac  the  velocity  divided  by  the  force 
applied  to  the  structure.  Ac  is  related  to  the  measured  quantity  by, 

Ae 


.4.  =  ^^  = 


(B.2) 


1  —  ZlAc 

From  this  equation  it  can  be  seen  that  if  the  product  of  the  measured  admittance 
of  the  beam  and  the  impedamce  of  the  mass  loading  is  much  less  thsui  one.  the 
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measured  amd  the  actual  admittance  are  the  seune.  This  product  will  generally 
be  much  lower  than  one  for  the  beam  and  the  shaher  used  in  tnis  thesis. 

The  second  quantity  of  interest  is  the  measured  transfer  admittance,  Ai  or. 


Ai  = 


Vi 


iB.3) 


Fc  +  Zlvc 

where  u,  is  the  velocity  at  the  ith  location.  The  actual  transfer  admittance  ,  .4^ 


IS 


(BA) 


where  both  the  measured  and  corrected  drive  point  admittances  have  been  used 
in  the  correction  of  the  transfer  admittance. 

The  final  quantity  to  be  corrected  is  the  normalized  force  mezisured  at  each 
of  the  shakers,  , 


Ft  +  ZiVi 


(B.5) 


Fc  +  ZlVc 

The  desired  quantity  is  Ft/Fc  or  F-^  which  is  obtained  using  the  previous 
correction  for  /!,  .  The  corrected  normalized  force  is 


Ff^  =  [F^/Ac-ZL]At  , 


(5.6) 


which  may  be  verified  by  substitution. 

In  Figure  B.l,  the  difference  between  the  corrected  and  uncorrected  drive 
point  admitteince  is  shown  for  Experiment  3.  This  experiment  was  chosen  because 
the  admittance  of  the  beam  rises  above  the  admittance  of  the  meiss  loading.  This 
means  that  the  product  of  the  beam  admittance  and  the  mass  impedzince  will  be 
around  one.  The  plot  shows  that  only  when  besun  admittance  is  above  the  mass 
loading  line  does  the  corrected  admitteince  \'ary  from  the  measured.  However, 
the  general  trend  of  the  data  is  not  significantly  altered  for  the  purposes  of  this 
study.  Hence,  the  mass  correction  algorithm  has  not  been  implemented  to  correct 
the  entir'*  data  set. 


167 


APPENDIX  C 

CITED  COMPUTER  CODE 

In  this  appendix,  the  computer  code  cited  in  the  text  of  the  thesis  is  listed 
in  order  of  citation.  The  first  program  listed  is  FORCE-FILTER. 

1  PROGRAM  Force.Filter 

2  c...  This  program  caculates  the  seighting  coellicients  lor  the 

3  c...  drive  voltages. 

4  CHARACTER  *20  name 

5  REAL  PLT(512) .PHS(512) 

6  COMPLEX  W(512)  .etapdOO) 

7  COMPLEX  WW(256) ,SUM.c_i 

s  c_i  =  (0.0, 1.0) 

9  pi  =  3.1415926S4 

10  beam.length  =  37.9 

11  lUM.WAV  =  128 

12  c...  input  location  lor  the  number  ol  shakers. 

13  12121  urite(6,*)  ’  input  number  ol  points  in  lilter' 

14  READ(5,*)num 

15  c...  linding  spatial  locations  lor  the  shakers  in  real  non-dimensional 

16  c..  space  (length  ol  the  bar  is  one.) 

17  c . .  this  method  places  one  shaker  in  the  middle  and  symmetrically  out 

18  c . .  Irom  this  shaker  locates  the  rest  ol  the  shakers 

19  il(num.eq.9)then 

20  do  i  =  l.num 

21  etap(i)  =  (2.9S  +  lloat(i-l)*4. )/beam_length 

22  enddo 

23  elseil  (num.eq.17)  then 

24  do  i  =  l.nuffl 

25  etap(i)  =  (2.95  +  lloat(i-l)*2. )/beam_length 

26  enddo 

27  else 

28  write (6,*)  '  input  either  9  or  17  shakers’ 

29  goto  12121 

30  ENDIF 

31  C.  .  . 

32  c...  lorcing  window  options  are: 

33  c...  1)  rectangular  window 

34  c...  2)  Kaiser-Bessil  window 

35  c. . .  3)  Hannning  window 

36  c...  4)  Sell  input  ol  window  coelicients 

option  held  in  variable  i.llag 


37  c. . 
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38  C  .  .  . 

39 

write(6,*)  ’  input  rect  (1),  kaiser  (2) 

or  hanningO)  self(4) 

40 

read  (5,*)  i_llag 

41  C.  .  . 

42  C.  .  . 

rectanguleir  uindos 

43 

il  (i_flag  .  eq  .  1)  then 

44 

DO  I  =  l.KUM 

45 

W(I)  =  cmplxd  .0,0. ) 

46 

ENDDO 

47  C.  .  . 

Kaiser-Besael  sindos 

48 

elseif (i.llag  .  eq  .  2)  then 

49 

vrite(6.e)  ‘  input  alpha’ 

SO 

read(S.*)  alpha 

51 

CALL  KAISERCU, MUM. ALPHA) 

52  C.  .  . 

Hannning  vindov 

53 

elseif (i.llag  .  eq  .  3)  then 

54 

DO  I  =  l.KUM 

55 

arg  =  etap(i)*2.*pi 

56 

M(i)  =.5e(l-  cos(2irg)) 

57 

EHDDO 

58  C  .  .  . 

salt  input  option 

59 

al8e].f(i_llag.eq.4)  then 

60 

do  i  «  l.nun 

61 

vrite(6, 1111)1 

62 

raad(S,4>)u(i) 

63  dd 

write(6,1112)i 

64  dd 

read(6,*)etap(i) 

65 

enddo 

66 

endif 

67  1111 

lonnatC’  input  lorca  number  ’,12) 

68  1112 

formate  Input  location  number  ’,12) 

69  C  .  .  . 

70  C  .  .  . 

In  this  section,  the  wavenumber  content  Is 

shifted  using  a 

71  C.  .  . 

symmetrical  cosine  weighting.  The  number  of  waves  across 

72  C  .  .  . 

the  length  of  the  beam  Is  used  as  the  means 

for  determining 

73  C .  .  . 

the  drive  wavenumber. 

74 

write (6,*)  '  Input  cos  weight (1)  or  not’ 

75 

read  (5,*)  l.cos 

76 

If  (i.cos.eq.l)  then 

77 

wrlte(6,*)  '  Input  desired  waves 

over  the  bar(k/2pi)’ 

78 

read  (5,e)  x_k 

79 

do  1  s  i.num 

80 

arg  =  (etap(i)-.5)*2.*x_k*pi 

81 

W(i)=  WCDecosCarg) 

82 

write  (6,*)w(I) 

83 

84 

85  C .  .  .  . 

86  C  .  .  .  . 

87  C  .  .  .  . 

88 

89  1512 

90 

91 

92 

93  1513 

94 

95 

96 

97  C.  .  . 

98  C .  .  . 

99  C .  .  , 
100 

101 

103 

103 

104 

105 

106 

107  C.  .  . 

108  C  .  .  . 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 

130 

131 

133  C.  .  . 

123  1222 

124  10 

125  C.  .  . 

126  C.  .  . 


enddo 

endil 

output  files 
type  1512 

fonnat(’  input  file  for  force  mag  uave-freq  plot  file — >',$) 
read(5, 10)name 

open(unit=ll .f ile  =  name  ,  status  =  'nes') 
type  1513 

format(’  input  file  for  force  phase  uave-freq  plot  file — >’,$) 
read(S,10)name 

open(unit=12.file  =  name  ,  status  =  'new') 
savmm  &  -1000. 

This  section  performs  the  Fourier  transform  as  per  Equation  2.27. 

do  k  =  -num_wav,  num^uav 

sum  =  cmplzCO. ,0. ) 
x.wav.num  =  float ((k))/2. 
do  j  s  l,num 

sum  *  sum  ♦  w(j)*cexp(-x_sav_aum* 

>  e_i*etap(j)) 

enddo  !  enddo  for  j 

set  up  for  plotting  routines.  Hormalize  data  to  maximum  point 

temp  s  cabs (sum) 
if  (temp. eq. 0.0)  then 

temp  =  .0000001 

sum  =  cmplx(0. 000001. 0.0000001) 
urite(6,*)  '  zero' 
endif 

plt(k)  =  20*logl0(temp) 
phs(k)  =  atan2d(aimag(sum) .real(sum)) 
if  (plt(k) .gt.vavmm)  then 
vavmm=plt (k) 
i_rem  =  k 

endif 

enddo  !  enddo  for  k 

format (20a4) 
format(a20) 

erite  transformed  data  to  data  files. 


127  c. . . 


170 

128 

DO  I  =  -nuiii_Bav,iiuin_Bav 

129 

plt(i)  =  plt(i)  “  vavDin 

130 

plis(i)  =  phs(i)  -  ph8(i_rem) 

131 

it  (phs(i)  .  It  .  -180.)  phs(i)  =  phsCi)  +360. 

132 

XX  =  lloat(i)/(2.*bea]n_length) 

133 

writeCll ,*)xx,plt(i) 

134 

urit«(12,*)xx,ph8(i) 

135 

ENDDO 

136  C.  .  . 

137  C  .  .  . 

138 

close(unit=ll) 

139 

close(unit=10) 

140 

OPEH  (UHIT=10.FILE='B«ight.DAT* ,STATUS=*HEU') 

141 

OPEH  (UHIT=1 1 . FILES 'wphase . DAT ’ , STATUS= ' HEW ' ) 

142  C.  .  . 

write  spatial  weighting  coefficients 

143 

DO  I  =  l.HUM 

144 

WRITEC 10 ,*)f loat(I)-l . . cab8(W(I) ) 

145 

WRITEdl  ,*)float(I)-l.  ,atan2d(ainag(W(I))  ,real(W(I) )) 

146 

EIDDO 

147 

CL0SE(UHIT=11) 

148 

CLQSE(UHIT=10) 

149 

STOP 

150 

EID 

151  C 

152  C 

153  C 

154 

SUBROUTIHE  KAISER  (W, I, ALPHA) 

155  C.  .  . 

Kaiser-Bessel  weights.  Algorithia  from  Carter  Ackerman 

156 

COMPLEX  V(512) 

157 

F  =  (FL0AT(I)+l.)/2. 

158 

PIA  s  3 . 1415927*ALPHA 

159 

BD  s  AIHO(PIA) 

160 

DO  1  I  »  1.1 

161 

BI  =  AlIO(PIA*SqRT(l.-(FLOAT(I)/F-l.)**2)) 

162  1 

W(I)=CMPLX(BH/BD,0.0) 

163 

RETURH 

164 

EID 

165  C.  . 

166 

FUHCTIOH  AIIOCX) 

167  C 

X  BETVEEI  0.  AID  20. 

168 

Y  *  X/2. 

169 

T  =  l.E-8 

170 

E  =  1. 

171 

DE  s  1. 

172 

DO  1  I  *  1.2S 

173 

DE  =  DE*Y/FLQAT(I) 

174 

SDE  =  DE**2 

175 

E  =  E+SDE 

176 

IF(E*T-SDE)1.1,2 

177  1 

CONTINUE 

178  2 

AINO  =  E 

179 

RETURN 

180 

END 

The  next  program  listed  is  VALTIM,  the  program  that  compares 
measured  to  the  predicted  admittance. 

1  program  valtim 

2  c - - 

3  c  This  program  calculates  the  adittance  and  phase  oi  a 

4  c  Timoshenko  bar  in  cm-dyne~sec  units  and  compares  it  to 

s  c  measured  values.  The  measured  data  is  in  g's  and  pounds-lorce 

6  c  and  must  be  converted  to  cm-dyne-sec. 

7  c  This  program  is  a  combination  of  the  vork  of  P.K.  Butto 
s  c  and  K.  Grosh. 

9  c  This  specific  program  is  for  the  first  three  nine  shaker 
10  c  experiments . 

13  REAL  E.G.RHO.IN.AREA.AX.OBE.TVO.TVEI.LN.Al.FSV 

13  REAL  PI,ZERO,BBT,BAL.FREq 

14  REAL  0HG(0:40),PHI(0:40),B(0-.40).AL(0:40),CETA 

15  REAL  BT(0;40) ,DEL(0:40),CHE(0:40) .LANDA(0;40) 

16  REAL  DF,MB,q(2).qy(2) 

17  REAL  ETAP(10}.BBTE.BALE.FRQ(450) 

le  REAL  YR(0:40.20) 

19  REAL  Y(201,2),XR,YI.PHS{201,2),fa 

20  REAL  0RT(0:40).PEIS(2),PITS(2) 

31  IITEGER  I.J.N.I.Z.ind 

22  COMPLEXES  DISP(2S0),CY, SUM, DI0I(250. 0:40) ,K0(10),SSUM 

23  C0MPLEXS8  FC(250.20).SUMS 

24  CBARACTERS20  BANE, 11. 12 

25  CHARACTER*!  P.ff.esc 

26  c  .  .  . 

27  C - 

38  c  The  values  used  are: 

39  c  E  (elastic  modulus)  -  210E10  dynes/cm*2 

30  c  G  (shear  modulus)  -  77EtO  dynes/cm~2 

31  c  IN  (area  moment  of  inertia)  =  bb‘3/12 

33  c  vd  (sidth  of  beam)  *  1.2T  cm. 
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33  c  ht  (height  of  beam)  s  1.92  cm. 

34  c  RHO  (density  of  material)  =7.8  grams/cm'S 

35  c  LH  (length  of  beam)  =37.9.  cm. 

36  c  AK  (cross  sec.  constant)  =.833 

37  c  DF  (damping  factor)  =.05 

3g  c - 

39 

40  DATA  E,G,0HE,RH0,LM,wd.ht  /  210E10  ,  77E10  , 

41  11.,  7.8  .  37.9  ,1.27.1.91/ 

42  DATA  ZERO, AK, TVER .TWO  /  0.0  .  .833  ,  19.  ,  2.  / 

43  PI=4eATAH(0HE) 

44  df  =  .05 

45  PEIS(1)=1.0 

46  PEHS(2)=2.0 

47 

48  PNTS(1)=201. 

49  PHTS(2)=201. 

50  area  =  adeht 

51  im  =  Hd*ht**3/12. 

53  MB=(RH0eAREA4>LM) 

53  P=*\' 

55  C  There  are  nine  noint  forces: 

56  C 

57  M=9 

58  C 

60  C 

61  c  reading  in  the  eigen-frequencies 

62  c 

63  C 

65  OPEI  (  UIIT=10  .  FILE=’ Cgrosh.kgl.barlbsfreq.DAT’  ,TYPE=’0LD'  ) 

66  DO  1=2,21 

67  READ  (10.*)  QNG(I) 

68  0NG(I)=0MG(I)*2*PI 

69  EKD  DO 

70  CLOSE  (  UIIT=10  ) 

71  0MG(0)=0.0 

72  0MG(1)=0.0 

73  c .  .  . 

74  c . . .  Reading  in  the  orthoganality  constants 

75  c . . . 

76  OPEI ( UIIT= 1 0 , FILE= ' [grosh . kg 1 . bar] bsort . dat ' , STATUS= ' OLD ' ) 

77  DO  1=2,21 
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78  READClO,*)  ORT(I) 

79  END  DO 

80  CLOSE (UNIT= 10) 

81  0RT(0)=1.0 

82  0RT(1)=(LM**2)/12 

83  c. . .  reading  in  forcing  function  data  —  this  is  assumed  to  be  a  constant 

84  c . . .  no  matter  where  the  accelerometer  is 

85  c -  opening  file  of  force  data  - 

86  type  1511 

87  1511  formate  input  force  file — >',$) 

88  read(5, 10)name 

89  OPEN ( UNIT= 1 0 , FILE=name , TYPE= ' OLD ' ) 

90  DO  J  =  l.M 

91  IF(J.NE.(M+l)/2)  THEN 

92  DO  I=1.PNTS(1) 

93  READdO.iaOO)  FRQ(I),FA.PHASE 

94  FA  =  10** (FA/20.) 

95  FCd.J)  =  FA*CMPLX(COSD(PHASE),SIND(PHASE)) 

96  END  DO 

97  ELSEIF(J.Eq.(M+l)/2)THEN 

98  DO  I=1,PNTS(1) 

99  FCd.J)  *  CMPLXd. 0,0.0) 

100  ENDDO 

101  endif 

102  END  DO 

103  CL0SE(UIIT=10) 

104  1300  F0RMAT(1X,FB.2,F12.8,F7.2) 

105  c - closing  file  of  force  data - 

106  C. . 

107  R*SQRTdM/(AREA*LM**2)) 

108  S*SqRT(E*IN/(AK*AREA«G*LH**2)) 

no  DO  1=2,21 

111  Bd)=SqRT(RH0*AREA*0MGd)**2*LM**4/(E*IM)) 

112  ALd)  =  (((R**2-S**2)**2+4/Bd)**2)**.6-R**2-S**2)**.5 

113  1  /SqRT(TWO) 

114  BTd)  =  <((R**2-S**2)**2+4/Bd)**2)**.6+R**2+S**2)**.S 

115  1  /SqRT(TWO) 

116  LAMDAd)=ALd)/BTd) 

117  CHEd)  =  (BTd)**2-S**2)/(ALd)**2+S**2) 

118  BAL=Bd)*ALd) 

119  BBT=Bd)*BTd) 

120  IF  (BAL.GT.TWEN)  GOTO  170 

121  DELd)  =  (COSH(BAL)-COS(BBT))/(LAMDAd)*SIHH(BAL)- 

122  1  CHEd)*SII(BBT)) 
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123  GOTO  171 

124  170  DEL(I)=1/LAMDA(I) 

125  171  COHTINUE 

126  END  DO 

127  C - 

128  C  CALCULATIONS  DUE  TO  THE  FORCING  FUNCTION. 

129  C 

130  do  i  =  l,m 

131  etap(i)  =(2.95  +  (i-l)*4.  )/lni 

132  enddo 

134  DO  J=1,M 

135  YR(0,J)=1.0 

J36  YR(l,J)=ETAP(j)-.S 

137  DO  1=2,21 

138  BBT=B(I)*BT(I) 

139  BAL=B(I)*AL(I) 

140  BBTE=BBT*ETAP(J) 

141  BALE=BAL=ETAP(J) 

142  IF  (BALE.GT.TVEN)  GOTO  200 

143  YR(I,J)=(COSH(BALE)+COS(BBTE)/CHE(I)-LAMDA(I)*DEL(I) 

144  1  *SINH(BALE)-DEL(I)*SIN(BBTE)) 

145  GOTO  201 

146  200  YR(I,J)=(COS(BBTE)/CHE(I)*DEL(I)*SIH(BBTE)) 

147  201  CONTINUE 

148  END  DO 

149  END  DO 

150  c . . .  loop  point 

151  c. . . 

153  C  SET  FOR  READING  II  DATA 

154 

155  VRITE(5,*)  'VHAT  IS  TEE  NAME  OF  THE  DATA  FILE?' 

156  READ(6.10)  NAME 

157  10  FORMAT  (A20) 

158  OPEN (UIIT=2 1 . FILE=nan« , STATUS= ' OLD ' ) 

159  DO  I=1,PITS(2) 

160  READ  (21.*)ii,  Irq(i) ,Y(i,2) .PHS(I,2) 

161  ph8(i,2)=ph8(i,2)-180. 

162  if  (phs(i,2).gt.l80.)  phs(I.2)=ph8(i,2}-360 

163  if  (ph8(I,2).lt.-180.)  ph8(I,2)=ph8(i,2)+360 

164  EID  DO 

165  CL0SE(UNIT=21) 

166  naB«(20:20}='$* 

TYPE  1313 


167 


175 


168  1313  format  ('  input  ceta — >’,$) 

169  READ(5,*)  ceta 

170  C - 

171  C  CALCULATIBG  EICEHFUNCTION  FOR  A  PARTICULAR  POIHT  OH  THE  BEAM 

172  C  FOR  THE  FIRST  21  MODES 

173  C 

174  C 

175  C 

176  C - 

177  C - 

178  C 

179  C  CALCULATIQH  OF  THE  EIGENFUNCTIONS  FOR  DIFFERENT  FREQUENCIES 

180  C  AT  A  SINGLE  POINT  "CETA".  THE  EIGENFREQUENCIES  ARE  FOR  A  FREE-FREE 

181  C  BEAM  OF  THE  TIMOSHENKO  TYPE  TAKING  INTO  ACCOUNT  SHEAR  ARE  ROTARY 

182  C  INERTIA  EFFECTS 

183  C 

184  C- 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197  C 

198  C 

199  C 
300 
201 
202  C 


DO  Z=1,PHTS(1) 

FREQ=FRq (Z) *1000 . *2*PI 
SSUM=CNPLX(ZERO .ZERO) 

DO  1=0.21 

DI0I(Z.I)=1/(MB*0RT(I)*(0MG(I)**2*CMPLX(0NE.DF)-FREQ**2)) 
SUMS=CMPLX(0. 0.0.0) 

DO  J=1.M 

SUMS=SUMS+FC(Z . J)*YR(I . J) 

END  DO 

SSUM=SSUM+DIOI (Z . I ) *SUMS*PHI ( I ) 


203 

304 

205 

306 

207 

208 
309 
210 

311 

312 


DO  1=2,21 
BBT=B(I)*BT(I) 

BAL=B(I)*AL(I) 

SALE=BAL*CETA 
BBTE=BBT*CETA 
IF  (BALE.GT.TWEN)  GOTO  150 

PHI ( I ) =COSH (B ALE) -LAMDA (I ) *DEL (I ) *SINH (B ALE) + 
1  COS(BBTE)/CHE(I)-DEL(I)*SIH(BBTE) 

GOTO  ISl 

160  PHI(I)=COS(BBTE)/CHE(I)-DEL(I)*SIN(BBTE) 

151  CONTINUE 
END  DO 


THE  RIGID  BODY  NODES.  BOTH  ROTARY  AND  TRANSLATIONAL 

PHI(0)=1.0 

PHI(l)».6*(CETA-.5) 
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213  EHD  DO 

214  cy  =  SSUM*cmplx(0.0,lreq) 

215  y(2,l)  =  20*logl0{cabs{cy)) 

216  phs(z,l)  =  atan2d(aimag(cy)  .reziKcy) ) 

217  END  DO 

219  C 

220  C 

221  C  PLOTTING  ROUTINE  GOES  HERE 

222  C 

223  C 

224  C 

226  C 

227  STOP 

228  END 

229 

230 

The  source  code  for  the  FORTRAN  program  ARRAY -EFFECT  is  given  next. 

1  program  azTay.ellect 

2  c. . . 

3  c...  This  program  allova  lor  tha  comparison  of  the  analytic  and 

4  c...  discrata  Fouriar  transform  of  tha  valocity  rasponsa  of  a 
sc...  Timostianko  beam  to  tha  azparimantally  maasurad  forcing  function. 

6  c...  Also,  tha  azparimantally  obtained  valocity  response  may  be 

7  c...  plotted  as  vail.  Tha  number  of  sensors  used  to  measure  tha 

8  c...  theoretical  response  may  be  varied  as 

9  c...  vail  as  tha  measurement  frequency.  These  results  axe  plotted  as 

10  c...  a  function  of  vavenumber  at  a  puticular  frequency. 

11  real  x{64,4) .y(64,4) ,fftfor(64, 101) 

12  complex  fc(20,201) 

13  character  namee40 

14  common  frequency, irun,num_accel,ihan,m 

15  *>’  input  forcing  file' 

16  accept  100, name 

17  100  format (a40) 

18  type  a, 'number  of  shakers' 

19  accept  a,m 

30  type  *,  ’run  number' 

31  accept  * , irun 

22  if  (irun.le.3)  num_rd_fieq  =  201 

23  if  (irun.gt.3)  num.rd.freq  =  101 
34  c...  Reading  in  tha  forcing  function. 

2S  0PEH(UIITsl0,FILE*name,TYPE='OLD') 

DO  J  =  1,M 


36 


27 


28 


29 


30 

31 

32 

33 

34 

35 

36 

37 


IF(J.NE.{M+l)/2)  THEM 
DO  1=1  ,nuiB_rd_freq 

READ(10,133)  XXX. FA. PHASE 
FA  =  10** (FA/20.) 

FC(J.i)  =  FA*CMPLX(COSD(PHASE).SIHD(PHASE)) 

END  DO 

ELSEIF(J.Eq. (M+1)/2)THEN 
DO  1=1  .iiuin_rd_freq 

FCCJ.i)  =  CMPLXd. 0.0.0) 

ENDDO 

endil 


36 

39  C  .  .  . 

40  133 

41 

42  C.  .  . 

43  C  .  .  . 

44 

45 

46 

47 
46 

49 

50 

51  C.  .  . 

52  C .  .  . 

53  C .  .  . 

54  11000 

55 

56 

57 

58 

59 

60  C .  .  . 

61  C.  .  . 

62 

63 

64 

65 

66 

67  C .  .  . 
66  C.  .  . 

69  C .  .  . 

70 


END  DO 

fonnat(lx.f5.2.fl2.6.17.2) 

CLOSE (UNIT= 10) 

Construction  ths  nams  of  the  experinentally  measursd 
valocity  rasponsa  in  casa  this  option  is  dasirad 
if  (ffl.eq.9)nama  =’fft_'//nama(3:S)//’a’//naaa(8:8)//'n' 
if  (m.aq. 17)naDa  =*fft_ V/naaa(3:6)//’a'//naina(9:9)//'m' 
typa  a.naaa 

opan  (unit=19,forms'ttnforaattad* .status='old* .fila=nana) 

raad(19)fftfor 

closa(ttnit=19) 

call  proapt(idav)  !  plotting  initialization  subroutina. 

Which  fraquancy  and  than  dacida  in  which  fraquancy  bin  does 

tha  chosan  fraquancy  falls. 

typa  * , '  which  f raquancy  ' 

accapt  a. fraquancy 

fain  ~  2000. 

faax  s  12000. 

ibin  -  int(float(nuB_rd_fraq-l)a(froquancy-fmin)/10000.)al 
iibin=int ( 100 . a (f raquency-f ain) /lOOOO . ) al 

Call  tha  analytic  Fourier  transfora  of  tha  velocity  response, 
call  axftfbCfraquancy, fed, ibin)  .a.xCl.l)  .yCl.l)) 
do  i  =  1,64 

x(i,3)=x(i,l) 

y(i,3)afftfor(i,iibin) 

anddo 

Input  tha  nuaber  of  sensors  for  tha  theoretical  FFT  of  the 
velocity . 

typa  a,'  input  nuaber  of  accel' 
accapt  a,nuB_accal 


71 


ITS 


72  C  .  .  . 

73 


74 

75 

76 

77 

78  C  .  .  . 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93  C .  .  . 

94  C.  .  . 

95 

96  C .  .  . 

97  C  .  .  . 

98  C .  .  . 

99  C .  .  . 
100 

101 

102 

103 

104 

105 

106 

107 

108 
109 
no 
111 
112 

113 

114 

115 


116  C™ 


Call  the  theoretical  FFT  of  the  velocity, 
call  theo_llt(lrequency,ih2ui,lc(l  ,ibin)  ,nu]ii_accel.m,z(l  ,2) 
>  .y(i.2)) 

t3rpe  ■»,*  enter  numcurves* 
accept  4,nuB_curves 

call  zyplot(z,y.num_curves,idav)  !  plotting  routine, 
begin  again? 

type  another  frequency ( 1) ' 
accept  *,  iopt 
if  (iopt.eq.l)then 
do  i  =  1.64 

y(i.l)=0.0 
y(i.2)=0.0 
y(i.3)  =  0.0 

enddo 

goto  11000 
endif 

call  donepl(O) 

stop 

end 


subroutine  ezf tf b(f requeney >f c .m.hold . vout) 
this  program  reads  in  the  forcing  function  ezperimentally  measured 
and  uses  it  in  a  model  for  the  ezact  spatial  fourier  transform 
of  the  timoshenko  beam. 

REAL  E.G.RHO.IM, AREA. AK. DIE, TWO. TVEI.LM.Al.FSI 

REAL  PI.ZERO.BBT.BAL.FREq 

REAL  0KG(0:40),B(0:40}.AL(0:40).CETA 

REAL  BT(0:40) .DEL(0:40) .CHE(0:40) .LA)a>A(0:40) 

REAL  DF.HB.resp(64.101) 

REAL  ETAP(20}.BBTE.BALE,FRq(460) 

REAL  YR(0:40,20} 

REAL  XR,YI.fa.hold(64) 

REAL  0RT(0:40).PEHS(2).PITS(2) 

IITEGER  I.J.M.H.Z.ind.kik 

C0NPLEX*8  CY. SUM. DNOI (0:40) ,SSUN.f_temp(20) 

COMPLEXeS  FC(20) .SUMS, tempi, temp2,PHI{0:40) .c_i 
real  vout (64) 

CHARACTER4>20  BAME .  H 1 . 12 
CBARACTERel  P 
CBARACTERe21  lAMEP 


117 

118 

119 

120  C.  .  . 

121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134  C  .  .  . 

135 

136 

137 

138 

139  C .  .  . 

140  C.  .  . 

141  C.  .  . 

142 

143 

144 

145 

146 

147 

148 

149 

150  C.  .  . 

151  C.  .  . 

152  C.  .  . 

153 

154 

155 

156 

157 

158 

159 

160  133 

161 


DATA  ONE  /  1.  / 

DATA  ZERO. TWEH, TWO  /  0.0  .  12.  ,  2.  / 

PI=4*ATAK(0NE) 

constants  the  units  are  cgs  system, 
e  =  210el0  !  young’s  modulus 

g  =77  elO  !  shear  modulus 

rho  =7.8  !  density 

ak  =  .833  !  kappa  the  shape  factor 

df  =  .05  !  damping  factor 

num_eig=21  !  number  of  modes  used 

Im  =  37.9  !  length  of  the  beam 

sd  =  1 . 27  !  midth  of  the  beam 

ht  =  1.91  !  height  of  the  beam 

nl  =  'bsfreq*  !  eigenfrequency  file 

n2  =  'bsort*  !  orthoganality  constant  file 
c_i  =  cmplxCO. 0,1.0) 
wav_min  =0.0 

scaling  of  the  variables  back  to  inverse  centimeters, 
eav.maz  =  pielm 
area  =  vdeht 
ia  «  Bd*ht*e3/12. 

NB=(RB0*AREA) 

reading  in  the  eigen-frequencies 

OPEN  (  UNIT=10  ,  FlLE=»CgroBh.kgl.bar]'//nl  ,TyPE='0LD'  ) 
DO  I=2,nuffl_eig 
READ  (10,*)  ONG(I) 

0MG(I)=0MG(I)*2*PI 
END  DO 

CLOSE  (  UNIT* 10  ) 

0NG(0)=0.0 

0MG(i)=0.0 

reading  in  the  orthoganality  constants. 

OPEN (UIIT= 10 , FILE= • [grosh . kg 1 . bar] * //n2 , STATUS=  *  OLD  * ) 

DO  I=2,nua_eig 
READdO,*)  ORT(I) 

END  DO 

CL0SE(UIIT=10) 

0RT(0)=1.0 

0RT(1)=(LM**2)/12 

format(lx,fS.2,fl2.6,f7.2) 

CLOSE (UNIT* 10) 
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162  C.  .  . 

163  c...  location  oi  the  point  sources  (shakers)  depends  on  the  number  of 

164  c . . .  shakers 

165  il  (m.eq.9)delta=4. 

166 

167  if  (m.eq. 17)delta=2. 

168 

169  do  i  =  l,m 

170  etap(i)  =(2.95  +  (i-l)*delta  )/lm 

171  enddo 

172  c. . . 

173  c. . . 

174  c...  Begin  the  calculation  of  the  necessary  eigen  constants 

175  C.  .  . 

176  R=SQRT(IM/(AREA*LM**2)) 

177  S=SQRT(E*IM/(AK*AREA*G*IJI**2)) 

178  DO  Is2,nuffl_eig 

179  B(I)=SQRT(RHO*AREAeOMG(I)**2*LM**4/(E*IM)) 

180  AL(I) = ( ( (R**2-See2) **24-4/B(I) **2) *• . 5-R**2-S**2) . 5 

181  1  /SqRT(TWO) 

182  BT(I) =( ( (R**2-S**2) ••2+4/B(I) **2) ** . B+B**24-S**2) ** . B 

183  1  /SQRT(TWO) 

184  UMDA(I)-AL(I)/BT(I) 

185  CHE(1)=(BT(I)**2-S*e2)/(AL(1)**2+S**2) 

186  BAL=B(l)eAL(I) 

187  BBT=B(I)*BT(I) 

188  IF  (BAL.GT.TWEI)  GOTO  170 

189  DEL(I)=(COSH(BAL)-COS(BBT))/(LAMDA(I)*SIIH(BAL)- 

190  1  CHE(I) *311(881)) 

191  GOTO  171 

192  170  DEL(I)=1/LAMDA(I) 

193  171  corriiuE 

194  EID  DO 

195  C.  .  . 

196  C .  .  . 

197  c...  Calcualtlon  of  the  modal  expansion  coefficients,  a  .  -  tion 

198  c...  of  mode  shapes  and  the  forcing  function. 

199  c. . . 

200  DO  J=1,N 

201  YB(0.J)=1.0 

202  YR(1,J)=  .5e(ETAP(j)-.6) 

203  DO  I=2,num_eig 

204  BBT=B(I)*BT(I) 

205  BAL=B(I)*AL(I) 

BBTE*BBT*ETAP(J) 


206 
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207  BALE=BAL*ETAP(J) 

208  IF  (BALE.GT.TWEH)  GOTO  200 

209  YR(I,J)=(COSH(BALE)+COS(BBTE)/CHE(I)-LAMDA(I)*DEL(I) 

210  1  *SIHH(BALE)-DEL(I)*SIH(BBTE)) 

211  GOTO  201 

212  200  YR(I.J)=(COS(BBTE)/CHE(I)-DEL{I)*SIH(BBTE)) 

213  201  COITIHUE 

214  EHD  DO 

215  EHD  DO 

216  c...  at  each  fraquency  dnoi  holds  the  information  about  th« 

217  c...  orthoganality  of  th«  forcing  function  with  any  of  the  modes 

218  c . . .  at  each  frequency . 

219  FREQ  s  2.ePIefreqnency 

220  DO  I  =  0,num_eig 

221  SUM=CMPLX(0.0.0.0) 

222  DO  Jal.M 

223  SUM=  SUM  ♦  FC(J)*YR(I,J) 

224  EHD  DO 

225  TEMP1=1/(MB*0RT(I)*(0MG(I)**2*CMPLX(0HE,DF)-FREQ**2)) 

226  DHOI(I)  s  TEMPl*SUMeCMPLX(0.0,FREQ) !last  part  changes  to  velocity 

227  enddo 

228  c  . . 

230  C  CALCUUTE  THE  VAVEHUKBER  RESPOHSE  OF  EACH  NODE  SHAPE 

231  C  PHI  IS  THE  COMPLEX  ARRAY  OP  THE  VAVEHUMBER  RESPOHSE. 

233  c  do  loop  over  gammas - — 

234  c...  vavenumber  response  of  the  two  rigid  body  modes 

235  do  1010  kik  e  1,64 

236  hold(kik)  s  (vav_miii  +  («av_maz-vav_min)*(kik-l)/64.)/lm 

337  gamma  »  vav_ain  4-  (sav_maz-vav_ain)*(kik-l)/64. 

238  c . . . 

239  if  (  gamma. ne. 0.0)  then 

240  phi(O)  s  c_i/gammae{cezp(-c_iegamiaa)-l . ) 

241  phl(l)  -  (l./gammae*2)*(cexp(-c_iegamma)  -  1.) 

242  >  4-c_i/(2.egamma)e(cezp(-c_iegamma)4-l  ) 

243  elseif (gamma. eq. 0.0)  then 

244  phi(O)  *  CMPLXd.O.O.) 

245  phi(l)  »  CMPLX(.5,0.) 

246  endif 

247  C. .  THE  TRAHSFORN  OF  THE  lOH-RIGID  BODY  NODES. 

248  DO  Is2,num_eig 

249  BBT-B(I)*BT(I) 

250  BAL-B(I)eAL(I) 
call  transform(bal,bbt,lamda(i) ,del(i) ,che(i) ,1.0, gamma 


251 
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3S2  >  ,tamp2) 

253  call  trazislorm(taal,bbt,lainda(i)  ,del(i)  ,che(i)  ,0.0,gaiuna 

254  >  .tempi) 

255  phi(i)  =  temp2-templ 

256  enddo 

257  C - 

25S  SSUM=CMPLX(0.0.0.0) 

259  DO  Z=0,aam_eig 

260  SSUM=SSUM+DHOI(I)*PHI(I) 

261  EID  DO 

262  vottt(kik)  -  20elogl0(cabs(SSUM)) 

263  1010  continue  !  enddo  lor  uavenumber  loop 

264  1  00  00  lo]:mat(20a4) 

265  return 

266  end 

267  c. . . 

268  aubroutine  translormCbal.bbt, lambda, delta, zeta.xsi, gamma, veil) 

269  c...  Evaluation  ol  the  analytic  translorm  of  the  mode  shapes. 

270  real  lambda 

271  complez*8  val, tempi ,temp2,c_i 

272  c_i=cmplx(0. 0,1.0) 

273  bale  s  bal*zsi 

274  bbte  s  bbt*zsi 

275  gane  »  gamma^zsi 

276  c . . . 

277  tempi  =  cezp(-c_i'*game)/(gaiiimaee2  ♦  bal**2) 

278  tempi  -  templ*((c_i*gaiiima*co8h(bale)  4-  beil*sinh(bale)) 

279  >  -lambda«delta«(c_iegamma*sinh(bale)  4-  bailecosb(beile))) 

280  c . . .  il  argument  ol  the  hyberbolic  lunctions  is  too  high 

281  il(bale.gt.l2)  tempi  -  0.0 

282  il  (gamma  .  ne  .  bbt)  then 

283  c . . . 

284  temp2  =  cezp(-c_i*game)/(gamma*a2  -  bbt**2) 

285  temp2  =  temp2*(l./zeta*(c_i*gamma*cos(bbte)  -  bbtesinCbbte)) 

286  >  -deltae(c_i*gamma*8in(bbte)  4-  bbt*cos(bbte))) 

287  C .  . 

288  elseil (gamma  .  eq  .  bbt)  than 

289  temp2  =  l./zetae(z8i  -  l./(2.*c_i*gamma)*cezp(-2.*c_i*game)) 

290  >  -delta*(z8i  ♦  l/(2.*c_i*gamma)*cezp(-2.ec_i*game)) 

291  endil 

292  val  s  tempi4temp2 

293  return 

294  end 

295  c . . . 

396  subroutine  theo_llt(lrequancy.ihan.lc.num.accel,m.hold,vout) 
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297  c...  calc  1ft  ol  respones*  theoretically 

298  c...  First  the  spatial  response  of  the  bean  is  obatained  at  num.accel 

299  c . . .  points .  Then  the  FFT  is  taken  of  the  respose 

300  C — 

301  C - 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321  C.  . 

322  C .  . 

323  C .  . 

324  c...  The  sane  physical  constants  are  used  in  this  routine  as  in  EXFTFB 

325  c . . . 

326  C .  .  . 

327  DATA  E.G.OHE.REO.LN.sd.ht  /  210E10  ,  77E10  . 

328  1  1.  .  7.8  ,  37.9  ,1.27,1.91/ 

329  DATA  ZERQ,AK,DF,TVEI,TVa  /  0.0  ,  .833  ,  .05  ,  12.  ,  2.  / 

330  PI«4*ATAI(0IE) 

331  PElS(l)sl.O 

332  PElS(2)-2.0 

333  area  ■  «d*ht 

334  i*  *  ad*ht**3/12. 

335  nun_fftsl28 

336  NB-fRH0*AREA«LM) 

337  P* '  \  ’ 

338  nua_eigs21 

339  c . .  before  things  get  started  calculate  f ft  weights 

340  call  cafft(a,ww,nua_fft,0)  !  FFT  alogorithm 
banning  window  option. 


REAL  E. G,RH0,IM, AREA, AK.OHE, TWO. TWEH.LM.Al, FSB 
REAL  PI, ZERO, BBT, BAL, FREQ 

REAL  0MG(0;40) ,PHI(0:40) .B(0:40) .AL{0:40) ,CETA 
REAL  BT(0:40) ,DEL(0:40) .CHE(0:40) ,LAMDA(0:40) 
REAL  DF,MB 

REAL  ETAP(20),BBTE,BALE,FRq(450) 

REAL  YR(0:40.20),x(64) 

REAL  XR , YI . PHS ( 500 ) , f a . vout ( 64 ) . hold( 64 ) 

REAL  0RT(0;40),PEHS(2).PHTS(2) 

IHTEGER  I,J.M.N,Z.ind 

COMPLEXES  CY,SUM,DH0I(0:40) .SSUM. tempi 

COMPLEXES  FC(20},SUMS,a(1024),ww(512},resp(303) 

CEARACTERe20  BANE, 11. 12 

CHARACTER*!  P.cone 

CHARACTER*2  ctwo 

CHARACTER«21  lAMEP 

charaeter*6  base 


341  c. . . 
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342  ^7P*  *'•'  banning  optionCl)' 

343  accept  * , ihan 

344  c . . . 

345  c...  option  to  use  spacing  exactly  the  same  as  per  the  experiments 

346  c...  or  a  spacing  that  evenly  spreads  the  accelerometers  over  the 

347  c. . .  beam  ior  any  spacing. 

348  c . . . 

349  141  type  same  as  exp.(l)  or  even  dist.(2)’ 

350  accept  e.imatch 

351  if (imatch.ne.l.and.imatch.ne.2)goto  141 

352  c 

353  c  reading  in  the  eigen-frequencies  from  the  file 

354  c  BSFREQ.DAT 

355  c 

356  OPEI  (  UIIT=10  .  FILE=' Cgrosh.kgl.barlbsfreq.DAT'  .TYPE='0LD'  ) 

357  DQ  l=2.num_eig 

358  READ  (10,*)  QMGd) 

359  0MG(I)*0NG(I}*2*PI 

360  EHD  DO 

361  CLOSE  (  URIT*10  ) 

362  0NG(0)«0.0 

363  0MG(1)«0.0 

364  c . .  . 

365  c . . .  reading  in  the  orthoganality  conatants . 

366  c. . . 

367  OPEI (UIIT*10 . FILE* ' Cgrosh .kgl . bar] bsort . dat ' , STATUS* ’ OLD ' ) 

368  DO  I*2,num_eig 

369  READClO,*}  ORT(I) 

370  EID  DO 

371  CL0SE(UIIT*10) 

372  0RT(0)al.0 

373  0RT(l)a(LM**2)/12 

374  1222  format(lx,f6,2,fl2.6,f7.2) 

375  CL0SE(UIIT*10) 

376  c . . .  Eigenf nnciton  constants 

377  c . . . 

378  R*SQRT(IN/(AREA*LM**2)) 

379  S«SQRT(E*IN/(AK*AREA*C*LN**2)) 

380  DO  I*2,nnm_eig 

381  B(I)  aSqRT(RHO*AREA*OMG(I)  **2*LI!**4/  (E*H!)  ) 

382  AL(I)»(((R**2-S**2)**2+4/B(I)**2)**.6-R**2-S**2)**.5 

383  1  /SqRT(TWO) 

384  BT(I) * ( ( (R**2-S**2) **2*4/B(I) **2) •* . 64-R**2+S**2) ** . 5 

385  1  /SqRT(TWO) 

386  LAI(DA(I)*AL(1)/BT(I) 


387  CHE(I)=(BT(I)**2-S**2)/(AL(I)**2+S**2) 

388  BAL=B(I)*AL(I) 

389  BBT=B(I)*BT(I) 

390  IF  (BAL.GT.TWEH)  GOTO  170 

391  DEL(I)=(COSH(BAL)-COS{BBT))/(LAMDA(I)*SIHH(BAL)- 

392  1  CHE(I)*SIH(BBT)) 

393  GOTO  171 

394  170  DEL(I)=1/LAMDA(I) 

395  171  COHTIHUE 

396  END  DO 

397  C - 

398  C  CALCULATIONS  DUE  TO  THE  FORCING  FUNCTION.  EVALUATING  THE  POINT  FORCES 

399  C  AT  THE  DIFFERENT  POINTS 

400  C 

401  if  (m.«q.9)  delta  >  4. 

403  if  (m.eq.l7)  delta  =  2. 

403  do  i  =  l.m 

404  etap(i)  =(2.9S  4-  (i-l)*delta)/lin 

405  enddo 

407  DO  J=1,M 

408  YR(0,J)*1.0 

409  YR(l,J)=.6*(ETAP(j)-.5) 

410  DO  Is2,xiua_eig 

411  BBT=B(I)eBT(I) 

413  BAL=B(I)eAL(I) 

413  BBTE=BBTeETAP(J) 

414  BALE^BALeETAPCJ) 

415  IF  (BALE.GT.TVEN)  GOTO  200 

416  YR(I,J)*(C0SH(BALE)4-C0S(BBTE)/CHE(I)-LAMDA{I)*DEL(I) 

417  1  *SIIH(BALE)-DEL(I)*SIN(BBTE)) 

416  GOTO  201 

419  200  YRCl.J)=(COS(BBTE)/CHE(I)-DEL(I)eSIN(BBTE)) 

420  201  CONTINUE 

431  END  DO 

422  END  DO 

423  c...  at  each  frequency  dnoi  holds  the  information  about  the 

424  c...  orthoganality  of  the  forcing  function  vith  any  of  the  modes 

425  c . . .  at  each  frequency . 

426  FREQ  *  2.*PI*frequeacy 

427  DO  I  a  0,num_eig 

436  SUN=CMPLZ(0.0.0.0) 

429  DO  J>1,N 

430  SUM«  SOM  +  FC(J)*YR(I,J) 

END  DO 


431 


432 

433 

434 

435  C  .  .  . 

436  C  .  .  . 

437  C  .  .  . 

438  C  .  .  . 

TEMP1=1/(MB*0RT(I)*(0MG(I)**2*CMPLX(0NE,DF)-FREQ4-4.2)) 
DNOI(I)  =  TEMP1*SUM*CMPLX(0.0. FREQ) Hast  part  changes  to 
enddo 

At  this  point  all  of  the  information  concerning  the 
forcing  function  has  been  obtained.  Hov,  the  sensor 
location  is  varied. 

1S6 

velocity 

440  10 

FORMAT  (A20) 

441  C.  .  . 

442  C  .  .  . 

443  C  .  .  . 

begin  acclerometer  loop 

444  C .  . 

445 

do  kkj  =  l,nun_accel 

446  1000 

foxmat(il) 

447  2000 

format (i2) 

448  C  .  .  . 

449  C  .  .  . 

Bhat  is  the  location  of  the  accel.  for  this  kkj . 

450 

if (nuffl_accel.eq.17)  then 

451 

r.off=2.96/lm 

452 

ceta  =  (2.96+float(kXj-l)*2. )/lm 

453 

delta.fft  =  2. 

454 

elseif (nuffl_aceel.eq.9)  then 

455 

x_off*2.9S/lB 

456 

ceta  =  (2.96+float(kkj-l)*4.)/la 

457 

delta.x  =  4. 

458 

delta.fft  =  4. 

459 

else 

460 

if ( imat  ch.eq.l)then 

461 

x_off=0.95/lm 

462 

delta_x=  (Iffl-i . 9 ) /float (num.accel- 1 ) 

463 

ceta  =x_ofl+ (float (kX j -1 ) ) *delta_x/lm 

464 

delta_fft=delta_x 

465 

elseif  (iDatch.eq.2  )then 

466 

ceta  s  (float(kkj)-.E)/float(nuffl_accel) 

467 

delta.f f t=lffl/f loat (num.accel ) 

468 

endif 

469 

endif 

470  C  .  .  . 

471  C - 

472  C 

473  C  .  .  . 

evaluate  the  mode  shapes  at  ceta. 

474  C  .  .  . 

475 

DO  ls2,nuffl_eig 

476 

BBT=B(I)*BT(I) 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487  C 

488  C 

489  C 

490 

491 

492  C- 

493 

494 

495 

496 

497  C. 

498  C. 

499 

500 

501 

502 

503 

504  C. 

505  C. 

506  C. 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

530 

531 


BAL=B(I)*AL(I) 

BALE=BAL*CETA 
BBTE=BBT*CETA 
IF  (BALE.GT.TWEN)  GOTO  150 

PHI(I)=C0SH(BALE)-LAMDA(I)*DEL(I)4SI1IH(BALE)+ 
1  C0S(BBTE)/CHE(I)-DEL(I)4.SIK(BBTE) 

GOTO  151 

150  PHr(I)=COS(BBTE)/CHE(I)-DEL(I)*SIH(BBTE) 

151  COHTIHUE 
END  DO 


THE  RIGID  BODY  MODES.  BOTH  ROTARY  AND  TRANSLATIONAL 

PHI(0)=1.0 

PHI(l)=.6*(CETA-.5) 


SStIM=CMPLX(0. 0,0.0) 

DO  l=0.tium_«ig 
SSUM=SSUM+DNOI(I)*PHI(I) 
END  DO  !  enddo  for  kkj 


typ«  j j j ,c«ta,c«ta*l0,cab8(8SUD) 

aritedO,  *) j j j ,c«ta.cab8(88U]n) 

Cys881UI 
r«8p(kkj}  =  cy 
enddo  I  enddo  Xkj 

. .  Next,  perform  FFT  on  the  velocity  data  and  return  to  main 
. .  routine. 

conat  =  2 . *pi*delta_lf t 

z_norm  =  lm/float(num_accel)  !  proper  nozmaJ-ization 
do  j  -  l.num_accel 

if  (iban.eq.l)  then 

weights. 5*(l.-co8{(x_off+float(j-l))*const)) 

else 

weights 1. 
endif 

a(j)  s  eeighterespCj) 

enddo 

do  j  s  num_accel4l,num_fft 
a(j)scmplz(0. 0,0.0) 

enddo 

call  C5ifft(a,ww,nuffl_fft,l) 
do  ik  «  1,64 


522 


> 


voutCik)  =  20.*alogl0(cab8(a(ik)*x_nonn)) 
hold(ik)  =float(ik-l)*pi/(lloat(niim_llt/2) 

•delta_fft) 
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523 

524 

525 

526 

527 
52S 


enddo 

return 

end 


Finally,  the  FORTRAN  program  that  compares  the  finite  and  infinite  bezim 
admittances  at  a  particular  wavenumber,  titled  COMPARE,  is  listed. 


1 


2 

C  .  .  . 

3 

C  .  .  . 

4 

C.  .  . 

5 

C.  .  . 

6 

7 

6 

C 

9 

10 

U 

12 

13 

14 

c 

IS 

16 

17 

18 

19 

20 

211 

21 

112 

22 

190 

23 

24 

25 

c. . . 

26 

c .  .  . 

27 

28 

29 

30 

31 

1804 

32 

1806 

33 

10 

program  compare 

This  program  compares  the  predicted,  experimentally  measured 
finite  beam  admittance  to  the  theoretically  derived  infinite  beam 
admittance . 

CHARACTER  FHAME*30 
LOGICAL*!  TEXT(80) 

real  *8  gamma 

DIMEISIOH  X(101),q(2),qy(2) 
real  mag(101,4) 

CHARACTER  FILEII*40 .texpl*40,header«40 
CHARACTER*30  LIHE_TyPE(4) ,cone*l ,ctuo*2 

pi  =  3.141592654 

^7P*  * • '  input  number  of  shakers  ' 

accept  211,ctvo 

type  input  ezper  no  in  the  series  (ie,<=5)  ' 

accept  112, cone 

format  (a2) 

format (al) 

format(a40) 

ncurves  =  3 

npts  =  101 

variables  for  the  plotting  routine 
theses  are  the  line  descriptions 
line.typed)  ^  'exact  transform!' 
line_type(4)  *  'fft  of  theoretical  results!' 
line_type(2)3'fft  of  experimental  results!' 
line_type(3)* 'admittance  of  infinite  beam!' 
formate  input  filename  ',il,'  — >',!) 

formate  input  LIIETITLE' ,il. '  -->',!) 

format (a30) 
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34  C  .  .  . 

35  C  .  .  . 

36  C .  .  . 

37 
3d 

39 

40 

41 

42 

43  9  99 

44  10000 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 

62  C.  .  . 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 


begin  frequency  loop 

do  i  *  1,101 
x(i)  =  2.  +lloat(i-l)*.l 
enddo 

CALL  PRQMPT(IDEV)  !  plotting  initialization  subroutine 
type  fill  letters  (1)* 
accept  999,ifil 
fonnat(il) 

type  <*.'  automatic  option  or  self  input  (1,2)’ 
accept  999,iopt 
if  (iopt.eq.l)  then 
if  (ctuo .eq. '9* )  then 
if  (cone.eq. '1' .or. cone. eq. ’5’)then 
ibin  *  1 
gammas . 0 

elseif  (cone.eq. ’2')then 
ibin  s  16 
gamma  =  .73631 
elseif  (cone.eq. '4*)then 
ibin  s  4 
gammas . 14726 

elseif  (cone.eq. ’3’)then 
ibin  s  8 
gamma  =  .34361 

endif 

endif 

if  (ctvo.eq. ’17’}  then 

if  (cone.eq. ’1 ’)tbea 
ibin  s  1 
gammas . 0 

elseif  (cone.eq. ’4’ }tben 
ibin  s  16 
gasima  =  .73631 
elseif  (cone.eq. ’2’)tben 
ibin  s  4 
gammas . 14726 
elseif  (cone.eq. ’3’)tben 
ibin  *  8 
gamma  *  . 34361 
elseif (cone.eq. ’6’)  then 
ibin  s  31 
gamma  «  1.4726 


78 


79 


«ndif 

endif 
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ao 

61 

62 

63 

64 

65 

66 

67  C  .  .  . 
66  C.  .  . 
69  C  .  .  . 

90 

91  C.  .  . 
93  C  .  .  . 

93 

94 

95 

96  C 

97  C  .  .  . 
96 

99 

100 

101 

103 

103 

104 

105 

106  C 

107  60 
106  C  .  .  . 

109  C.  .  . 

110  C.  .  . 

111  C.  .  . 
113  C.  .  . 

113  C.  .  . 

114  C.  .  . 

115  C.  .  . 

116  C.  .  . 

117  C.  .  . 
116  C.  .  . 


119 


130 


131 


133 


elseil  (iopt.eq.2)then 

*>'  input  bin  number  (calculates  gamma)’ 
accept  *,ibin 

gamma  =  pi*(ibin-l)/64.*l .dO 

endil 

read  in  calculated  and  measured  velocity  responses, 
these  values  are  stored  in  a  data  file 

call  indat (mag, cone. ctuo.ibin) 

program  that  calculates  infinite  beeun  admittance 
call  inl.admit(mag(l,3) .x.gamma) 
zmin  =  z(l) 
zmax=  x(lOl) 

scaling  lor  plot 
header=cone// '  ’//ctBo//'$» 

il  (ibin.le.20)  then 
ymin«-160 
ymaxa-60 

elseil  (ibin.gt.20)  then 
ymins-200 
ymax=-100 

endil 

continue 


plotting  routine  resides  in  this  space 


do  you  vish  to  continue  questions, 

type  TYPE  1  lor  another  plot,  0  to  ask  lor  another  run’ 
accept  999,  iplag 
il(iplag.eq. 1)  go  to  60 

type  TYPE  1  lor  another  run,  0  to  exit’ 
accept  999,  irllag 
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124  if  (irllag.eq.l)  goto  10000 

125  44  CALL  DONEPL(O) 

126  Stop 

127  end 

128  C.  .  . 

129  C.  .  . 

130  C  .  .  . 

131  subroutine  ini_adinit(xinag,lrq,ganiiia) 

132  c. . . 

133  c...  this  subroutine  calculates  the  savenuinber-lrequency 

134  c...  admittance  ol  an  infinite  beam. 


135  c . . . 

136  implicit  realms  (a-h.o-z) 

137  real*8  rho.im, area, ak, one, two, toen.lm 

138  realms  pi, zero, bbt.bal, Iraq 

139  realms  dl 

140  re2J.*4  xmagClOl) ,lrq(10l) 

141  re2il*8  first 

142  real*8  kappa 

143  complex*16  resp,argl,arg2,g,e,bomega,zeta, theta 

144  complex  *16  wav_free_4 

145  integer  i, j ,m,n,z,ier,IWEIGH 

146  cccccccccccccc  Input  material  and  dimension  parameters  ccccccccccccccc 

147  data  real_e,real_g,one,rho  /  210el0  ,  77el0  , 

148  1  1 .  ,  7 . 8  / 

149  data  zero, ak, df ,twen, two  /  0.0  ,  .833  ,  .08  ,  12.  ,  2.  / 

150  pi=4.*datan(one) 

151  212  format  (al5) 

152  ht=1.92d0 

153  wd=  1.27d0 

154  area  =  ht  *  wd 

155  kappa  =  ak 

156  im  *  wd*ht**3/12. 

157  g  =  real_g*dcmplx(1.0d0,dl) 

158  e  *  real_e*dcmplx(1.0d0,df ) 

159  write(6,*)g,e,ak,df 

160  C...  THIS  PART  OF  THE  PROGRAM  CALCULATES  THE  wavenumber 

161  C...  admittance  AT  VARIOUS  frequencies 
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163 

164 

165 

166 

167 

168 


zeta  =  kappa  *  g  *  area 

DO  jj  «=  l.lOl 
freq  =  frq(j j)*1000. 
hold_max=- 1000 . 

bomega  =  rho  *  area  •  (2.*pi  *freq)**2 
theta  =  rho  *  im  •  (2.*pi  •freq)**2 
wav_free_4  ~  bomega  /e/im 
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169 

170 

171 

172 

173 

174 

175 

176 

177 


argl  =  ga]iuiia«*2/zeta-*-l/e/iin*(l-theta/zeta) 
arg2  =  gaiUBa**4  +  g2unBa**2  •  (-th«ta/«/ijn  - 
>  bomaga/zata)  +  ®av_frae_4*(thata/zata  -  1) 

raap  =  -argl/arg2 

nnag  =  20adlogl0(  cdabs(rasp)alraqa2.api) 
phasa=datan2d(diinag(rasp)  ,  draal(rasp)} 
xinag(jj)  =  sngKnnag) 
anddo 
raturn 
EHD 
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APPENDIX  D 

EXPERIMENTAL  EXCEPTIONS 

In  this  appendix,  any  problems  associated  with  collecting  the  data  amd  the 
remedies  are  delineated.  In  some  of  the  experiments,  accelerometer  data  were 
lost.  When  FFT’s  of  the  experiments  were  performed,  these  data  locations  were 
taken  to  be  zero.  Obviously,  this  induces  error  into  the  measured  wavenumber 
velocity  response.  The  accelerometer  data  at  location  one  were  missed  in 
Experiments  1  and  4.  Additionally,  locations  34  and  37  were  lost  in  Experiment  4. 
There  was  one  problem  associated  with  the  force  gage  records.  In  Experiment  6, 
the  last  force  record  was  lost.  After  examining  the  force  records  at  the  other 
locations,  a  strong  symmetry  was  noticed  (fortunately).  Using  this  symmetry, 
an  approximation  was  made  which  was  to  use  the  first  shaker  force  record  for  the 
seventeenth  shzdcer.  As  can  be  seen  by  the  excellent  agreement  between  theory 
and  experiment  for  the  Experiment  6,  this  was  a  good  approximation. 

Also,  a  note  should  be  made  concerning  the  storage  of  the  data  on  the 
computer.  For  the  second  and  fourth  experiments,  the  data  are  stored  in  different 
subdirectories  on  the  computer  tham  the  experiment  number  listed  in  this  thesis. 
Experiment  4  is  stored  in  the  second  subdirectory  and  Experiment  2  is  stored 
in  the  fourth  subdirectory.  Also,  experimentadly  there  is  a  difference  in  the 
number  of  frequencies  measured  from  the  first  three  experiments  and  the  last 
six  experiments.  In  the  first  three  experiments  the  data  was  taken  using  50  Hz 
increments  while  the  final  six  were  taiken  using  100  Hz  increment.  The  range  of 
frequencies  measured  was  same  for  all  experiments  (2  kHz  to  12  kHz). 


194 

Finally,  it  should  be  noted  that  the  180  degree  phase  factor  present  in  the 
force  gage  calibrations  was  not  subtracted  out  in  the  data  acquisition  programs. 
This  phase  factor  is  removed  by  any  FORTRAN  program  that  made  use  of  the 
measured  admittzmce  data. 
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